Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,8170GC34.4% 38.5 / 27.6 32intergenic (‑166/+63)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/10);  new base C (6/5);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAG  >  NC_002947/2282675‑2282966
                                                                                                                                              |                                                                                                                                                      
gAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt                                                                                                                                                 <  2:556214/150‑1 (MQ=255)
gAATCGCGCTTTTGCTGGCGCGGATGAGTTATTTTAATGGCTACCGGAAACAGATCTTTGCGACAAGTCTTGGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCt                                                                                                                                                 <  1:255367/150‑1 (MQ=38)
  aTGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgt                                                                                                                                               <  2:405935/150‑1 (MQ=255)
      cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                                                           <  2:779044/150‑1 (MQ=37)
        ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                                         <  1:596739/150‑1 (MQ=255)
                      cAGGAGTCATTTTAATGGCCACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCg                                                                                                                           <  2:569091/150‑1 (MQ=25)
                              aTTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGc                                                                                                                            <  1:764633/141‑1 (MQ=25)
                              aTTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGc                                                                                                                            >  2:764633/1‑141 (MQ=25)
                                     tGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                            <  1:308057/150‑1 (MQ=255)
                                         tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                            >  1:439829/1‑114 (MQ=255)
                                         tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCa                                                                                                                                            <  2:439829/114‑1 (MQ=255)
                                            cGGAAACAGATTTTTTCGCTAGGTCGGTTAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCgg                                                                                                     >  2:496617/1‑150 (MQ=25)
                                                                      gtgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACAt                                                                                        >  1:47670/1‑136 (MQ=255)
                                                                      gtgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACAt                                                                                        <  2:47670/136‑1 (MQ=255)
                                                                                               gCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGa                                                 >  1:327947/1‑150 (MQ=255)
                                                                                                cccctTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTgg                                                                                                                                >  1:279640/4‑71 (MQ=255)
                                                                                                cccctTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTgg                                                                                                                                <  2:279640/68‑1 (MQ=255)
                                                                                                    ttCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCCGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAATTACCGGGCGGCTGGCGCCCGGCAGCAGCCACCGCAAacac                                            >  2:770735/1‑150 (MQ=21)
                                                                                                                     tCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTAcc                                                                                >  1:31722/1‑97 (MQ=255)
                                                                                                                     tCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTAcc                                                                                <  2:31722/97‑1 (MQ=255)
                                                                                                                        cacaGTCCAACTGCACAGTCTTGGAAGGCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGaa                                                                                                             <  2:463552/66‑1 (MQ=14)
                                                                                                                        cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGaa                                                                                                             >  1:463552/1‑66 (MQ=255)
                                                                                                                         acaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAggg                                                                    >  2:277889/1‑105 (MQ=18)
                                                                                                                         acaGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAggg                                                                    <  1:277889/105‑1 (MQ=18)
                                                                                                                               cAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                 >  1:514944/1‑150 (MQ=37)
                                                                                                                                aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCGGTGGCTGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc                 >  1:389771/1‑150 (MQ=18)
                                                                                                                                   tGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa             >  1:701461/1‑150 (MQ=255)
                                                                                                                                   tGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACCTCAGTTACCGAGCAGAGAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa             >  1:696526/1‑150 (MQ=37)
                                                                                                                                      acaGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCa          >  1:504846/1‑150 (MQ=255)
                                                                                                                                        aGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAgg        >  1:788361/1‑150 (MQ=255)
                                                                                                                                          tCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGcg      >  1:289754/1‑149 (MQ=255)
                                                                                                                                             tGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCa   <  2:701461/150‑1 (MQ=255)
                                                                                                                                              ggAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAg  >  2:279224/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
GAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGA‑GCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATCAGGCACCAG  >  NC_002947/2282675‑2282966

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