Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,780,211 A→T 7.7% intergenic (‑374/‑148) thrS ← / → arpB threonyl‑tRNA synthetase/hypothetical protein; ankyrin‑repeat protein B; b1720(b4494)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,780,2110A→T 7.7% 2.3 ‑1.4 ‑0.039intergenic (‑374/‑148)thrS/arpBthreonyl‑tRNA synthetase/hypothetical protein; ankyrin‑repeat protein B; b1720(b4494)
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (19/17):  total (22/17)
Fisher's exact test strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.99e-01

GCAAAGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAACGCTCATGAGATTAATGACGAAGTGGTCATATC  >  REL606/1780176‑1780246
                                   |                                   
gCAAAGTTCGTTGCGTGATTTGGCGGTAGTATTTTa                                     <  1:1476666/36‑1 (MQ=255)
gCAAAGTTCGTTGCGTGATTTGGCGGTAGTATTTTa                                     <  1:373923/36‑1 (MQ=255)
 cAAAGTTCGTTGCGTGATTTGGCGGTAGTATTTTaa                                    <  1:3240833/36‑1 (MQ=255)
   aaGTTCGTTGCGTGATTTGGCGGTAGTATTTTtaaa                                  >  1:1180730/1‑35 (MQ=255)
   aaGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAAc                                  <  1:5542225/36‑1 (MQ=255)
   aaGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAAc                                  <  1:1599166/36‑1 (MQ=255)
   aaGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAAc                                  <  1:7507496/36‑1 (MQ=255)
    aGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAACg                                 <  1:6011862/36‑1 (MQ=255)
      ttCGTTGCGTGATTTGGCGGTAGTATTTTTAACGCt                               >  1:6000121/1‑36 (MQ=255)
      ttCGTTGCGTGATTTGGCGGTAGTATTTTAAACGCt                               <  1:1496799/36‑1 (MQ=255)
      ttCGTTGCGTGATTTGGCGGTAGTATTTTAAACGCt                               <  1:6728430/36‑1 (MQ=255)
      ttCGTTGCGTGATTTGGCGGTAGGATTTTTAAAGCt                               >  1:4696180/1‑36 (MQ=25)
       tCGTTGCGTGATTTGGCGGTAGTATTTTAAACGCTc                              <  1:7503574/36‑1 (MQ=255)
        cGTTGCGTGATTTGGCGGTAGTATTTTAAACGCTCa                             <  1:5060875/36‑1 (MQ=255)
           tGCGTGATTTGGCGGTAGTATTTTAAACGCTCATga                          <  1:6688257/36‑1 (MQ=255)
            gCGTGATTTGGCGGTAGTATTTTAAACGCTCAttat                         >  1:436072/1‑33 (MQ=255)
            gCGTGATTTGGCGGTAGTATTTTAAACGCTCATgag                         <  1:2322458/36‑1 (MQ=255)
             cGTGATTTGGCGGTAGTATTTTAAAACCTCATgaga                        >  1:7787928/1‑36 (MQ=38)
               tGATTTGGCGGTAGTATTTTAAACGCTCATGAGAtt                      <  1:4760621/36‑1 (MQ=255)
                gATTTGGCGGTAGTATTTTAAACGCTCATGAGATTa                     <  1:8208047/36‑1 (MQ=255)
                 aTTTGGCGGTAGTATTTTAAACGCTCATTAGAttta                    >  1:7465093/1‑34 (MQ=38)
                 aTTTGGCGGTAGTATTTTAAACGCTCATGAGATTaa                    >  1:8335440/1‑36 (MQ=255)
                 aTTTGGCGGTAGTATTTTAAACGCTCATGAGATTaa                    <  1:8110375/36‑1 (MQ=255)
                  tttGGCGGTAGTATTTTTTAAACTCATTAGAttttt                   >  1:6683860/1‑33 (MQ=9)
                     ggCGGTAGTATTTTAAACGCTCATTAGATTAATTAc                >  1:8006006/1‑36 (MQ=38)
                     ggCGGTAGTATTTTAAACGCTCATGAGATTAATGAc                <  1:6840934/36‑1 (MQ=255)
                       cGGTAGTATTTTAAACGCTCATGAGATTAATGACGa              <  1:7119908/36‑1 (MQ=255)
                       cGGTAGTATTTTAAACGCTCATGAGATTAATGACGa              >  1:7873818/1‑36 (MQ=255)
                        ggTAGTATTTTAAACGCTCATGAGATTAATGACGaa             >  1:113982/1‑36 (MQ=255)
                        ggTAGTATTTTAAACGCTCATGAGATTAATGACGaa             >  1:2794627/1‑36 (MQ=255)
                          tagtaTTTTAAACGCTCATGAGATTAATGACGAAGt           >  1:865386/1‑36 (MQ=255)
                            gtaTTTTAAACGCTCATGAGATTAATGACGAAGTgg         >  1:3051712/1‑36 (MQ=255)
                             taTTTTAAACGCTCATGAGATTAATTACGAAGTGGt        >  1:5732155/1‑36 (MQ=255)
                             taTTTTAAACGCTCATGAGATTAATGACGAAGTGGt        >  1:4521538/1‑36 (MQ=255)
                             taTTTTAAACGCTCATGAGATTAATGACGAAGTGGt        >  1:142124/1‑36 (MQ=255)
                             taTTTTAAACGCTCATGAGATTAATGACGAAGTGGt        >  1:600342/1‑36 (MQ=255)
                               ttttAAACGCTCATGAGATTAATGACGAAGTGGTCa      >  1:292023/1‑36 (MQ=255)
                                tttAAACGCTCATGAGATTAATGACGAAGTGGTCat     <  1:6894942/36‑1 (MQ=255)
                                tttAAACGCTCATGAGATTAATGACGAAGTGGTCat     <  1:409542/36‑1 (MQ=255)
                                tttAAACGCTCATGAGATTAATGACGAAGTGGTCat     >  1:3995524/1‑36 (MQ=255)
                                tttAAACGCTCATGAGATTAATGACGAAGTGGTCat     >  1:205433/1‑36 (MQ=255)
                                 ttAAACGCTCATGAGATTAATGACGAAGTGGTCata    >  1:4267963/1‑36 (MQ=255)
                                  tAAACGCTCATGAGATTAATGACGAAGTGGTCatat   >  1:6772189/1‑36 (MQ=255)
                                   aaaCGCTCATGAGATTAATGACGAAGTGGTCATATc  <  1:2916020/36‑1 (MQ=255)
                                   aaaCGCTCATGAGATTAATGACGAAGTGGTCATATc  <  1:2894573/36‑1 (MQ=255)
                                   |                                   
GCAAAGTTCGTTGCGTGATTTGGCGGTAGTATTTTAAACGCTCATGAGATTAATGACGAAGTGGTCATATC  >  REL606/1780176‑1780246

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG