Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,7960CG44.7% 29.2 / 40.8 39A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base G (10/7);  total (25/14)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGG  >  NC_000913_3_bme_pgi/526648‑526917
                                                                                                                                                    |                                                                                                                         
cTTGCGCAGTGCGGGCCGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgcgc                                                                                                                           <  2:571446/149‑1 (MQ=255)
          gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                 <  1:93224/149‑1 (MQ=255)
              gCTGACGTTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                             <  2:422097/149‑1 (MQ=11)
              gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                             <  1:250421/149‑1 (MQ=255)
                tGACGGTGGAGCAACTGGCCAGAAAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATTAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                           <  1:569677/149‑1 (MQ=11)
                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACTGGCGTGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg                                                                                                          <  1:365417/149‑1 (MQ=255)
                 gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                          <  1:323684/149‑1 (MQ=255)
                    ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTAACAGCGCAGACAGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCa                                                                                                       >  1:262723/1‑149 (MQ=255)
                      tGGAGCAACTTGCCAGAGAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAc                                                                                                     <  2:477779/149‑1 (MQ=255)
                         aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAg                                                                                                  <  2:567251/149‑1 (MQ=255)
                          gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCAGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAAGATGATGAATGGGGCGACCAGc                                                                                                 >  1:555215/1‑149 (MQ=255)
                                  acAGAAAGGTGGAGCCGTAATACACACCGGCTCTAAAAGCTTATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                         <  1:3220/148‑1 (MQ=11)
                                   cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACTGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGAGAGCATGATGAAAGGGGCAACCAGCTTAATgggg                                                                                        >  2:38096/1‑149 (MQ=11)
                                   cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAGGATGAATGGGGCAACCAGCTTCATgagg                                                                                        >  2:578403/1‑149 (MQ=11)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGTGAAACAAGCTTAAggagga                                                                                       >  1:74037/1‑149 (MQ=255)
                                    aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                                       >  1:420564/1‑149 (MQ=255)
                                      aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATCAATGGGGAAACCAACGTAATGAGgaga                                                                                     >  1:190257/1‑149 (MQ=11)
                                      aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCTCGGAATATAAAGAATGGAGCAACCAGCTTTAtaataaaa                                                                                     >  1:259419/1‑142 (MQ=17)
                                      aCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAAACGGCAATACCGCGTGGCGCGGGGAGTAAGATGAATGGGGCAACCAGCTAAAAgaggaaa                                                                                     >  1:190648/1‑147 (MQ=11)
                                          gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAAccc                                                                                 >  1:423245/1‑149 (MQ=255)
                                           tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGCGTATGATGAATGGGGCAACCAGCGTGAGGAGGAGAACCCg                                                                                >  2:555787/1‑149 (MQ=11)
                                             gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGTTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATAGGGAGAACCAGca                                                                              >  1:107352/1‑149 (MQ=11)
                                                  ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGACTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                                                         >  2:463110/1‑149 (MQ=255)
                                                         acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGAAGCGCGGAATATGAGGAGTCGGGCAACCAGCTCAATGAGGAGGACACGCGTCAACTGTAtac                                                                  >  2:358508/1‑147 (MQ=17)
                                                           acCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCAGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAACGGGGCAACCGCCTCAATGAGGAGAACCCGCATCACGGGTATAAga                                                                >  2:339413/1‑148 (MQ=17)
                                                               gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTa                                                            >  2:114904/1‑149 (MQ=255)
                                                               gcgcGAAAAGATCATCTATATAAATGAGAAAACAGGGGACTGCAGGTGGAGATTATCAGCGAAGACGAAAATACGGAGTGGAGAGAGGAATATGATGAATGGGGCAAACAGCTTAATGAGGAGAACACGAATCATGTGTATCAGGCGTa                                                            >  2:371580/1‑149 (MQ=255)
                                                                  cgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCg                                                         >  1:467137/1‑149 (MQ=255)
                                                                        gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                           <  2:369556/93‑1 (MQ=11)
                                                                        gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                           >  1:369556/1‑93 (MQ=14)
                                                                          tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                       >  2:13251/1‑111 (MQ=255)
                                                                          tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                                       <  1:13251/111‑1 (MQ=11)
                                                                                  aTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGCGGAGCGCGGGATATGATGAATGGAGCAACCAGCTTAATGCGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcaa                                         >  2:555434/1‑148 (MQ=17)
                                                                                     aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCat                                      <  1:114904/149‑1 (MQ=255)
                                                                                                cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                                            >  2:337535/1‑132 (MQ=255)
                                                                                                cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcag                                            <  1:337535/132‑1 (MQ=255)
                                                                                                       tGCCGCTGGCGCTTATCAGCGAAGACGGCACTCCGAAATGTCGCGGGGAGAATGATGAATGGGGCAAACAGCTTAATGAGGAGAAACCGCATCACCTGAACCAGACGTAACGTATGCAAGGGCAGCAGCAAGATGAGGAGTCGGGGCTg                    >  1:591544/1‑149 (MQ=255)
                                                                                                               gcgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCTGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTCCTATaa            >  1:71290/1‑149 (MQ=255)
                                                                                                                               aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                   >  2:73563/1‑110 (MQ=255)
                                                                                                                               aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                   <  1:73563/110‑1 (MQ=255)
                                                                                                                                                  cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg  <  2:500418/124‑1 (MQ=255)
                                                                                                                                                  cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg  >  1:500418/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                         
CTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGG  >  NC_000913_3_bme_pgi/526648‑526917

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