Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,439,141	C→T	*123* (TAG→TAA)	yhdN ←	DUF1992 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,439,141	0	C	T	100.0%	115.3 / NA	38	*123* (TAG→TAA)	yhdN	DUF1992 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/17);  total (21/17)

TAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATC  >  NC_000913/3439014‑3439275                                                                                                                                |                                                                                                                                       taGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGtt                                                                                                                             >  2:987276/1‑139 (MQ=255)     cGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                         <  2:867405/139‑1 (MQ=255)                 tCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCaacaac                                                                                                             <  1:987276/139‑1 (MQ=255)                  cAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACt                                                                                                            >  2:274616/1‑139 (MQ=255)                       cGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCa                                                                                                       >  2:334381/1‑139 (MQ=255)                               aTGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATc                                                                                               >  2:607327/1‑139 (MQ=255)                                      atcatcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGt                                                                                        >  1:715314/1‑139 (MQ=255)                                               tgTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                               <  1:334381/139‑1 (MQ=255)                                               tgTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                               >  1:346564/1‑139 (MQ=255)                                                 tttttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                         <  2:1076674/95‑1 (MQ=255)                                                 tttttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                         >  1:1076674/1‑95 (MQ=255)                                                  ttttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTa                                                                            <  1:677152/139‑1 (MQ=255)                                                              cGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCt                                                                >  1:567297/1‑139 (MQ=255)                                                                       tCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGt                                                                                                                                    >  2:543423/1‑62 (MQ=255)                                                                       tCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGt                                                                                                                                    <  1:543423/62‑1 (MQ=255)                                                                       tCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGc                                                       >  2:51889/1‑139 (MQ=255)                                                                          ggTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTc                                                    >  1:617565/1‑139 (MQ=255)                                                                              ttACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAtt                                                <  1:51889/139‑1 (MQ=255)                                                                                     cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaacaa                                                                                                              >  1:574672/1‑70 (MQ=255)                                                                                     cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaacaa                                                                                                              <  2:574672/70‑1 (MQ=255)                                                                                              tGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCg                                <  2:346564/139‑1 (MQ=255)                                                                                                    cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                          >  1:978864/1‑139 (MQ=255)                                                                                                    cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                          <  1:274616/139‑1 (MQ=255)                                                                                                    cTCACCAATGCGATACATACATACTCCATTACTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                          <  2:567297/139‑1 (MQ=255)                                                                                                         cAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACct                     >  1:955316/1‑139 (MQ=255)                                                                                                            gacagtacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCtt                  <  1:119287/134‑1 (MQ=255)                                                                                                              cGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAAt                                                                                 <  2:910405/74‑1 (MQ=255)                                                                                                              cGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAAt                                                                                 >  1:910405/1‑74 (MQ=255)                                                                                                                atacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGata              <  2:1156502/139‑1 (MQ=255)                                                                                                                 tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCg                                <  1:1210208/120‑1 (MQ=255)                                                                                                                 tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCg                                >  2:1210208/1‑120 (MQ=255)                                                                                                                           cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGa                                      >  2:29949/1‑104 (MQ=255)                                                                                                                           cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGa                                      <  1:29949/104‑1 (MQ=255)                                                                                                                            tCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAGAATCGGTACTCAATCCAGCTTATCGAAGCGTCACGTCCCGTAATGACAATCGACGACAAACCTCTTGATATTGCgtactgac  >  2:119287/1‑132 (MQ=255)                                                                                                                            tCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGatcatc  >  2:155428/1‑139 (MQ=255)                                                                                                                            tCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGatcatc  >  2:205797/1‑139 (MQ=255)                                                                                                                             ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAGTCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                          <  2:107698/114‑1 (MQ=255)                                                                                                                             ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                          >  1:107698/1‑114 (MQ=255)                                                                                                                                |                                                                                                                                       TAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATC  >  NC_000913/3439014‑3439275

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