Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913307,8210GT48.5% 36.2 / 13.9 33Y504* (TAC→TAA) ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (3/14);  new base T (16/0);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCG  >  NC_000913/307684‑307969
                                                                                                                                         |                                                                                                                                                    
ccTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTg                                                                                                                                          <  2:251868/150‑1 (MQ=255)
          cGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCGTTATTGTAACTTGTGATGGATATGcccc                                                                                                                                >  1:324725/1‑149 (MQ=255)
             tGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTGGTGCTGGATGTGccc                                                                                                                                 >  2:35299/1‑146 (MQ=255)
              gCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTc                                                                                                                           <  2:523329/151‑1 (MQ=255)
                      cGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTTGGAGATGCTGAATGTGTCCCGCTTCGAcc                                                                                                                       >  2:415270/1‑147 (MQ=255)
                      cGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGGCATCATTTTAGCTGATGCTGAAAGTGCCCCGCTTCGAcc                                                                                                                       >  2:231827/1‑147 (MQ=255)
                                  gcTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTTACTGATTCTGAATGTGCCCCGCCTCGACCACAATGGGGTCAgg                                                                                                         >  2:538839/1‑149 (MQ=255)
                                    tGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCCTTTTAGCTGATGGTGAATGTGCCCAGCTTCGACcccaaggaggttagggt                                                                                                       >  2:457897/1‑136 (MQ=255)
                                         aCCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTTTTGCTGAAAGGGCCCCGCCTTGACCCCAGTGAGGTCCGGGTGAgtg                                                                                                  >  2:369570/1‑149 (MQ=255)
                                               ccGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGccgccg                                                                                          <  1:350724/151‑1 (MQ=255)
                                                      ccAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGGCATCATTTTAGCTGATGCTGAATGTGCCCAGCTTCCACCACAAGGGGTTTAGGGTGAGGGTGGCGCCGGCTGCAc                                                                                   >  1:263841/1‑151 (MQ=255)
                                                             aaCCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGGTGATGCTGAATGTGGCCCGCTTCGACCACAGGGAGGTTAGGGTGGGGGTGGCGCCGAGTTCACGGgtgt                                                                             >  2:249749/1‑150 (MQ=255)
                                                              aCCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCg                                                                           <  2:81055/151‑1 (MQ=255)
                                                                     gTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTTGCTGATGCTGAATGTGCCCCGCTTTTCCCACAGTGAGTTCAGGGTTAGGGTGTCGCCGCTTTCACCGTTGTTGGCATCGc                                                                    >  2:406792/1‑151 (MQ=255)
                                                                       aCCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATTATTTTTGCTGATTCTTAATGTGCCCAGCTTCGGCCACCCGGAGATCAAGTTTAGTGTGCCGCCCATTGCACgggttttgggctc                                                                      >  2:441594/1‑136 (MQ=255)
                                                                        ccGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCCTTTTTGCTGATGCTTAATGTGCCCCCCTTCGCCCCCCATGAGTTTCGGTTTAGTGTGCCCCCCCTTTCACCGTTTTCGGCATCGCTAc                                                                 >  2:123093/1‑151 (MQ=255)
                                                                        ccGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTa                                                                  <  2:462792/150‑1 (MQ=255)
                                                                        ccGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTAc                                                                 <  1:457807/151‑1 (MQ=255)
                                                                         cacTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGTCGCCGATTGCACGGTTGTCGGCATCGCTACg                                                                <  1:414084/149‑1 (MQ=255)
                                                                               aGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTTTAGCTGATGCTGAATGTGCCCCGCTTCCACCCCAATGAGGTCAGGGTGAGTGTGGCGCCGATTGCACGGGTTTTGGCATCGCTACggccc                                                            >  2:438335/1‑145 (MQ=255)
                                                                                 aCATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAatt                                                        <  2:344965/151‑1 (MQ=255)
                                                                                     tttGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTTATCATTTTAACTGATGCTGTATGTGGCCCGCCTTGACCACCGTTGGTTTCGGGTGGGGGTGGCGCCGGCTGCACGGGTGGTGGGACCGCCACCTTGCAatttatt                                                    >  1:282718/1‑151 (MQ=255)
                                                                                                 cTGCCGTGTAATAATGGCTGTTGTAACGGCGGGTATTCTTTTAGCCGGTGCTGAAAGTGCCCCCCGTCCACCCCAGGGGGGTTAGGGTGAGGGGGGCGCCGCATGCACGGGTGTCGGCAACCCCACGTCGCGatttattgtccaagtgggg                                        >  1:578787/1‑140 (MQ=255)
                                                                                                          aataatGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACTGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGTGTTTTTTCCtt                               <  1:178894/151‑2 (MQ=255)
                                                                                                           ataatGGCTGTTGTAACGGCGGTCATCATTTTAGCTGATGCTGAATGTGCCCCGCTTCGACCACAGGGAGGTCAGGGTGGGGGTGCCGCCGGTTTCACGGGTGTCGGCATCGCTACGTCGCAAATGATTGCCAATGCGGGgttttt                                   >  2:502444/1‑146 (MQ=255)
                                                                                                              atGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTa                           <  1:542057/151‑1 (MQ=255)
                                                                                                              atGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTa                           <  2:168794/151‑1 (MQ=255)
                                                                                                                ggCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTaa                          <  2:61527/150‑1 (MQ=255)
                                                                                                                  cTGTTGTAACGGCGGTCCTCATTTTAGCTGATGCTGAATGTGCCCCGCTTCGACCACAGTGGGTTCAGGTTGAGTGGGCCGCCGGTTGCACGGTTGTCGGCATCGCTACGTCGCCATTGATTGCCAATGCGGGGTTTTTTCTGATTTAccc                       >  1:366983/1‑151 (MQ=255)
                                                                                                                      tgtAACGGCGGTCAGCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCcaca                    <  1:186667/150‑1 (MQ=255)
                                                                                                                       gtAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCcactg                   <  1:56661/150‑3 (MQ=255)
                                                                                                                                      atTGTAACTGATGCTGAATGTGCCCAGCTTCGACCACAGTGCGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTccc   >  1:518317/1‑151 (MQ=255)
                                                                                                                                       ttGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCg  >  1:238699/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTGATAGTAATCTGCCGTGTAATAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAGTGAGTTCAGGTTGAGTGTGCCGCCGATTGCACGGTTGTCGGCATCGCTACGTCGCAATTGATTGCCAATGCGGGTTTTTTCCTGATTAACCCACAGCGAACTAAAGCCTCCCG  >  NC_000913/307684‑307969

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