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%	168.8 / NA	52	*123* (TAG→TAA)	yhdN	DUF1992 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/24);  total (28/24)

TCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCG  >  NC_000913/3439006‑3439276                                                                                                                                        |                                                                                                                                        tCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAg                                                                                                                                      >  1:656355/1‑139 (MQ=255) tCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAg                                                                                                                                      >  1:327630/1‑139 (MQ=255)        ataGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATTGt                                                                                                                               <  1:1070700/139‑1 (MQ=255)                 aCCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGt                                                                                                                      >  1:276638/1‑139 (MQ=255)                 aCCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGt                                                                                                                      >  1:170509/1‑139 (MQ=255)                      ccTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaac                                                                                                                 <  2:170509/139‑1 (MQ=255)                        ttCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaacaa                                                                                                               <  2:656355/139‑1 (MQ=255)                          cAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACt                                                                                                             <  2:956416/139‑1 (MQ=255)                                    ttCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaac                                                                                                                 <  2:514854/125‑1 (MQ=255)                                    ttCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaac                                                                                                                 >  1:514854/1‑125 (MQ=255)                                              atcatcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCg                                                                                                                                <  1:1005641/100‑1 (MQ=255)                                              atcatcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCg                                                                                                                                >  2:1005641/1‑100 (MQ=255)                                              atcatcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGt                                                                                         >  1:348697/1‑139 (MQ=255)                                                 atcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCGTAATCGCCGCGTaaa                                                                                      <  2:404406/139‑1 (MQ=255)                                                          ttttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTa                                                                             <  2:848058/139‑1 (MQ=255)                                                                          tCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTc                                                             >  1:188876/1‑139 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCa                                                          <  2:276638/139‑1 (MQ=255)                                                                                    tgtTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCaa                                                   >  2:852440/1‑139 (MQ=255)                                                                                         cTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCa                                              <  1:951237/139‑1 (MQ=255)                                                                                            ccGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTa                                           <  2:471788/139‑1 (MQ=255)                                                                                            ccGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTa                                           >  2:586693/1‑139 (MQ=255)                                                                                             cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGc                                                        >  2:1009583/1‑125 (MQ=255)                                                                                             cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGc                                                        <  1:1009583/125‑1 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAg                                                                                                                                      <  2:52058/46‑1 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAg                                                                                                                                      >  1:52058/1‑46 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTAc                                                                            <  2:1090009/104‑1 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTAc                                                                            >  1:1090009/1‑104 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAAt                                         >  2:595752/1‑139 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           <  1:672593/139‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           <  1:921715/139‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTATGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           >  1:863984/1‑139 (MQ=255)                                                                                                                      cGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGa                 >  2:141894/1‑139 (MQ=255)                                                                                                                         tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTc                                                     <  2:688146/100‑1 (MQ=255)                                                                                                                         tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTc                                                     >  1:688146/1‑100 (MQ=255)                                                                                                                         tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATc                                  >  2:235526/1‑119 (MQ=255)                                                                                                                         tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATc                                  <  1:235526/119‑1 (MQ=255)                                                                                                                         tacatacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGatat              >  2:846158/1‑139 (MQ=255)                                                                                                                             tacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCg          <  2:1177333/139‑1 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACt                                                                           >  1:1087272/1‑73 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACt                                                                           <  2:1087272/73‑1 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCgg         >  2:617637/1‑139 (MQ=255)                                                                                                                                 tacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGAc                               <  1:514720/114‑1 (MQ=255)                                                                                                                                 tacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGAc                               >  2:514720/1‑114 (MQ=255)                                                                                                                                   cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           >  1:1143464/1‑116 (MQ=255)                                                                                                                                   cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           <  2:1143464/116‑1 (MQ=255)                                                                                                                                   cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGatcat    >  1:94024/1‑139 (MQ=255)                                                                                                                                    tCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           <  1:194350/115‑1 (MQ=255)                                                                                                                                    tCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                           >  2:194350/1‑115 (MQ=255)                                                                                                                                     ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaaaa                                                                                     <  2:725409/56‑1 (MQ=255)                                                                                                                                     ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCg  >  2:215327/1‑139 (MQ=255)                                                                                                                                     ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCg  >  1:662386/1‑139 (MQ=255)                                                                                                                                     ccATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCg  >  1:1078092/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        TCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGGATCATCG  >  NC_000913/3439006‑3439276

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