Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913115,7970CA41.9% 33.8 / 10.2 31E435* (GAA→TAA) hofBT2SE secretion family protein; P‑loop ATPase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base A (0/13);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 1.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACAC  >  NC_000913/115649‑115946
                                                                                                                                                    |                                                                                                                                                     
ctGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTc                                                                                                                                                       <  2:578543/149‑1 (MQ=255)
 tGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaa                                                                                                                                                     >  2:44112/1‑150 (MQ=255)
    ccATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCaaaa                                                                                                                                                   >  1:866660/1‑149 (MQ=255)
        aGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTa                                                                                                                                             <  1:2995011/151‑1 (MQ=255)
                     gcgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                                                                                                 <  2:153325/150‑1 (MQ=255)
                      cgcATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGcc                                                                                                                                 <  1:2606192/149‑1 (MQ=255)
                            gCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTcgtgc                                                                                                                         >  2:3640807/1‑151 (MQ=255)
                                   ccGGTTATGCCATGCCAGCGCCAGGGTTGTTTACTCGCCAGGCGGCACCCCCAGTACGCGGATTAACTCTTAAAAGGTGGTAAAGCTTTGCTCTCCGGCCGGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGtt                                                                                                                  <  2:204206/151‑1 (MQ=255)
                                     ggTTATGCCATGCCAGCGCAAGAGTTGCTTCCCCGCCATGCGGCCCCCCCAGTCCGGGGATTACCTCTTAAAGGGTGGTTAAGCCTTGCTCTGCGGCCGGGCGGCTGTTTAAAAAAGGCGTACGCATCCCCGCCTGTCGTGCGTGCGttt                                                                                                                 <  1:44109/150‑1 (MQ=255)
                                     ggTTATGCCAGGCCAGCGCCAGGGTTGTTTCCTCGCCAGGCGGCACCCCCAGGACGGGGATTAACTCTTCAAAGGTGGTTAAGCCTGCCTCTACGGCCAGGCGGCCGTTTACAAAAGGCGTCCGCATCCCCGCCTGTCGTGCGTGCGttt                                                                                                                 <  1:4047004/150‑1 (MQ=255)
                                     ggTTATGCAAGGCCAGCGCCAGGGAAGATTACTCGCCCTTTGGTATCCCCAGTACGCGGATTAACTCTTTAAAAGAAATAAAGCCTTGCTCTACGGCCAGGCAGACGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTc                                                                                                                >  1:1515990/1‑151 (MQ=255)
                                         atgccatgccAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTACCTCTTAAAAGGTGTTTAAGCCTTGCTCTACGGCCGGGCAGCCTTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAg                                                                                                             <  1:235458/150‑1 (MQ=255)
                                         atgccatgccAGCGCCAGAGTTGCTTACTCGCCATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCa                                                                                                              <  2:4093477/149‑1 (MQ=255)
                                                        cAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCg                                                                                                >  2:1040281/1‑148 (MQ=255)
                                                               gCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCgg                                                                                      >  1:3416272/1‑151 (MQ=255)
                                                                      cGCAATGCGGCACCCCCAGTACGCGAATTACCTCTTAAAAGGTGGTTAAGCTTTGCTCTACGGCCGGGCGGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATaa                                                                                <  1:3033999/150‑1 (MQ=255)
                                                                        ccATGCGGCACCCCCAGTACGCGGATTAACTCTTCAAGGGTGGTTACGCTTTGCTCTACGGCCAGGCAGCCTTTTAAAAAAAGCGTACGAATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCt                                                                             <  1:2515989/151‑1 (MQ=255)
                                                                                   ccccAGTACGCGGATTAACTCTTCAAAGGTGGTTAACCCTTCCTCTACGGCCGGGCAGCCTTTTTAAAAAAGCGTCCGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGa                                                                    <  2:3245981/149‑1 (MQ=255)
                                                                                           cggggATTAACTCTTCAAAGGTGGTTAAGCTTTGCTCTGCGGCCGGGCGGCTGTTTTAAAAAAGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGtt                                                            <  2:2933481/146‑1 (MQ=255)
                                                                                             cgGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATgg                                                        >  2:90336/1‑151 (MQ=255)
                                                                                              gaaTTACCTCTTAAAAGGTGGTTAAGCTTTGCTCTCCGGCCGGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTAt                                                          <  2:1614202/146‑1 (MQ=255)
                                                                                                           aaaGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCTTTTTAAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCaa                                            <  2:938653/149‑1 (MQ=255)
                                                                                                            aaGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAAt                                          >  2:369840/1‑150 (MQ=255)
                                                                                                                     aaGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTAc                                <  2:2890691/151‑1 (MQ=255)
                                                                                                                          ttGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGAcc                            >  2:616347/1‑150 (MQ=255)
                                                                                                                                  cGGCCAGGCGGCCGTTTTAAAAAAGCGTACCCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATaaa                       <  1:1768700/147‑1 (MQ=255)
                                                                                                                                     ccgggcggcttttttCAAAAGGCGTCCGCAAATCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCg                  <  2:2773082/138‑1 (MQ=255)
                                                                                                                                     ccAGGCAGCCGTTTTAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGt                 <  2:3049495/150‑1 (MQ=255)
                                                                                                                                        ggCAGCTGTTTAAAAAAAGCGTACGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTa             <  1:3206253/151‑1 (MQ=255)
                                                                                                                                        ggCAGCCTTTTACAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTa             <  2:2989353/151‑1 (MQ=255)
                                                                                                                                                 tttAAAAAGGCGTCCGCATCCCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTa     <  2:2274938/150‑1 (MQ=255)
                                                                                                                                                   tCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTacac  >  1:1706407/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCCCATAGCATCCCATCTTGCGCATTGCCGTCGCCGGTTATGCCATGCCAGCGCCAGAGTTGCTTACTCGCCATGCGGCATCCCCAGTACGCGGATTAACTCTTCAAAGGTGGTTAAGCCTTGCTCTACGGCCAGGCAGCCGTTTTCAAAAAGCGTACGCATACCCGCCTGTCGTGCGTGCGTTTCCAGCGATTCAACGTCGGTATTAGCGGAAATAAGCTGACGAATGACCGGCGTTATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACAC  >  NC_000913/115649‑115946

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