Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913730,7960GC56.1% ‑4.2 / 96.4 57G405A (GGC→GCC) rhsCRhs protein with putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (13/12);  new base C (19/13);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGG  >  NC_000913/730667‑730933
                                                                                                                                 |                                                                                                                                         
aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGaa                                                                                                                                  <  2:583805/139‑1 (MQ=255)
     gACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCggg                                                                                                                             <  2:464709/139‑1 (MQ=255)
      aCGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGggt                                                                                                                            <  2:265201/139‑1 (MQ=255)
      aCGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGggt                                                                                                                            <  2:1127628/139‑1 (MQ=255)
      aCGGGCGGGAGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGCAGCCGGGCTGAAGCgggg                                                                                                                            >  2:818503/1‑138 (MQ=255)
                      aCAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAacac                                                                                                            <  1:818503/139‑1 (MQ=255)
                                  ggCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGc                                                                                                >  2:840246/1‑139 (MQ=255)
                                      ggCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCa                                                                                            <  1:1095779/139‑1 (MQ=255)
                                                 cacGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCgg                                                                                                                              <  1:860932/94‑1 (MQ=255)
                                                 cacGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCgg                                                                                                                              >  2:860932/1‑94 (MQ=255)
                                                 cacGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAg                                                                                 >  1:480882/1‑139 (MQ=255)
                                                      aTCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGa                                                                            >  2:1071279/1‑139 (MQ=255)
                                                      aTCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGa                                                                            >  2:206771/1‑139 (MQ=255)
                                                          gTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcgcc                                                                        <  1:840246/139‑2 (MQ=255)
                                                                gaAGGACCGCATCACCGTCACCGACAGCGTGAACCGCCGTAAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGCAGCGGGTGGTGAAACAGGAACGCGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGc                                                                  <  1:766923/139‑1 (MQ=11)
                                                                gaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGc                                                                  >  2:105940/1‑139 (MQ=255)
                                                                  aGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAg                                                                >  1:918809/1‑139 (MQ=255)
                                                                  aGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAg                                                                >  1:312183/1‑139 (MQ=16)
                                                                   ggACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAgg                                                               >  2:685560/1‑139 (MQ=18)
                                                                     aCCGCGTCACCATCACCGACAGCCTGAACCGCCATGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCa                                                                                  <  1:145273/118‑1 (MQ=255)
                                                                     aCCGCGTCACCATCACCGACAGCCTGAACCGCCATGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCa                                                                                  >  2:145273/1‑118 (MQ=11)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGGACACGCAGGGTGAAGGCGGGCTGAAGc                                                                                                                                >  1:862175/1‑72 (MQ=255)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGc                                                                                                                                <  2:862175/72‑1 (MQ=255)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCt                                                             >  2:701350/1‑139 (MQ=12)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCt                                                             >  2:482619/1‑139 (MQ=17)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCt                                                             <  1:1071279/139‑1 (MQ=16)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCt                                                             >  2:1131345/1‑139 (MQ=16)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTG‑AAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCa                                                                                  <  1:1105190/117‑1 (MQ=255)
                                                                     aCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTG‑AAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCa                                                                                  >  2:1105190/1‑117 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACATGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAgg                                                         >  2:24783/1‑139 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAgg                                                         >  2:998533/1‑139 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAgg                                                         >  2:413404/1‑139 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCggg                                                                                                                             <  1:315177/71‑1 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCggg                                                                                                                             >  2:315177/1‑71 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGAGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTc                                                            >  2:475306/1‑136 (MQ=255)
                                                                         catcaccatcacCGACAGCCTGAACCGCCGAGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTc                                                            <  1:475306/136‑1 (MQ=255)
                                                                           tcaccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGc                                                       >  1:649387/1‑139 (MQ=255)
                                                                           tcaccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGc                                                       >  2:990166/1‑139 (MQ=12)
                                                                              ccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGcaca                                                    <  1:998533/139‑1 (MQ=255)
                                                                               catcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAg                                                   >  2:49118/1‑139 (MQ=12)
                                                                                       aCAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGc                                           >  2:416779/1‑139 (MQ=255)
                                                                                       aCAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGt                                                                                    >  2:521212/1‑98 (MQ=37)
                                                                                       aCAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGt                                                                                    <  1:521212/98‑1 (MQ=37)
                                                                                        cAGCCTGAACCGCAGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGca                                          <  1:338405/139‑2 (MQ=255)
                                                                                          gCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAGCACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcgccg                                                                       <  1:494332/108‑3 (MQ=255)
                                                                                          gCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcgccg                                                                       >  2:494332/1‑106 (MQ=255)
                                                                                           ccTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCgca                                       >  2:400910/1‑139 (MQ=16)
                                                                                                  cgccgTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGa                                <  2:767519/139‑1 (MQ=40)
                                                                                                     cgTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACgg                                               <  1:440704/121‑1 (MQ=255)
                                                                                                     cgTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACgg                                               >  2:440704/1‑121 (MQ=255)
                                                                                                     cgTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGa                                <  2:197320/136‑1 (MQ=40)
                                                                                                     cgTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGa                                >  1:197320/1‑136 (MQ=40)
                                                                                                                      gCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACg                                                                                          <  1:223222/61‑1 (MQ=255)
                                                                                                                      gCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACg                                                                                          >  2:223222/1‑61 (MQ=255)
                                                                                                                             gAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGa     <  1:413404/139‑1 (MQ=255)
                                                                                                                              aaGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGAt                                             >  2:700182/1‑98 (MQ=39)
                                                                                                                              aaGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGAt                                             <  1:700182/98‑1 (MQ=39)
                                                                                                                                ggCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACgg  >  2:1068347/1‑139 (MQ=39)
                                                                                                                                 |                                                                                                                                         
ACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAATACAGTCCGGATGTGGTGACGG  >  NC_000913/730667‑730933

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