Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,8640TA20.0% 84.5 / 17.2 40L541F (TTA→TTTydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (17/15);  new base A (4/4);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGG  >  NC_000913/1584726‑1585005
                                                                                                                                          |                                                                                                                                             
cacaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccgg                                                                                                                                     <  2:179593/149‑3 (MQ=255)
               aTTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAatca                                                                                                                      >  8:100168/1‑149 (MQ=255)
                       tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAtt                                                                                                              <  4:1822/149‑1 (MQ=255)
                          gATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccgg                                                                                                                                     >  8:23892/1‑121 (MQ=255)
                          gATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccgg                                                                                                                                     <  7:23892/123‑3 (MQ=255)
                                aaTCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTc                                                                                                      >  1:89013/1‑148 (MQ=255)
                                  tCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGaa                                                                                                                         >  3:69242/1‑127 (MQ=255)
                                  tCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGaa                                                                                                                         <  4:69242/127‑1 (MQ=255)
                                  tCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                                                                                                   >  3:48329/1‑149 (MQ=255)
                                      ccaccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAAc                                                                                                   <  2:89013/145‑1 (MQ=255)
                                         ccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCCGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                    <  1:89405/109‑4 (MQ=255)
                                         ccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                    >  2:89405/1‑106 (MQ=255)
                                                             aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                                                                                                    >  5:22086/1‑105 (MQ=255)
                                                             aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcatc                                                                                                                    <  6:22086/105‑1 (MQ=255)
                                                                ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgc                                                                                      >  4:73442/1‑132 (MQ=255)
                                                                ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgc                                                                                      <  3:73442/132‑1 (MQ=255)
                                                                 tctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                    >  8:72471/1‑82 (MQ=255)
                                                                 tctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                    <  7:72471/85‑4 (MQ=255)
                                                                     gggTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGt                                                                   <  7:100168/146‑1 (MQ=255)
                                                                               gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGt                                                                   <  5:139017/136‑1 (MQ=255)
                                                                               gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGt                                                                   >  6:139017/1‑136 (MQ=255)
                                                                                  tgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTcc                                                   >  6:104761/1‑149 (MQ=255)
                                                                                    cgcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCa                                                 <  4:48329/149‑1 (MQ=255)
                                                                                             cgcgACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCt                                                         <  2:68103/132‑1 (MQ=255)
                                                                                             cgcgACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCt                                                         >  1:68103/1‑132 (MQ=255)
                                                                                                     cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCa                                                                      >  4:155065/1‑111 (MQ=255)
                                                                                                     cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCa                                                                      <  3:155065/111‑1 (MQ=255)
                                                                                                      tGCACACTCAGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                               >  8:53192/1‑149 (MQ=255)
                                                                                                             tctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGt                                                                                                >  6:60964/1‑77 (MQ=255)
                                                                                                             tctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGt                                                                                                <  5:60964/77‑1 (MQ=255)
                                                                                                             tctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCt                              <  3:87138/143‑1 (MQ=255)
                                                                                                             tctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCt                              >  4:87138/1‑143 (MQ=255)
                                                                                                              ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCa                                                 <  6:113410/123‑1 (MQ=255)
                                                                                                              ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCa                                                 >  5:113410/1‑123 (MQ=255)
                                                                                                               tGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTAc                                                                                             <  6:77480/78‑1 (MQ=255)
                                                                                                               tGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTAc                                                                                             >  5:77480/1‑78 (MQ=255)
                                                                                                                          ttACACCGGCGGGTTTTAACACGCCACGAGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCaa               >  8:59169/1‑145 (MQ=255)
                                                                                                                              acCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGAGGTCAACAGAt          >  8:143784/1‑146 (MQ=255)
                                                                                                                               gtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTCATGTCCATCTCGCTACGTCCCAGGACCGGCAGAATATAg                                >  8:128543/3‑123 (MQ=255)
                                                                                                                                   cggGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCgg  <  7:53192/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
CACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGG  >  NC_000913/1584726‑1585005

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