Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,436,833	C→A	A21S (GCA→TCA)	ompN ←	outer membrane porin N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,436,833	0	C	A	95.7%	158.7 / ‑0.4	47	A21S (GCA→TCA)	ompN	outer membrane porin N
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (18/27);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAG  >  NC_000913/1436700‑1436966                                                                                                                                      |                                                                                                                                      gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTgc                                                                                                                                 <  1:663830/140‑1 (MQ=255)        cgcgaaaaCCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCt                                                                                                                          >  1:457659/5‑140 (MQ=255)          ttAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTgc                                                                                                                        >  2:686327/1‑140 (MQ=255)          ttAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTgc                                                                                                                        >  1:564888/1‑140 (MQ=255)           tAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTgcg                                                                                                                       <  1:684185/140‑1 (MQ=255)           tAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTgcg                                                                                                                       <  2:428740/140‑1 (MQ=255)               aCCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGa                                                                                                                   >  1:280143/1‑140 (MQ=255)                cccAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAg                                                                                                                  <  1:403168/140‑1 (MQ=255)                cccAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAg                                                                                                                  <  2:222177/140‑1 (MQ=255)                cccAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAg                                                                                                                  <  1:369758/140‑1 (MQ=255)                cccAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAg                                                                                                                  <  1:287417/140‑1 (MQ=255)                    gACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAgcagg                                                                                                              >  1:544375/1‑140 (MQ=255)                       gcgcATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAgcagggca                                                                                                           <  1:572507/140‑1 (MQ=255)                        cgcATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAgcagggcag                                                                                                          >  1:277073/1‑140 (MQ=255)                           aTAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGaa                                                                                                       <  1:390642/140‑1 (MQ=255)                                tctGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTaaa                                                                                                  <  2:549609/140‑1 (MQ=255)                                tctGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTaaa                                                                                                  <  1:234827/140‑1 (MQ=255)                                tctGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTaaa                                                                                                  <  1:714216/140‑1 (MQ=255)                                     tCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGc                                                                                             <  2:163186/140‑1 (MQ=255)                                              ttCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACttt                                                                                    >  1:143181/1‑140 (MQ=255)                                                 gCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCt                                                                                 <  2:615437/140‑1 (MQ=255)                                                       ttATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCAt                                                                           <  2:457659/140‑1 (MQ=255)                                                          tCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGa                                                                        >  2:592744/1‑140 (MQ=255)                                                              aaaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATaa                                                                    <  1:13090/140‑1 (MQ=255)                                                                aaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAAt                                                                  <  2:277073/140‑1 (MQ=255)                                                                  ataatGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATcc                                                                >  1:482299/1‑140 (MQ=255)                                                                    aatGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCtt                                                              <  2:143181/140‑1 (MQ=255)                                                                      tGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTa                                                            >  2:192542/1‑140 (MQ=255)                                                                                   cTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTg                                               <  2:390536/140‑1 (MQ=255)                                                                                        ccATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCttt                                          <  2:295511/140‑1 (MQ=255)                                                                                             cAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATcc                                     <  2:350373/140‑1 (MQ=255)                                                                                              aGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATccc                                    >  2:545644/1‑140 (MQ=255)                                                                                                     aTTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCg                             <  2:564888/140‑1 (MQ=255)                                                                                                     aTTTGTTGCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCg                             >  2:392453/1‑140 (MQ=255)                                                                                                       ttgttgCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCggg                           <  2:706033/140‑1 (MQ=255)                                                                                                       ttgttgCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCggg                           <  1:384547/140‑1 (MQ=255)                                                                                                       ttgttgCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCggg                           <  2:280143/140‑1 (MQ=255)                                                                                                          ttgCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAg                        >  2:431364/1‑140 (MQ=255)                                                                                                          ttgCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAg                        >  1:285589/1‑140 (MQ=255)                                                                                                            gCCGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTg                      >  1:484137/1‑140 (MQ=255)                                                                                                              cGTCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGaa                    <  2:285589/140‑1 (MQ=255)                                                                                                                tCTTTATTATAAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGaaaa                  <  1:686327/140‑1 (MQ=255)                                                                                                                        ataAACTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTaa          <  2:484137/140‑1 (MQ=255)                                                                                                                            aCTTCGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGga      >  2:324976/1‑140 (MQ=255)                                                                                                                                cGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTGGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAg  >  1:9764/1‑140 (MQ=255)                                                                                                                                cGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAg  >  2:539178/1‑140 (MQ=255)                                                                                                                                cGGCTGAATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGGGCCAGTACTTTGCTTTTCATTGAACAAATCCTTTAGTTATTTTATTTACCTTTTATCCCCATATCGGGGAGTGAAAACATCCTAAAGgagatg  >  2:524728/1‑138 (MQ=255)                                                                                                                                      |                                                                                                                                      GTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCATCTTTCGCTGAATTATCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGTGCAGCACCTGCGGCGAGCAGGGCAGGAATTAAAAGTGCCAGTACTTTGCTTTTCATTGAATAAATCCTTTAGTTATTTTATTTGCCTTTTATCCCCATTTCGGGGAGTGAAAACATCCTAAAGGAGAAG  >  NC_000913/1436700‑1436966

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