Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,494	0	T	G	41.2%	82.5 / 12.2	51	intergenic (‑512/+91)	ompF/asnS	outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (5/25);  new base G (21/0);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 6.85e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  NC_000913/987349‑987637                                                                                                                                                  |                                                                                                                                                 ttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTgta                                                                                                                                              >  1:71734/1‑148 (MQ=255)     gACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAg                                                                                                                                         <  1:511786/151‑1 (MQ=255)          ccAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTgggaaaaaa                                                                                                                                        >  2:728752/1‑139 (MQ=255)            aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑tt                                                                                                                                 <  2:60921/151‑1 (MQ=255)            aGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑tt                                                                                                                                 <  1:674064/151‑1 (MQ=255)                 gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACtt                                                                                                                            <  2:224834/151‑1 (MQ=255)                 gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACtt                                                                                                                            <  2:595075/151‑1 (MQ=255)                    aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATaaaattttttgctt                                                                                                                             >  1:505392/1‑135 (MQ=255)                     aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAGAATGC‑TTTACTTaa                                                                                                                          >  2:70092/1‑149 (MQ=255)                        aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACtt                                                                                                                            <  1:491719/144‑1 (MQ=255)                        aaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAATAGaatttttttctttta                                                                                                                          >  2:491771/1‑133 (MQ=255)                         aTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTaacaaaca                                                                                                                    <  1:641260/151‑1 (MQ=255)                                 aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACtt                                                                                                            <  2:387150/151‑1 (MQ=255)                                    aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAAATGT‑TTTTCTTAACAAAAAAATaaaat                                                                                                            >  2:258834/1‑145 (MQ=255)                                       gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAg                                                                                                       <  2:283656/150‑1 (MQ=255)                                       gtgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGGAAAAAATTGC‑TTTGATTAAAAAAAAGaataat                                                                                                             >  1:630619/1‑139 (MQ=255)                                        tgCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGAGGGGGAAATAAAATGC‑TTTATTTTATAAAAAAAAAACTTTAc                                                                                                         >  2:49968/1‑147 (MQ=255)                                         gCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGTGTTGCAAATAGATTGC‑TTGACTTAACACACAGATAACTTGACAGaaa                                                                                                    <  1:267635/151‑1 (MQ=255)                                            aTTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGa                                                                                                 <  1:741625/151‑1 (MQ=255)                                               aGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGAAATAAATTGA‑TTTACTTTAAAAAAAAAAAAATTGAGCGAAAAGATa                                                                                               >  2:690573/1‑150 (MQ=255)                                                 cTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTGC‑TTGACTTTAAAAACAGATAAATTGACAGaaaaaaaaaa                                                                                             >  2:177783/1‑150 (MQ=255)                                                   aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaa                                                                                          <  2:565834/151‑1 (MQ=255)                                                       aaaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAATTTT‑TTTACTTTAAAAACAGATAAATTTACAGAAAAAATaaaaaaaaa                                                                                       >  2:132637/1‑148 (MQ=255)                                                          aaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCa                                                                                   <  1:36225/151‑1 (MQ=255)                                                           aCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGAATGCCAAAAAAATTGA‑GTTACTTAACAAAAACAAAAATTGACAGGAAAGAAAAAAAAaaaggc                                                                                    >  2:335455/1‑147 (MQ=255)                                                              aaTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc                                                                                 <  2:290823/149‑1 (MQ=255)                                                               aTCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATTT‑TTGTCTTTAAAAAAAAAAAAAATGACaaaaaaaaaaaaaaaaaaacaaccgg                                                                               >  1:396217/1‑140 (MQ=255)                                                                tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGGAAATAAAATTC‑TTGTCTTTAAAAAAAAATAAATTTACAGAAAAAAAAAAAAAAAAGACAGCggg                                                                              >  1:10803/1‑150 (MQ=255)                                                                 ccGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcg                                                                            <  1:487113/151‑1 (MQ=255)                                                                      aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                       <  2:230728/151‑1 (MQ=255)                                                                      aTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAATTGC‑TTGTCTTTACAAAAAAAAAAATTTACAGAAAAAATAAAAAAAAAAACAGCGGGCGCTgg                                                                        >  2:185887/1‑150 (MQ=255)                                                                              gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                               <  2:235423/151‑1 (MQ=255)                                                                                ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGAc                                                             <  1:343364/151‑1 (MQ=255)                                                                                 gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTTGAAAAAGATTGC‑TTTAATTAAAAAAAAAAAAACTTTAAAGAAAAAATAAAAAAAAAAACAGCGGGGGCTggggccttta                                                                >  2:687977/1‑139 (MQ=255)                                                                                   ttAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGAAAAAAAATGT‑TTTTCTTTACAAAAAAAAAAATTGACAGAAAAAAAAAAAAAAAAAAAGACGGGGGGTGGCTCTTGaaaa                                                              >  1:29191/1‑145 (MQ=255)                                                                                     aGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaa                                                        <  1:453323/151‑1 (MQ=255)                                                                                         tttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAAATGC‑TTTACTTTAAAAAAAGAAAAATTTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCCGGCTCTTTAAGAAGCCaaa                                                       >  2:303464/1‑148 (MQ=255)                                                                                          ttGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAATAAATTGCTTTGATTTAA‑AAAAAAAAAAAATGACAGAAAAAAAAAAAAAAAAAACAGCGGGCGCTGGCTCTTTAAAACGCTaaa                                                       >  1:285405/1‑147 (MQ=255)                                                                                             cgcgTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGGGGAAAAAAAATTC‑TTGTCTTTACAAAAAAAAAAAATGATAGAAAAAAAAAAAAAAAAAACAGCGGGGGCCGGGTCTTTGAGAAGCCAAATTaaaa                                                 >  2:207354/1‑150 (MQ=255)                                                                                                      gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                        <  2:10831/150‑1 (MQ=255)                                                                                                      gTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATGG‑TTGACTTAACAAAAAAAAAAATTTACAGAAAAAAAAAAAAAAAAACCAACGGGGGCCGGGGCTTTAAAAAGCtaaaataaaaacagggggt                                        >  2:702619/1‑141 (MQ=255)                                                                                                          cAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAAAATGG‑TTGGCTTAACAAACAAATAAAATGAAACAAAAGATAAAAAAAAAAACCGCGGGCGGTGGCTCTTTAAGAAGCTAAATTAGAAGATTGCGTTAcg                                     >  2:29003/1‑149 (MQ=255)                                                                                                                aaTATTCCGTAACAGGATGATCGTTATGGCGGGGGTGGAAAAAAATTGT‑TTTTATTAAAAAAAAAAAAAATATAGAGAAAAAAAAAAAAAAAAAACAGCGGGGGGCGGCTCTGTTAGACACtaaaataaaaaaaggagggacg                                     >  1:370846/1‑131 (MQ=255)                                                                                                                aaTATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAAAAAAATTT‑TTTAATTTAAAAAAAGATAAAATGGAAGAAAAAAAAAAAAAAAAAACACAGGGGGCTGGGTCTTGAAAAAGCTAAAATCAAAAAAGGCGTTTCGCGCAGTa                              >  1:193588/1‑150 (MQ=255)                                                                                                                   atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                          <  2:776854/151‑1 (MQ=255)                                                                                                                       cGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                      <  2:237024/151‑1 (MQ=255)                                                                                                                                         aTGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgta     <  2:51292/150‑1 (MQ=255)                                                                                                                                         aTGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgta     <  2:282200/150‑1 (MQ=255)                                                                                                                                         aTGGCGGGGTTGCAAAAAGATTTC‑TTGACTTAACAAACAAATAACTTTGCAGAAAAAATAAAAAAAAAGCCAGCGGGCGCTGGCTCTTGAAGAAGCTAAAATAAAAGCTGGCGGTACGCGGAGTACGTGGGGAAGGAATCACATCacgta     >  1:531907/1‑150 (MQ=255)                                                                                                                                           ggCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacgt  >  2:126389/1‑151 (MQ=255)                                                                                                                                           ggCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacg   <  1:132604/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGC‑TTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGT  >  NC_000913/987349‑987637

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