Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,8180GC55.2% ‑4.1 / 36.9 29intergenic (‑167/+62)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base C (11/5);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAG  >  NC_002947/2282673‑2282966
                                                                                                                                                  |                                                                                                                                                    
cTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGc                                                                                                                                                  <  1:509518/150‑1 (MQ=255)
  gAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt                                                                                                                                                <  1:207574/150‑1 (MQ=255)
      ggCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                                                                                            >  1:942737/1‑150 (MQ=255)
       gcgcTCTTGCTGGCACTTAGAAGGCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                           <  2:942737/150‑1 (MQ=255)
          ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTG‑CAAACACGTGGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                                        <  1:796134/150‑1 (MQ=38)
                  ggCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACtgtg                                                                                                                                <  1:963252/150‑1 (MQ=255)
                                           tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgtg                                                                                                                                             <  1:518636/112‑1 (MQ=25)
                                                                                              tccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccac                                                     >  1:77057/2‑150 (MQ=14)
                                                                                                         cGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGTCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAccc                                          <  1:449055/150‑1 (MQ=255)
                                                                                                                           cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGCACAAGATATATTTGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg                        >  2:624864/1‑150 (MQ=37)
                                                                                                                                  cAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGa                                                                         >  2:575685/1‑94 (MQ=255)
                                                                                                                                  cAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGa                                                                         <  1:575685/94‑1 (MQ=255)
                                                                                                                                  cAACTGCACAGTCTTGCAAGCCTGTGCAGTTGGACTGTGTGAGCGCGCGTGCCCGCGAAGGGGCCCGAGCAGGCGACACCAGTTACCGGGCAGATAAGGCCAGTCATCACCCACCAGAAACACCCGTTCTCCCTCCTTCCCGGCACCATc                 >  2:933747/1‑150 (MQ=255)
                                                                                                                                   aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                >  2:73336/1‑150 (MQ=255)
                                                                                                                                   aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                >  2:575567/1‑150 (MQ=37)
                                                                                                                                   aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                >  1:648814/1‑150 (MQ=37)
                                                                                                                                   aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCg                >  1:44354/1‑150 (MQ=37)
                                                                                                                                   aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTACCAGGCACCATCg                >  1:408583/1‑150 (MQ=25)
                                                                                                                                    aCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccacc                                                    <  1:166003/113‑1 (MQ=21)
                                                                                                                                    aCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATcacccacc                                                    >  2:166003/1‑113 (MQ=21)
                                                                                                                                    aCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCgg               >  2:335819/1‑150 (MQ=37)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCa                                                                 >  1:292193/1‑98 (MQ=18)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCa                                                                 <  1:325692/98‑1 (MQ=18)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCa                                                                 <  2:292193/98‑1 (MQ=18)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCa                                                                 >  2:325692/1‑98 (MQ=18)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcc                                  <  1:135072/129‑1 (MQ=25)
                                                                                                                                      tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcc                                  >  2:135072/1‑129 (MQ=25)
                                                                                                                                       gctcaGTCTTGGAAGCCTGTGCAGTTGGAATGTGGGACCGGGCGTGCCCGCGAAGGGGCGTTATCAGGCGACATCAGTTACCTGGCATATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCAACATTGGCAt            <  2:127750/147‑1 (MQ=255)
                                                                                                                                               ttGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCAcc    >  2:651603/1‑150 (MQ=255)
                                                                                                                                                 ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                           <  2:167060/45‑1 (MQ=255)
                                                                                                                                                 ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                           >  1:167060/1‑45 (MQ=255)
                                                                                                                                                 ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCTGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg  >  1:813770/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGT‑GTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAG  >  NC_002947/2282673‑2282966

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