Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,595,3160CA40.5% 56.1 / 10.3 38F129L (TTC→TTA) PP_3171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/20);  new base A (15/0);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGA  >  NC_002947/3595181‑3595449
                                                                                                                                          |                                                                                                                                     
gTTGGTGGGGCCCGAACAGTTTGCCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGGG‑CGCCAGGCACGCCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTc                                                                                                                                       <  2:274723/136‑1 (MQ=255)
    gtggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TTGGCGGGGCGCGCCCCCCGGGGGGGATTTGGGc                                                                                                                                    >  1:436088/1‑136 (MQ=255)
     tggtggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCa                                                                                                                                  <  2:823446/136‑1 (MQ=255)
     gggtggCCGTACAGTTTGGCACCCTCAATACTTTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTT‑CGCCAGGCACGTCCGGGCCGACCGTTCGGTCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCa                                                                                                                                  <  1:230080/135‑1 (MQ=255)
          gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCCCGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTACGCCAggtcc                                                                                                                             >  2:141361/1‑133 (MQ=255)
          gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GGGCCGCGCCGCGCCCGCGCGGGCTTTAGGCCAGGACa                                                                                                                             >  2:288108/1‑136 (MQ=255)
          gCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCCCGGCCCGGCGGGGTTTAGGCCCgggcc                                                                                                                             >  2:693635/1‑133 (MQ=255)
               aCAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGcc                                                                                                                        <  2:400754/136‑1 (MQ=255)
                ctgTTTGGCACCCTCAATACATTGCACCCAGGGCGTCTCGACCTGGGCCTGGGCCGGG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCg                                                                                                                       <  1:2795/134‑1 (MQ=255)
                cattttGGCACCCTAATTACATTGCACCCAGGGCGTATCGACGGGGGCCTGGGCGGTG‑CGCCTGGAACGTCCGGGCCGACGGTTCGGGCT‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCg                                                                                                                       <  1:46012/133‑1 (MQ=255)
                cAGTTTGCCACCCTCAATCCATTGCACCCAGGTCGTAGCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCg                                                                                                                       <  1:464010/136‑1 (MQ=255)
                 cgtctGGCCCCCTCAAACCATTGCACCCAGGTCGTAGCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGa                                                                                                                      <  1:375770/132‑1 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCCGGCCGGGCTTTAGGCCCCGGCAAAGCCGCGCTGCGGGCTTAACAGGCCGACaa                                                                                              >  1:821011/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCCGCGCGGGCTTTAGGCCAGGCCCTCGCCGAGCTGCGTGATTACCTGGCCGACaa                                                                                              >  2:793614/1‑136 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGGGGGGGATTTAGGCCCGGGCCTCGCCGCGCGGCGGGGTTTCCTGGCCgcccc                                                                                              >  2:821419/1‑132 (MQ=255)
                                         aCCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCCCCAAAGCGGGATTTAGGCCAGGACCTCGCCGCGCTGGGGGATTAACTGGCCGCCaa                                                                                              >  1:350785/1‑136 (MQ=255)
                                          cccAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCCCGGCGCCGCGCCAGAGCGGGATTTAGGCCAGGAAAAAGCCGAGCAGCGGGTTTACCTGGGCGACAAc                                                                                             >  2:564195/1‑136 (MQ=255)
                                              gggcGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGCCCGTTCGGGCA‑TT‑GGGCGGCGCCGCCGCCGTGCGGGCTTTAGGGCAGGACAACGGCGCGCTGGGGGATTTACTGGGCGCCCACGGcc                                                                                         >  2:564428/4‑135 (MQ=255)
                                                      cccccggggcccggggcGGGCCGCCCGGCCCGCCCGGGCCTACCGTGCGGGCT‑TG‑GCGGGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCg                                                                                  <  2:295946/130‑1 (MQ=255)
                                                           ggggcctgggccGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTgcg                                                                            <  1:56345/135‑1 (MQ=255)
                                                            gggcctgggccGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGGTTTAGGCCCGGGCATCGCCGCGCTGGGGGATTACCTGGGCGCCCACGGCAAACGCCgggggggc                                                                           >  2:302623/1‑131 (MQ=255)
                                                            gggcctgggccGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCACAGCGGGAATTAGGCCAGGACAACGCCCACCTGCGTGATTAACTTGCCGACAACCGCAAACCCCGGGTgccc                                                                           >  2:33153/1‑134 (MQ=255)
                                                                ctgggccGTG‑CGCCTGGCACGTCCGGGCCGACCGTTCGGGCA‑TG‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCgcggcg                                                                       <  1:190384/136‑1 (MQ=255)
                                                                  gggccGGC‑CCCCAGGCCCGCCCGGGCCTCCCGTGCGGCCA‑TG‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTa                                                                     <  2:720953/136‑1 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCCGAGCGGGTTTTCGGCCAGGACCTCGCCGAGCGGCGTGATTACCTGGCCGCCAACGGCCAACGCCCGGGGCGCGGCGTACCGGGTgcg                                                            >  1:763188/1‑136 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGCGCGGGATTTCGGCCAGGCCATCGCCGAGCTGCGTGATTACCTGGCCGCCACCGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:168932/1‑136 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGGGGGGTTTAGGGCAGGGCATCGCCGGGCTGCGTGATTTCCTGGCCGACAACGGCAAACCCCGGGTGCGCGGCGTACCGGGGgcg                                                            >  1:802931/1‑136 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGTTTTAGGCCCGGGCATAGCCGAGCGGCGTGATTTCCAGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTgcg                                                            >  1:724647/1‑136 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTAGGGCAGGAACTCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCCGGTgcg                                                            >  1:187552/1‑136 (MQ=255)
                                                                            gcCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCCGCGCCCCCCCAGAGCGGGATTTAGGCCAGGACAACCCCGAGCTGCCAGATTAAAAGGCCGACCACGCCAAAAGCCCGGTGACCCGCGTACCCGGTgcc                                                            >  2:731332/1‑135 (MQ=255)
                                                                                   aCGTCCGGGCCGACCGTGCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATATCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGAtg                                                     <  2:204501/136‑1 (MQ=255)
                                                                                     ggcggggCCGA‑CGTTGGGGGATTT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGAtgtg                                                   <  2:821011/132‑1 (MQ=255)
                                                                                           cccGACCGTTCGGGAT‑TG‑GCGGGGCGCGGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGgtg                                             <  1:397083/135‑1 (MQ=255)
                                                                                                ccGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGAt                                        <  2:642178/136‑1 (MQ=255)
                                                                                                ccGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGAt                                        <  2:436088/136‑1 (MQ=255)
                                                                                                     ggggCA‑TG‑GCGGGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCg                                   <  2:763188/135‑1 (MQ=255)
                                                                                                                                   gggATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAACCTGTTCGGTGCTGACCTggc       <  1:821419/136‑1 (MQ=255)
                                                                                                                                       tttCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGgcggcg   <  1:222894/136‑1 (MQ=255)
                                                                                                                                        ttCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGa  <  2:726650/136‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
GTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTG‑CGCCAGGCACGTCCGGGCCGACCGTTCGGGCA‑TT‑GCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACGCCGGGTGCGCGGCGTACCGGGTGCGCACGATGTGCCGGTGTGGATCCTCGGTTCCAGCCTGTTCGGTGCTGACCTGGCGGCGA  >  NC_002947/3595181‑3595449

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