Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	403,368	C→T	Q30* (CAA→TAA)	psiF →	phosphate starvation‑inducible protein PsiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	403,368	0	C	T	100.0%	65.7 / NA	24	Q30* (CAA→TAA)	psiF	phosphate starvation‑inducible protein PsiF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/12);  total (12/12)

ATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAA  >  NC_000913/403159‑403591                                                                                                                                                                                                                  |                                                                                                                                                                                                                                aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCACCCCGACTTACTCGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATACCATTACTGGTTACCTTGCTTTTCGGACTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACATCCTGTAATCAGCAGGCGACGgcgc                                                                                                                                                                                          >  1:494485/1‑249 (MQ=255)                                       ggtTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACAATCATGGTTAAAACATGAAGAGGGATGGTGCCATGAAAATAAAATTACCGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCCCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACACTAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg                                                                                                                                                   >  1:81266/1‑250 (MQ=255)                                                gCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAgg                                                                                                                                                                                       >  1:119432/1‑205 (MQ=255)                                                        cAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTGTTTAACCCCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGGATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGc                                                                                                                                              <  1:417941/238‑1 (MQ=255)                                                             aGNCTTNANNGNNGATNGCGGGCGCATACGCTCATGGTNNAAACATGAANAGGTATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTNTTTTTAACCNNCGNCGGCGNTNNNGNNNGANCTTTAACCCCACAATAACANCGNNTGACCTNCTGTAATCAGCAGGCGACGGCGCAgg                                                                                                                                                                                       <  2:119432/192‑1 (MQ=255)                                                                      ggCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTgg                                                                                                                    <  1:403708/250‑1 (MQ=255)                                                                                 ggCGCATACGCTAATGGTTAAAACATGAAGAGGGCTGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGa                                                                                                                                                                          >  1:427647/1‑185 (MQ=255)                                                                                         ngCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAg                                                                                                                               >  1:314712/2‑220 (MQ=255)                                                                                          ncTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATca                                                                                                                                                                                                        >  1:303771/2‑146 (MQ=255)                                                                                             cATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAggg                                                                                                                                                                            >  1:211194/1‑171 (MQ=255)                                                                                                    aanaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATca                                                                                                                                                                                                        <  2:303771/136‑1 (MQ=255)                                                                                                                        gCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGCGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGACGAACAGCAAGTCTGCGCCTGGCGaaaaaa                                                                                                            >  1:9669/1‑208 (MQ=255)                                                                                                                          tATGAAAATAACATTACGGGTTACCGGGCTTTTCGGTCTGGGGTGTTTAACCACCGTCTGCGTTGCCGAGAGAACTGGAACCCCACAATAACAGCGTAGGACCTGCTGTAATCAGCAGGCGACGGCGCAg                                                                                                                                                                                        <  1:332127/130‑1 (MQ=255)                                                                                                                                                  cTTGCTTTTCTGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTAGGACCTCTTGTAATCAGCAGGCGACGGCGCATTCGGTTAAAGGGGGTGCTCGGAAGTCCGACATGAGTGATTGTCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTg                                                                                                       <  1:330132/187‑1 (MQ=255)                                                                                                                                                       ttttCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGAcctcc                                                                                                                                                                                                                >  2:20566/1‑77 (MQ=255)                                                                                                                                                       ttttCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGAcctcc                                                                                                                                                                                                                <  1:20566/77‑1 (MQ=255)                                                                                                                                                               nnnnnnttnnnNNCCANCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAg                                                                                                                               <  2:314712/144‑1 (MQ=255)                                                                                                                                                                        ttAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTg                                                                                                       <  1:172410/165‑1 (MQ=255)                                                                                                                                                                                 gTCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCATGCGTTGAAAgg                                                                                                                                                                             <  1:414597/86‑1 (MQ=255)                                                                                                                                                                                  tCGGCGCTGCCGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTg                                                                                                       <  1:427176/155‑1 (MQ=255)                                                                                                                                                                                            cGAGAGAACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGaa  >  1:240532/1‑246 (MQ=255)                                                                                                                                                                                                 gaACTTTAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGAGTTGAAAGGGGATGCTCGTAGGACCTACCTGCGTCATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCACAGCAGCAAAAGAt                                                                                   >  1:328280/1‑160 (MQ=255)                                                                                                                                                                                                 gaACTTTAACCACACAATAACAGCGTATGACCTCCTGTAATCAGCAgg                                                                                                                                                                                                   >  1:382285/1‑48 (MQ=255)                                                                                                                                                                                                      ttAACCCCACAATAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCGCAGCAGCAAAAGAt                                                                                   <  1:171624/155‑1 (MQ=255)                                                                                                                                                                                                             nacaataacaGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGtgcctg           >  1:300125/2‑220 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                ATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGCAAGTCTGCGCCTGGCGAAAAAAGTTTGACGCCACAGCAGCAAAAGATGCGCGAATGCAATAATCAAGCAACACAACAATCTCTGAAAGGTGATGATCGTAATAAGTTTATGAGTGCCTGCCTCAAGAA  >  NC_000913/403159‑403591

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