Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_pgi526,7920AC44.7% 17.0 / 28.1 40S1178R (AGC→CGC) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (5/16);  new base C (11/6);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 2.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAC  >  NC_000913_3_vch_pgi/526646‑526908
                                                                                                                                                  |                                                                                                                    
ggCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAgc                                                                                                                    <  1:3134/149‑1 (MQ=255)
 gCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAATGCGGGGTGgcgcc                                                                                                                   >  2:145015/1‑148 (MQ=255)
   ttGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                 <  2:125977/149‑1 (MQ=255)
   ttGCGCAGTGCGGGCGGAGGGTGGAGACACGGGCCGCACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTTCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGcgcg                                                                                                                 <  2:341817/149‑1 (MQ=17)
   tggcgcAGTGCGGGCGGACGGTTGAGCAACTGGTCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATTTGTTTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                 <  2:55943/147‑3 (MQ=11)
     gcgcAGTGCGGGCTGACGGTGGAGCAGCTGCGCACACAGGTCGCGCCGGAATACACACCGGCGCGAAAAGCTCATCTTCATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGa                                                                                                               <  2:180829/149‑1 (MQ=11)
              gggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgataa                                                                                                      >  2:990/1‑147 (MQ=255)
                  tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                  <  2:137982/149‑1 (MQ=255)
                  tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCCAAGCCGGCAATACGACGTGGGGCGGGGAGTCCGATAAATgg                                                                                                  >  2:150493/1‑149 (MQ=11)
                  tGACGGTGGACCAACGGGCCAGACAGGTGGAGCCGGAATGCACACCGGCGCGAACAGCTCGTCTTTATCGCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAACACGGCAATACGGCGTGGAGCGCGGAATATGATGAATgg                                                                                                  <  1:338426/149‑1 (MQ=17)
                   gACGGTGGAGCAACTGGCCGGACAGGTGGGGCCGGAATACACATCGGCGCGAAAAGCTCATATTTATCACTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGCGGAATATGATGAATggg                                                                                                 <  2:345661/149‑1 (MQ=14)
                       gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCaa                                                                                             >  1:343195/1‑149 (MQ=255)
                           aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAg                                                                                         <  2:347984/149‑1 (MQ=255)
                           aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAg                                                                                         <  1:208871/149‑1 (MQ=255)
                            gCAACGGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGCTCCTCTATATCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc                                                                                        <  2:389337/149‑1 (MQ=17)
                               aCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATACGCCAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCtga                                                                                     >  1:90691/1‑147 (MQ=11)
                                 tGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAAt                                                                                   >  1:87150/1‑149 (MQ=255)
                                   gcacaaCAGGTGGAGCCGGAATCCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGCGGGGAGAGAGTCATGAAGGGGGCAACCAGCTTAATga                                                                                 >  2:374887/6‑149 (MQ=11)
                                   gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                 <  1:990/149‑1 (MQ=255)
                                      aggcAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATTTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagga                                                                              <  2:381914/146‑1 (MQ=255)
                                          cggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCGATACGGCTTGGCGCCGGTATTTTTATGAATGGGGCAACCAGCTTAATGAGGAGAAc                                                                          <  1:290850/148‑1 (MQ=11)
                                           ggTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATAGGGCAACAAGCTTAATGAGGAGAAcc                                                                         >  2:54786/1‑149 (MQ=11)
                                            gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCACCGAAGACGGCAATACGGCGTGGCGGGGGGAGTCGGAAGAATGGGGCAACCAGCCTAATGTGCAGAAccc                                                                        >  2:156022/1‑149 (MQ=11)
                                             tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                       >  1:362294/1‑149 (MQ=255)
                                              ggAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGTGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACAGCCTCGCGCGTGGTAAATTCTGAATAGGACAACCAGCTTCGTTAGGAGAAccctc                                                                      >  2:45656/1‑147 (MQ=14)
                                                 gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGGACCCGcatc                                                                   >  1:312935/1‑149 (MQ=11)
                                                 gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGAGGACTGCCGCTGGCGCTTATCAGCGAAGACGCCAATACGGCGTGGCGCGGGGAGAATGATGAATGGCGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                   >  1:159601/1‑149 (MQ=11)
                                                     gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                              >  1:330173/1‑134 (MQ=255)
                                                     gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga                                                                              <  2:330173/134‑1 (MQ=255)
                                                     gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCGTGTCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacct                                                               >  1:308247/1‑146 (MQ=11)
                                                           acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa                                                         >  1:140545/1‑149 (MQ=255)
                                                                 gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCTATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGGGGACCCGCATCATGTGTATCAGCCGTa                                                   >  1:74763/1‑149 (MQ=255)
                                                                  cgcgAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTAc                                                  <  1:375834/149‑1 (MQ=255)
                                                                            tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg                                        <  1:320609/149‑1 (MQ=255)
                                                                             cATCTTTGTCACTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCTGCATCATGTGTATCAGCCGTACCGTCTGCCAgg                                       <  1:36413/149‑1 (MQ=255)
                                                                                       aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                 <  2:89557/65‑3 (MQ=255)
                                                                                       aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGcgcggg                                                                                                                 >  1:89557/1‑63 (MQ=9)
                                                                                          gcgcCCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGACAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                          <  1:358659/149‑1 (MQ=255)
                                                                                                    gggACTGCCGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAg                <  1:305948/149‑1 (MQ=255)
                                                                                                          gccgccCGCCCTAATCCGCGAAGACGGCCAGACGGCCTGGGGCCCGTAATAAAATGAATTGCGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGt          <  2:336621/149‑1 (MQ=14)
                                                                                                           ccGCTGGCGCTTATCAGCGAAGACGACAATACGGCGTGGAGCGCGGGATATGATGAATGGGGCAACCGGCTTAATGAGGAGAACACGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGGAGCATGATGAGGAATAAGCGCTGTa         >  1:354111/1‑149 (MQ=17)
                                                                                                            cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTAc        <  1:63822/149‑1 (MQ=255)
                                                                                                                ggCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACtata    <  1:198504/149‑1 (MQ=255)
                                                                                                                  cgcTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAc  <  1:213663/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                    
GGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAC  >  NC_000913_3_vch_pgi/526646‑526908

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: ‑