Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,501	0	T	A	45.2%	33.4 / 13.4	31	intergenic (‑519/+84)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (4/13);  new base A (14/0);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.91e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACG  >  NC_000913/987367‑987632                                                                                                                                        |                                                                                                                                    cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTGGGAAAAAAAATGTTTTACttt                                                                                                                      >  1:446684/1‑149 (MQ=255)    aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAAAAAAATGTTTTACTTaaaa                                                                                                                   >  1:482778/1‑148 (MQ=255)      aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTaacaaac                                                                                                                <  2:235913/151‑1 (MQ=255)      aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTaacaaac                                                                                                                <  1:566176/151‑1 (MQ=255)                  aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAAAAAAATTGTTTAATTTAA‑AAACAAATAAAtt                                                                                                       >  1:522963/1‑148 (MQ=255)                     gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAAATTCTTTTATTAACAAAAAAATAAATTGAc                                                                                                    >  2:399805/1‑148 (MQ=255)                       gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaa                                                                                                <  2:529335/150‑1 (MQ=255)                              gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaa                                                                                        <  1:602283/151‑1 (MQ=255)                                tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTTCTTTTCTTTACAAACAAATAACTTTACAGAAAAAATaaaa                                                                                       >  2:33689/1‑150 (MQ=255)                                   gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAAATTCTTTACTTTAAAAAAAAATAAATTTACAGAAAAGAaaaaaaaa                                                                                    >  2:477165/1‑150 (MQ=255)                                    aaaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAaga                                                                                  <  1:250576/151‑1 (MQ=255)                                              tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGGAATAGATTGTTTGTCTTAACAAACAAATAACTTTTAAGAAAAGAAAAAAAAAAAGCCCACggg                                                                         >  2:662406/1‑150 (MQ=255)                                               ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                       <  1:616234/151‑1 (MQ=255)                                                cGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAAATGTTTGACTTTAAAAACAAAAAAAATTACAGAAAAGAAAAAAAAAAAACCAGCGGgcg                                                                       >  1:518407/1‑150 (MQ=255)                                                     tAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTTTAAAAATAATTTTTTTCTTTACAAAAAAAAAAATTTAACGAAAAAAAAAAAAAAAAAACCGCGGGCGCtcg                                                                   >  2:509982/1‑147 (MQ=255)                                                          aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                            <  1:109852/151‑1 (MQ=255)                                                           tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAGAATGCTTGACTTAACAAAAAAAAAAATTTACAGAAAAAATAAAAAAAAAAACAGCGGGGGCTGGGTCtt                                                              >  2:553914/1‑148 (MQ=255)                                                           tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGCAAAAAAAATGCTTTTCTTTAAAAACAAATAAAATGACAGAAAAGAAAAAAAAAAAAACAGCGGGCGCTGGGTCTTgt                                                            >  1:159547/1‑149 (MQ=255)                                                            gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                           <  1:207297/150‑1 (MQ=255)                                                              ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAATTGTTTTTCTTTAAAAACAAAAAAATTTACAGAAAAGATAAAAAAAAAAGCAGCGGGGGCTGGCTCCTGaa                                                           >  2:438718/1‑148 (MQ=255)                                                                       tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                <  2:323439/150‑1 (MQ=255)                                                                         tGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGTAAAAAAAATTTTTGTCTTTAAAAACAAGAAAAATGAAAGAAAAAATAAAAAAAAAAGCAGCGGGGGCCGGGCCTTTTAGAAGCTaaaat                                                >  2:225438/1‑146 (MQ=255)                                                                           cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGGGGAAAAAAAATATTTTATTTAACAAAAAAATAAATTGGTCGGAAAAAAAAAAAAAAAAACAACGGGGGCTGGGGCTTGAAGAAGACAAAATAGaa                                            >  1:666821/1‑150 (MQ=255)                                                                             cgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCt                                         <  1:460010/151‑1 (MQ=255)                                                                                aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAAAATTCTTTAGTTTACAAAAAAATAAATTGACAGAAAAAAAAAAAAAAAAGACAGCGGGCGCCGGGGCCTGAAGAAGcccaaataaaagatg                                        >  2:204941/1‑135 (MQ=255)                                                                                  aTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGt                                    <  2:468040/151‑1 (MQ=255)                                                                                                 atTCCGTAACAGGATGATCGTTATGGC‑GGGGGTGGAAAAAGATTGCTTTACTTAAAAAAAAAAAAAATTTACAGAAAAGAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTGTAGACGCTAAAATAAAAAATGGCGTTACGCGGAGGACGTg                     >  2:367094/1‑151 (MQ=255)                                                                                                       tAACAGGATGATCGTTATGGC‑GGGGGTGTAAATAGATTTCATGACTTTACAAAAAAAAAAAATGTCAGAAAAAAAAAAAAAAAAACCCGCCGGCGCTGGCTCCTTTAGAACCCCAAATAGAAACCGGCGTTAAGCggcggaaatggg                   >  2:222142/1‑137 (MQ=255)                                                                                                          cAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa            <  2:31593/151‑1 (MQ=255)                                                                                                          cAGGATGATCGTTATGGC‑GGGGGTTTAAATAGATATCTTTTTTTAACAAACAGATAAATTTGCAGAAAAAAAAAAAAAAAAAACAGAGGGGGCTGGCTCTTGAAGAAACTAAAATATAAAGTGCCGGTAAGCCGAGAACGTggg                   >  1:391867/1‑144 (MQ=255)                                                                                                                    gTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacg  <  1:553909/151‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    CAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGC‑GGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACG  >  NC_000913/987367‑987632

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