Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,0170AT63.9% 78.0 / 86.2 72I125I (ATT→ATA) ‡nohANohA
Reads supporting (aligned to +/- strand):  ref base A (7/19);  new base T (24/22);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 4.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGC  >  NC_000913/1633925‑1634112
                                                                                            |                                                                                               
attatatCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAAtt                                                                                              <  2:3253644/93‑1 (MQ=1)
            ttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATc                                                                                  >  1:9440987/4‑96 (MQ=0)
            ttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATc                                                                                  >  2:12815345/4‑96 (MQ=0)
            ttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATc                                                                                  >  1:13955104/4‑96 (MQ=0)
            ttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATc                                                                                  >  1:1567031/4‑96 (MQ=0)
            ttgaGGAAATCAATATGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATc                                                                                  >  2:10407757/4‑96 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGACCTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  2:9645304/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       >  2:11996787/4‑96 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       >  1:8073481/4‑96 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  1:7813398/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  2:7190038/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  2:8013257/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  1:349425/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       <  2:9086962/93‑1 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCAc                                                                       >  1:13641332/4‑96 (MQ=0)
                       aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTATTCGGCCATCCGCGACAGCAc                                                                       >  2:13529911/4‑96 (MQ=0)
                           tggcGGTTTACCAGTTCCGGAATACGCCGCTGCACCGACAGGGAGATCCCATCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGaaa                                                                   <  1:7426460/93‑1 (MQ=1)
                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCTTCCGCGACAGCACGAACGTGCAGAATGCg                                                       >  2:11875764/4‑96 (MQ=0)
                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCg                                                       <  1:350822/93‑1 (MQ=0)
                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCg                                                       <  2:14364516/93‑1 (MQ=0)
                                       agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCg                                                       >  2:13222784/4‑96 (MQ=0)
                                           cccgAAAACACCGCTGCACCGACAGGGCGATCCCGTCGAGAATTCTGGCAATTTCCCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGttt                                                   <  1:9227304/93‑1 (MQ=1)
                                               taagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               >  1:11698563/5‑96 (MQ=0)
                                               gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               >  2:1315531/5‑96 (MQ=0)
                                               gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               <  1:5646078/92‑1 (MQ=0)
                                               gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               >  2:7809199/5‑96 (MQ=0)
                                               gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               <  1:12998471/92‑1 (MQ=0)
                                               gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccac                                               <  2:1874918/92‑1 (MQ=0)
                                                aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccacc                                              >  1:5829660/4‑96 (MQ=0)
                                                aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccacc                                              >  1:13500242/4‑96 (MQ=0)
                                                aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccacc                                              >  1:12515180/4‑96 (MQ=0)
                                                aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccacc                                              <  2:13352121/93‑1 (MQ=0)
                                                aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTccacc                                              >  2:12237152/4‑96 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     >  2:6948116/4‑96 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     <  2:6150608/93‑1 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     >  1:5202562/4‑96 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     >  2:907486/4‑96 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     <  2:9834239/93‑1 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     <  1:13298159/93‑1 (MQ=0)
                                                         tgaaCCGACAGAGGAACTCCATCAAGAATGCTGACTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCg                                     >  2:14064122/4‑96 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          >  2:3089767/4‑95 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  2:11343814/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  2:8988568/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  2:98565/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  2:10263076/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  1:13780972/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          >  1:9958264/4‑95 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  1:6917716/93‑2 (MQ=0)
                                                                                    atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTCCtgg          <  1:6253554/93‑2 (MQ=0)
                                                                                       cTGGTAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       <  1:933021/96‑1 (MQ=14)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       <  2:6087022/96‑1 (MQ=31)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       <  1:9964010/96‑1 (MQ=31)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       <  2:2936890/96‑1 (MQ=31)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       >  2:1705572/1‑96 (MQ=31)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       <  1:4208379/96‑1 (MQ=31)
                                                                                       cTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg       >  1:316259/1‑96 (MQ=31)
                                                                                        tGGCAATTTCACCTGCGATCCGCGACAGCCCGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt      <  1:4670738/96‑1 (MQ=17)
                                                                                        tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt      <  1:9023633/96‑1 (MQ=32)
                                                                                        tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt      >  2:435074/1‑96 (MQ=32)
                                                                                        tGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt      <  2:13195430/96‑1 (MQ=32)
                                                                                         ggCAGTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCCGCTCCTGTGCGtt     <  1:13864976/96‑2 (MQ=12)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAGAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  1:12296508/96‑1 (MQ=18)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTTTGCGTc     <  2:13711337/96‑1 (MQ=18)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  1:2970495/96‑1 (MQ=33)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  2:12539864/96‑1 (MQ=33)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  2:944942/96‑1 (MQ=33)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  2:4685767/96‑1 (MQ=33)
                                                                                         ggCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc     <  2:1432304/96‑1 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    <  2:5857033/96‑1 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    <  2:3504097/96‑1 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    >  1:8307363/1‑96 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    <  2:2051770/96‑1 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    <  1:8396024/96‑1 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    >  1:9922125/1‑96 (MQ=33)
                                                                                          gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg    >  2:10553777/1‑96 (MQ=33)
                                                                                           cAATTTCCCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg   <  2:10240378/96‑1 (MQ=18)
                                                                                           cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg   >  1:5110131/1‑96 (MQ=33)
                                                                                           cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg   <  2:2644032/96‑1 (MQ=33)
                                                                                            aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc  <  2:3950470/96‑1 (MQ=33)
                                                                                            aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc  >  1:8125454/1‑96 (MQ=33)
                                                                                            aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc  >  1:4724844/1‑96 (MQ=33)
                                                                                            aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc  >  2:7819500/1‑96 (MQ=33)
                                                                                            |                                                                                               
ATGATATCCCGTTTCAGGAAATCAACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGC  >  NC_000913/1633925‑1634112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑