Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,195,870 T→A 20.4% I316I (ATT→ATA aceB → malate synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,195,8700T→A20.4% 2.5 ‑3.1 0.041948aceBmalate synthase
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (11/20):  total (11/30)
Fisher's exact test strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGA  >  REL606/4195835‑4195901
                                   |                               
aaaCCATTCCTGAATGCTTACTCACGCCTGTTGAtt                                 <  1:2956272/36‑1 (MQ=255)
aaaCCATTCCTGAATGCTTACTCACGCCTGTTGAtt                                 <  1:2534825/36‑1 (MQ=255)
aaaCCATTCCTGAATGCTTACTCACGCCTGTTGAtt                                 <  1:2519296/36‑1 (MQ=255)
aaaCCATTCCTGAATGCTTACTCACGCCTGTTGAtt                                 <  1:14919/36‑1 (MQ=255)
aaaCCATTCCTGAATGCTTACTCACGCCTGTTGAtt                                 <  1:4186189/36‑1 (MQ=255)
 tgccATGNNNGAAGGCTGACTCACGCCTGTTGATTa                                <  1:2026811/34‑1 (MQ=16)
 aaCCATTCCTGAATGCTTACTCACGCCTGTTGATTa                                <  1:3457093/36‑1 (MQ=255)
   ccATTCCTGAATGCTTCCTCACGCCTGTTGATTaaa                              <  1:4005577/36‑1 (MQ=255)
   ccATTCCTGAATGCTTACTCACGCCTGTTGATTaaa                              <  1:3062250/36‑1 (MQ=255)
    cATTCCTTAATGCTTACTCACCCCTGTTGATTaaaa                             <  1:4265783/36‑1 (MQ=37)
     aTTCCTGAATGCTTACTCGCGCCTGTTGATTaaacc                            <  1:1519203/36‑3 (MQ=37)
      ttCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                           >  1:2933909/1‑36 (MQ=255)
      ttCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc                           <  1:2909851/36‑1 (MQ=255)
      ggccTGCNNNCTGACTCCCGCCTGTTGATTAAAAcc                           <  1:2002104/34‑1 (MQ=12)
       tCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt                          >  1:1182692/1‑36 (MQ=255)
         cTGAATGCTTACTCACGCCTGTTGATTAAACCCTGc                        <  1:3685130/36‑1 (MQ=255)
         cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc                        <  1:4578997/36‑1 (MQ=255)
          tGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc                       <  1:3912371/36‑1 (MQ=255)
          tGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc                       >  1:5405071/1‑36 (MQ=255)
           gAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa                      >  1:272063/1‑36 (MQ=255)
           gAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa                      >  1:1228075/1‑36 (MQ=255)
           gAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa                      <  1:3263624/36‑1 (MQ=255)
           gAATGCTTACTCACGAGTGTTGATTAAAACCTGCCa                      >  1:5562134/1‑36 (MQ=38)
            aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCAt                     <  1:2536679/36‑1 (MQ=255)
              tGCTTCCTCACGCCTGTTGATTAAAACCTGCCATaa                   <  1:5721354/36‑1 (MQ=255)
              tGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa                   <  1:143481/36‑1 (MQ=255)
                 ttACTCACGCCTGTTGATAAAAACCTGCCATAAAcg                <  1:3908870/36‑1 (MQ=255)
                  tACCCACGCCTGTTAATTAAACCCTGCCATAAAcgc               <  1:5167951/36‑1 (MQ=25)
                   acccACGCCTGTTGATTAAAACCTGCCATAAAcgcg              <  1:5581834/33‑1 (MQ=255)
                   aCTCACGCCTGTTGATAAAACCCTGCCATAAAcgcg              <  1:155578/36‑1 (MQ=37)
                     tCACCCCTGTTGATAAAACCCTGCCATAAACGCGGt            <  1:1686676/36‑1 (MQ=25)
                     ccACGCCTGTTAATAAAAACCTGCCATAAACGCGGt            <  1:5162945/35‑1 (MQ=37)
                      cACGCCTGTTGATTAAAACCTGCCATAAACGCGGTg           >  1:2229391/1‑36 (MQ=255)
                       aCGCCTGTTGATTAAAACCTGCCATAAACGCGGTGc          <  1:5424288/36‑1 (MQ=255)
                        cGCCTGTTGATTAAAACCTGCCATAAACGAGGCGCt         >  1:830000/1‑36 (MQ=37)
                         gCCTGTTGATTAAACCCTGCCATAAACGCGGTGCtt        <  1:4582656/36‑1 (MQ=255)
                         gCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtt        >  1:4865213/1‑36 (MQ=255)
                         gCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtt        >  1:3496758/1‑36 (MQ=255)
                         cccTATTGATAGAACCCTGCCATAAACGCGGTGCtt        <  1:5745114/35‑1 (MQ=16)
                           ttGTGGATAAAACCCTGCCATAAACGCGGTGCtttt      <  1:4186857/35‑1 (MQ=25)
                           cttttGATTAAACCCTGCCATAAACGCGGTGCtttt      <  1:607867/33‑1 (MQ=37)
                           cTGTTAATAAAACCCTGCCATAAACGCGGTGCtttt      <  1:2573307/36‑1 (MQ=25)
                            tattGATCAAACCCTGCCATAAACGCGGTGCTTTTg     <  1:5471301/34‑1 (MQ=37)
                            tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTg     >  1:5188904/1‑36 (MQ=255)
                             ggggATAAAACCCTGCCATAAACGCGGTGCTTTTGc    <  1:1021753/33‑1 (MQ=25)
                               tGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGa  <  1:742173/36‑1 (MQ=255)
                                   |                               
AAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGA  >  REL606/4195835‑4195901

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG