Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,503	0	G	A	51.6%	15.7 / 40.4	32	intergenic (‑521/+82)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base A (16/0);  total (22/10)
Fisher's exact test for biased strand distribution p-value = 6.77e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-05

ATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGT  >  NC_000913/987363‑987650                                                                                                                                             |                                                                                                                                                    aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAAtttttt                                                                                                                                              >  1:348981/1‑144 (MQ=255)   gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGGAAATAGATTATTTGACtt                                                                                                                                         >  1:356832/1‑151 (MQ=255)         aaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATTCTTTACTTaacaa                                                                                                                                    >  2:316328/1‑150 (MQ=255)           aTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaaca                                                                                                                                 <  2:7339/151‑1 (MQ=255)                aTGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATa                                                                                                                             <  2:53341/150‑1 (MQ=255)                    aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATGCTTGACTTtaaaaaaaaaaaaatt                                                                                                                         >  1:17289/1‑140 (MQ=255)                          tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAAAGAATGTTTGACTTTAAAAAAAGAAAAAttttca                                                                                                                     >  1:602938/1‑144 (MQ=255)                            cTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAAAAATTGCTTTACTTTAAAAACAAAAAAATTgaaagaaa                                                                                                                 >  1:657365/1‑150 (MQ=255)                                tAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTCTTTTATTAAAAAAAAAAAAACTTGAAAGaaaa                                                                                                                >  2:278849/1‑147 (MQ=255)                                     aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                                                                                                        <  2:21479/150‑1 (MQ=255)                                       gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTTCTTGTCTTAACAAACAAATAAAATTACAGaaaaaaaaaaaaaa                                                                                                      >  2:27315/1‑150 (MQ=255)                                         aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAagag                                                                                                   <  2:633015/151‑1 (MQ=255)                                          aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGc                                                                                                  <  2:304425/151‑1 (MQ=255)                                          aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTTGTTTAATTTACAAAAAAATAAATTGACAGaaaaaaaaaaaaaaaa                                                                                                    >  2:339109/1‑147 (MQ=255)                                                 aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAATAAAATGCTTGTATTTACAAAAAAATAAATTTACAGAAAAAAAAAAAAAAAAACCAGCggg                                                                                           >  1:316010/1‑151 (MQ=255)                                                      cAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGTTTTTCTTTAAAAAAAAAAAAATTTAAAGAAAAAAAAAAAAAAAAACCAGCGGgcg                                                                                         >  2:662024/1‑148 (MQ=255)                                                          aaGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAAATTTCTTTTTTTAAAAAAAAAATAAATTTACAAAAAAGAAAAAAAAAAAAGCAGCGGGGGCTgg                                                                                     >  1:31450/1‑148 (MQ=255)                                                              aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTGGCAAAAAAAAATGTTGTATTTACAAAAAAAAAAATTTAAAGAAAAGAAAAAAAAAAAAACAGCgggggggggcgctt                                                                                >  2:86472/1‑138 (MQ=255)                                                              aTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTCAAAAAAAATGAGTGGCTTTAAAAAAAAATAAATTGACAGAAAAAAAAAAAAAAAAGACAGCGGGCGCTCGCTCtt                                                                                >  2:422989/1‑149 (MQ=255)                                                                              gcgcGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGa                                                               <  2:394991/150‑1 (MQ=255)                                                                                   tAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGCTTGATTTTAAAAAAAAAAAAATTGACAGAAAAGAAAAAAAAAAAAACAGCGGGGGCCGGCTCTTGTAGACGCCAAAATAAAAACTgg                                                         >  2:582708/1‑151 (MQ=255)                                                                                     aaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCCAATAAATTGGTTGAATTAACAAAAAAAAAAAATTACAGAAAAAATAAAAAAAAAAACAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAAAAACTgg                                                         >  1:161532/1‑149 (MQ=255)                                                                                            cAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATATATTGTTTTTCTTTACAAAAAAAAAAAATTTCAGAAAAAAAAAAAAAAAAGCCAGCGGGGGCTGGGTCTTTGAAAAGCTAAAATAAAAGATGGCGGTAc                                                   >  1:543111/1‑148 (MQ=255)                                                                                                gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                            <  2:31449/151‑1 (MQ=255)                                                                                                    tatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAAAGATTTCTTGACTTCATAAACAGAATACATGACAGAAAAAAGAAAAAAAAAGCCCACCGGGGCCGGGTCTTCAAGACGCCTAATTAAAAAATGGCGTTACACgggg                                             >  2:141142/1‑144 (MQ=255)                                                                                                            aaCAGGATGATCGTTATGGCGGGGGTGCAAATAAATTTCTTGACTTAACAAAAAAATAACATGACAGAAAAAAAAAAAAAAAAGGCAGCGGGCGCGGGCTCTTGTAGAAGCCAAATTTGAAGCTGGCGGTACGCGGAGGACGTggg                                     >  1:627391/1‑146 (MQ=255)                                                                                                              cAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGCTTGTCTTTACAAAAAGATAAATTTTTAGAAAAAAAAAAAAAAAAAACAGCGGGGGCTGGGTCTTTTAGAAGGTCAAATAAAAAATGGCGGTTAGCGGAGGACGgggggaggg                                >  2:318757/1‑144 (MQ=255)                                                                                                                     atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatc                       <  2:17289/151‑1 (MQ=255)                                                                                                                                      gCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc      <  1:77412/151‑1 (MQ=255)                                                                                                                                        aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc     >  1:490483/1‑150 (MQ=255)                                                                                                                                        aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa    <  2:542979/151‑1 (MQ=255)                                                                                                                                          aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGt  >  1:200145/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                    ATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGT  >  NC_000913/987363‑987650

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