Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	401,893	T→C	G49G (GGT→GGC)	phoA →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	401,893	0	T	C	97.5%	136.4 / NA	40	G49G (GGT→GGC)	phoA	alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/17);  minor base A (1/0);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCACTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGTGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGCCGCACGTAATTATGCCGAA  >  NC_000913/401763‑402010                                                                                                                                   |                                                                                                                      ttGCACTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGa                                                                                                                     >  2:450135/1‑133 (MQ=255)    cACTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCa                                                                                                                  >  1:431423/1‑133 (MQ=255)    cACTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCa                                                                                                                  <  2:463248/133‑1 (MQ=255)      cTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGa                                                                                                                >  2:107559/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGTCGTTTAACGGGCGATCAGACTg                                                                                                             >  1:104441/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTg                                                                                                             >  1:403695/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTg                                                                                                             >  2:102097/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTg                                                                                                             >  2:101090/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCAGTGCTCGACGTATAGCGGGCGGTCAGACTg                                                                                                             >  1:165540/1‑133 (MQ=255)         gCACTCTTACCGTTACTGTCTACCCCTGTGACAAAACCCCGGACACAAAAAACGCCTGTTCTGGAAAACCGGGCTGTCCTTGGCGATATTAATGCACCCGACGGACCCACCCGTTTAGCGGGTTATCAAAcgg                                                                                                             >  2:248688/1‑131 (MQ=255)            ctctTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCg                                                                                                          <  1:107559/133‑1 (MQ=255)              ctTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGct                                                                                                        >  1:60655/1‑133 (MQ=255)                 aCCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGctct                                                                                                      >  1:37931/1‑132 (MQ=255)                 aCCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGctct                                                                                                      <  2:37931/132‑1 (MQ=255)                 aCCGTTACAGTTTACGCCCGTGCCAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTg                                                                                                     <  1:172416/133‑1 (MQ=255)                 aCCGTTACAGTTTACCCCCGTGCCAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTg                                                                                                     <  1:133920/133‑1 (MQ=255)                       aCTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTCCTCGCCGTATAAAGGGGGATCAGACAGCCGTTCTGCGTGAt                                                                                               >  2:172416/1‑133 (MQ=255)                       aCTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTACTCGCCGTGTAAGGGGAGATCAGACAGCCGATCTGCGTGAt                                                                                               >  2:133920/1‑133 (MQ=255)                         tGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATtc                                                                                             >  1:342975/1‑133 (MQ=255)                           tttACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATtctc                                                                                           >  1:81927/1‑133 (MQ=255)                            ttACCCCTATGACAAAAGCCCGCACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTctct                                                                                          <  1:248974/133‑1 (MQ=255)                             tACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCtt                                                                                         >  2:410500/1‑133 (MQ=255)                             tACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCtt                                                                                         <  1:373270/133‑1 (MQ=255)                                cccTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGc                                                                                      >  1:11892/1‑133 (MQ=255)                                  cTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGa                                                                                    >  2:466903/1‑133 (MQ=255)                                   tgtgACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGAt                                                                                   >  1:211108/1‑133 (MQ=255)                                     tgACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATaa                                                                                 >  1:504467/1‑133 (MQ=255)                                     tgACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATaa                                                                                 >  2:435952/1‑133 (MQ=255)                                                           aaTGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTg                                                           <  1:450135/133‑1 (MQ=255)                                                               ccTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGa                                                       >  2:282817/1‑133 (MQ=255)                                                                cTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGAt                                                      <  2:60655/133‑1 (MQ=255)                                                                       ggAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCgat                                               <  2:11892/133‑1 (MQ=255)                                                                        gAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCgatg                                              >  1:87613/1‑133 (MQ=255)                                                                                   tgctCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGAc                                   <  2:504467/133‑1 (MQ=255)                                                                                        aGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGa                              <  1:282817/133‑1 (MQ=255)                                                                                             gATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTa                         <  2:374323/133‑1 (MQ=255)                                                                                               tatTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGTACTCGGAAATTACt                       >  1:27416/1‑133 (MQ=255)                                                                                                 ttACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGc                     <  2:59333/133‑1 (MQ=255)                                                                                                  tACTGCCCCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGAGTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGCGACTCGGAAATTACTGcc                    <  1:248688/133‑1 (MQ=255)                                                                                                  tACTGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGcc                    <  2:403695/133‑1 (MQ=255)                                                                                                     tGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGCCGCa                 <  1:410500/133‑1 (MQ=255)                                                                                                     tGCACCCGGCGGTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGCCGCa                 <  1:102097/133‑1 (MQ=255)                                                                                                            ggcggTGCTCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGACGATGGGAGGGGGGACTCGGAAATTACTGCCGCACGTAAtt          >  2:212765/1‑133 (MQ=255)                                                                                                                    tCGCCGTTTAACGGGCGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGGTGGGGGACTCGGAAATAACTGCCGCGCGTAATTATGCCGaa  >  2:169844/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                      TTGCACTGGCACTCTTACCGTTACTGTTTACCCCTGTGACAAAAGCCCGGACACCAGAAATGCCTGTTCTGGAAAACCGGGCTGCTCAGGGCGATATTACTGCACCCGGCGGTGCTCGCCGTTTAACGGGTGATCAGACTGCCGCTCTGCGTGATTCTCTTAGCGATAAACCTGCAAAAAATATTATTTTGCTGATTGGCGATGGGATGGGGGACTCGGAAATTACTGCCGCACGTAATTATGCCGAA  >  NC_000913/401763‑402010

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