Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi526,8000AG42.5% 29.1 / 33.5 41E1180E (GAA→GAG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (15/8);  new base G (10/7);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGG  >  NC_000913_3_bme_pgi/526658‑526917
                                                                                                                                              |                                                                                                                     
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                                                                                                           >  1:369556/1‑93 (MQ=14)
                                                              gCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                           <  2:369556/93‑1 (MQ=11)
                                                                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                                   <  1:73563/110‑1 (MQ=255)
                                                                                                                     aCGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                   >  2:73563/1‑110 (MQ=255)
                                                                                                                                        cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg  <  2:500418/124‑1 (MQ=255)
                                                                                                                                        cgcgGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg  >  1:500418/1‑124 (MQ=255)
                                                                                                                                              |                                                                                                                     
GCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGG  >  NC_000913_3_bme_pgi/526658‑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: ‑