Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5510GA48.3% 8.7 / 22.7 30intergenic (‑569/+34)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (6/9);  new base A (14/0);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 7.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACC  >  NC_000913/987405‑987689
                                                                                                                                                   |                                                                                                                                          
aaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGcc                                                                                                                                        <  1:300723/151‑1 (MQ=255)
 aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAATAAATTTCTTTATTTTAAAAAAAGATAACTTTACAGAAAAGAAAAAAAAAAAGCCa                                                                                                                                       >  1:358988/1‑151 (MQ=255)
    tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCg                                                                                                                                    <  1:520038/151‑1 (MQ=255)
           cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGGAAATAGAATGCTTGACTTAAAAAAAAGATAAATTGACAGAAAAAAAAAAAAAAAAGGCAGCGGgggccg                                                                                                                             >  2:595347/1‑149 (MQ=255)
           cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAATAAATTGCTTGACTTTACAAAAAGATAACTTTACAGAAAAGATAAAAAAAAAAACAGCGGgcgccg                                                                                                                             >  2:363223/1‑149 (MQ=255)
                 gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAGATTGCTTTTCTTTACAAAAAAATAACTTGTCAGAAAAAATAAAAAAAAAACCAGCGGGCGCCGGCTCtt                                                                                                                       >  1:326638/1‑151 (MQ=255)
                    tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                                                    <  1:505567/151‑1 (MQ=255)
                       ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAGATTGCTTGACTTTAAAAACAAAAAAATTTACAGAAAAGAAAAAAAAAAAACCAGCGGGGGCTGGCTCTTGAAGaa                                                                                                                 >  2:127996/1‑150 (MQ=255)
                        gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                                                                                 <  2:9440/150‑1 (MQ=255)
                        gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTTTTTTTCTTTACAAAAAGATAAAATTACAGAAAAGAAAAAAAAAAAGACAGCGGGCGCTGGCTCTTTAAGaag                                                                                                                >  1:74774/1‑149 (MQ=255)
                                    cgcgTAAATGTTGCAATGTGATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                                                    <  1:87674/151‑1 (MQ=255)
                                      cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAAATTCTTTACTTAACAAAAAGATAAATTGGCAAAAAAAAAAAAAAAAAAAGCAGCGGGGGCTGGCTCTTGGAGACGCtaaaataaaaaca                                                                                                  >  1:230356/1‑148 (MQ=255)
                                                   aTGTAATATTCCGTAACAGGCTGCTCGTTATGGC‑GGGGGTTTTTAAAGATTTACTCTCCTCATACACAGAATACTTAACAGAAAACAACAAAAAACAGGCACGGGGGGCGGGCTCTTGACTAAGCTTAATTAAAAGCCGGGGGTAGGCGGa                                                                                     >  2:636439/1‑151 (MQ=255)
                                                   aTGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAATAAATTGCTTGACTTAACAAACAGATAAATTTACAGAAAAGATAAAAAAAAAACCAGCGGGCGCTGGCTCTTTAAGACGCTAAATTAGAAGGTGGCGTTAAGCGGa                                                                                     >  2:395303/1‑151 (MQ=255)
                                                     gTAATATTCCGTAACAGGATGATCGTTATGGC‑GGTGATGCGAAAAAATTGATTTGATTAAAAAACAAATGAATTGAAAGGAAAGATAAAAAAAAAGACAGCGGGCGCGGGCTCTTAAAAACGGTACAATATAAACAGGCGTTACGCGGAgg                                                                                   >  2:154380/1‑150 (MQ=255)
                                                     gTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTTCTTGACTTAACAAAAAAATAACTTTACAGAAAAAATAAAAAAAAAGCCAGCGGGCGCTGGCTCTTGAAGAAGCTAAATTAGAAAATGGCGGTAAGCGGAgg                                                                                   >  2:840743/1‑150 (MQ=255)
                                                     gTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAAATTCTTTACTTTAAAAACAGATAACTTTACAGAAAAGATAAAAAAAAAGCCAGCGGGGGCTGGGTCTTGAAGACGCTAAATTAAAAAGTGGCGTTACGCGGAgg                                                                                   >  2:853957/1‑150 (MQ=255)
                                                          atTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTGTTTTATTTAA‑AAAAAGAAAAATTTTAAGAAAAGATAAAAAAAAAAACAGCGGGCGCGGGCGCCTGTAGAAGCCAAAAAAAAAAAAGGCGGTACGCGGGGGACgcg                                                                              >  1:87099/1‑149 (MQ=255)
                                                           ttCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAATAGATTTGTTGTCTTAACAAAAAGATAACTTTAAAGAAAAGATAAAAAAAAAGACCGAGGGCGCTGGCGCTTTAAAAAGCCAAAATAAAAACTGGCGGTACGCGGGGGACgggg                                                                             >  1:593005/1‑151 (MQ=255)
                                                           ttCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAATAAATTGCTTTACTTAAAAAAAAAAAAACTTTACAGAAAAGATAAAAAAAAAGCCAGCGGGGGCTGGCTCTTGAAGACGCTAAATTAAAAACTGGCGGTACGCGGAGTACGTgg                                                                             >  2:400457/1‑151 (MQ=255)
                                                                          atCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                                               <  2:342958/150‑1 (MQ=255)
                                                                                         gTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                                                >  2:466157/1‑151 (MQ=255)
                                                                                             cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                             >  2:788739/1‑150 (MQ=255)
                                                                                                     ttGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGt                                    >  2:822645/1‑151 (MQ=255)
                                                                                                             cTTAACAAACAGATAACTTGACAGAAAAGATAAAAACAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                             >  2:585093/1‑150 (MQ=255)
                                                                                                                aacaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGa                         <  1:350446/151‑1 (MQ=255)
                                                                                                                  caaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACg                       <  1:324903/151‑1 (MQ=255)
                                                                                                                                 aCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTCGAAGCCCCTAAATTAGAAGCTGGCGTTACGCGGAGTCCGTGGGAACGGAATCACATCACGTACGCTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCACAACCCAGAc         >  1:706763/1‑150 (MQ=255)
                                                                                                                                       aaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAAcc  <  2:136527/151‑1 (MQ=255)
                                                                                                                                         gATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAAc   >  1:134854/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
AAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACC  >  NC_000913/987405‑987689

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