Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_pgi526,7960CG30.9% 58.4 / 35.4 42A1179G (GCG→GGG) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (18/11);  new base G (5/8);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGG  >  NC_000913_3_vch_pgi/526655‑526937
                                                                                                                                             |                                                                                                                                             
aGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGGGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAttat                                                                                                                                        <  2:37426/149‑1 (MQ=255)
   gCGGGCTGACGGTGGAGCGACTGGACAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTAGCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat                                                                                                                                     <  2:376487/149‑1 (MQ=255)
         tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATgg                                                                                                                               <  1:414646/149‑1 (MQ=255)
          gACGGTGGAGTACGTGGCAATACAGGTGGAGCCGGAAGACACACCGGCGGGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg                                                                                                                              <  1:42052/149‑1 (MQ=11)
             ggTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCAATGCGAACACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCa                                                                                                                           >  1:77003/1‑149 (MQ=11)
                 gAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCa                                                                                                                       >  1:360761/1‑149 (MQ=255)
                  aGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCACATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCCATACGGCGTGGAGCGCGGAATATGATGAGTGGGGccaccat                                                                                                                      >  2:320951/1‑148 (MQ=17)
                   gCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGc                                                                                                                     <  1:369467/149‑1 (MQ=255)
                          ggcgcgcAGGTGGAGTCGGGATAAACACCGGCGCGAAAATCTCATCTTTATCCCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATga                                                                                                              <  2:75302/143‑1 (MQ=17)
                          gcaggACAGGTTGGGCCGGAAGACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                              <  2:91058/145‑1 (MQ=11)
                           ccAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgag                                                                                                             <  1:27287/149‑1 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                            >  1:172807/1‑149 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg                                                                                                            >  1:243463/1‑149 (MQ=255)
                            cAGACAGGTGGAGCCGGAATACACACCGGCGCCAAAACCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGGAGACGGCAATACGGCGTGGCGCGCGAACTAGCATGAAGAGGACAACCAGCTCAATgaag                                                                                                            >  2:398217/1‑147 (MQ=14)
                              gACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagggg                                                                                                          >  1:149254/1‑147 (MQ=255)
                                cAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGaa                                                                                                        <  1:240942/149‑1 (MQ=255)
                                    tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCg                                                                                                    >  1:361344/1‑149 (MQ=255)
                                       aGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCCGGGCGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcaa                                                                                                 >  2:135905/1‑148 (MQ=11)
                                        gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                              <  1:101148/135‑1 (MQ=255)
                                        gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga                                                                                                              >  2:101148/1‑135 (MQ=255)
                                        gCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatc                                                                                                <  2:215774/149‑1 (MQ=255)
                                          cGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCACTACGGCGTGGCGCGGGGAGTATGATGAATGGGGCGACCAGCTTAATGAAGGGACCCCGcctcca                                                                                              >  2:157480/1‑147 (MQ=11)
                                           ggAATACACACCGGCGCGAAAAGCTGATCTTTAGCACTGCGACCAGCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAtg                                                                                             <  2:26936/149‑1 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAata                                                                                                                                         >  1:192764/1‑104 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAata                                                                                                                                         <  2:192764/104‑1 (MQ=255)
                                            gAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGCGGAGAACCCGCATCAtgt                                                                                            >  1:41006/1‑149 (MQ=255)
                                                  acacCGGCGCGAAATGCTCAGCGTGATGCTCGGTACCACCGGGGACTGTCGGTGGCGCTTATCAGCGAAGACGGCAATACGTCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTGAATGAGGAGAACTCGCATCATGTGTATCa                                                                                      <  2:100709/149‑1 (MQ=17)
                                                            gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                            >  1:519038/1‑149 (MQ=255)
                                                            gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGt                                                                            >  1:511923/1‑149 (MQ=255)
                                                                         ttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGACAACCAGCTGAATGAGGAGGACCCGCATCGTGTGTATCAGCCGGACCCTCGACCAGGgccaa                                                               >  2:329460/1‑146 (MQ=17)
                                                                          tATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGTAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca                                                              <  1:78861/149‑1 (MQ=255)
                                                                            tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagc                                                            >  1:45833/1‑149 (MQ=255)
                                                                            tCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGACGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagcagc                                                            <  1:159274/149‑1 (MQ=255)
                                                                               cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                                                         >  1:203469/1‑149 (MQ=255)
                                                                                        ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAAt                                                >  1:78495/1‑149 (MQ=255)
                                                                                                  ccGCGGGCGGTTATCAGCGAAGACGGCAACACCGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGGGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTa                                      <  1:430804/149‑1 (MQ=17)
                                                                                                         cgcTTATCAGGGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAAc                               <  2:174065/149‑1 (MQ=255)
                                                                                                                 aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                             >  2:83446/1‑63 (MQ=255)
                                                                                                                 aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgag                                                                                                             <  1:83446/63‑1 (MQ=255)
                                                                                                                  gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg                      >  2:406638/1‑149 (MQ=255)
                                                                                                                          ggCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGa              >  1:159974/1‑149 (MQ=255)
                                                                                                                                      tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAggg  >  2:318100/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                             
AGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGG  >  NC_000913_3_vch_pgi/526655‑526937

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