Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,366,935	C→T	*326* (TAG→TAA)	pspF ←	psp operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,366,935	0	C	T	100.0%	128.1 / NA	40	*326* (TAG→TAA)	pspF	psp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

TGAGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCT  >  NC_000913/1366805‑1367065                                                                                                                                   |                                                                                                                                   tGAGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTg                                                                                                                            <  1:158730/139‑1 (MQ=255)    gCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTg                                                                                                                         >  2:510715/1‑139 (MQ=255)    gCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTg                                                                                                                         >  2:168323/1‑139 (MQ=255)               cccAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTcaaca                                                                                                              >  2:214280/1‑139 (MQ=255)                      cAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       <  1:168323/139‑1 (MQ=255)                         taatGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGaa                                                                                                    <  1:544840/139‑1 (MQ=255)                         taatGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGaa                                                                                                    <  1:510715/139‑1 (MQ=255)                                 gATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTgatgat                                                                                            >  2:782673/1‑139 (MQ=255)                                       ctgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       >  2:341794/1‑122 (MQ=255)                                       ctgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       <  1:341794/122‑1 (MQ=255)                                       ctgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTa                                                                                      >  2:385796/1‑139 (MQ=255)                                        tgcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTaa                                                                                     >  2:475193/1‑139 (MQ=255)                                         gcgctTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAc                                                                                    <  1:475193/139‑1 (MQ=255)                                              tATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAg                                                                               <  1:782673/139‑1 (MQ=255)                                                       tGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        >  1:397347/1‑89 (MQ=255)                                                       tGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        <  2:397347/89‑1 (MQ=255)                                                         cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAc                                                                                    >  1:68543/1‑123 (MQ=255)                                                         cAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAAc                                                                                    <  2:68543/123‑1 (MQ=255)                                                          aaaCGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGAT‑AGTTAACCCCAGTAATTCAGCcgcg                                                                  <  1:230148/139‑1 (MQ=255)                                                              gacgTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCGGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        >  2:592229/1‑82 (MQ=255)                                                              gacgTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGc                                                                                                                        <  1:592229/82‑1 (MQ=255)                                                              gacgTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       <  2:626159/115‑1 (MQ=255)                                                              gacgTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGtt                                                                                       >  1:626159/1‑115 (MQ=255)                                                                   aTTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcgcg                                                                >  2:712958/1‑133 (MQ=255)                                                                   aTTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCcgcgcg                                                                <  1:712958/133‑1 (MQ=255)                                                                         tatgtatgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTtata                                                    >  1:323258/1‑139 (MQ=255)                                                                           tgtatgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATAtt                                                  >  1:73671/1‑139 (MQ=255)                                                                               tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTccc                                              >  2:82121/1‑139 (MQ=255)                                                                               tgCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTccc                                              <  1:214280/139‑1 (MQ=255)                                                                                     aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       <  2:14368/76‑1 (MQ=255)                                                                                     aTAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAAcgcgcg                                                                                                       >  1:14368/1‑76 (MQ=255)                                                                                      tAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACCACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCACTTGTTGc                                       >  2:170839/1‑139 (MQ=255)                                                                                              tCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGt                                                                                                                          <  2:472681/48‑1 (MQ=255)                                                                                              tCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGt                                                                                                                          >  1:472681/1‑48 (MQ=255)                                                                                                     gcATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTgcagc                        <  2:323258/139‑1 (MQ=255)                                                                                                      cATCCGGCAAGTTGTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTgcagca                       <  2:43284/139‑1 (MQ=255)                                                                                                                   gTATTGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCctg          <  2:73671/139‑1 (MQ=255)                                                                                                                     aTTGCTCAACTTCGTAAAATCGGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTGGCAGCAACTCTTTTTCctgct        <  2:285323/139‑1 (MQ=255)                                                                                                                      ttGCTCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCctgctg       <  2:469683/139‑1 (MQ=255)                                                                                                                          tCAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATc   >  1:24046/1‑139 (MQ=255)                                                                                                                           cAACTTCGTTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCt  >  1:276943/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                   TGAGCGATTAATAGCCCAACGCAATAATGTCTGATGCGCTGCGCTTATCAGGCCTGCAAACGACGTATTGATTATGTATGCCGAATAAAGCATTCACGCCGCATCCGGCAAGTTGTATTGCTCAACTTCGCTAAATCTGGTGCTTTTTCAACAACGCGCGGAACTGATGATAAGTTAACCCCAGTAATTCAGCCGCGCGCTTCTGGTTATATTTCCCTTGTTGCAAACTGAGTTGCAGCAACTCTTTTTCCTGCTGCATCT  >  NC_000913/1366805‑1367065

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