Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,796,337	C→T	*150* (TAG→TAA)	ygaU ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,796,337	0	C	T	100.0%	160.6 / NA	50	*150* (TAG→TAA)	ygaU	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/24);  total (26/24)

GGAGGATAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCA  >  NC_000913/2796201‑2796467                                                                                                                                         |                                                                                                                                   ggaggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                  <  1:1171614/139‑1 (MQ=255)   aggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAattact                                                                                                                                >  2:235877/1‑139 (MQ=255)   aggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAattact                                                                                                                                >  1:152836/1‑139 (MQ=255)    ggaTAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTActc                                                                                                                               <  2:1534596/139‑1 (MQ=255)       tAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTc                                                                                                                            >  2:315023/1‑139 (MQ=255)         aaCACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCg                                                                                                                          >  2:838404/1‑139 (MQ=255)           cacGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                        <  2:781630/139‑1 (MQ=255)               aTCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATAc                                                                                                                    >  2:736216/1‑139 (MQ=255)                        aGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACtt                                                                                                           >  2:708400/1‑139 (MQ=255)                         gACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTg                                                                                                          >  2:38338/1‑139 (MQ=255)                          aCCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                         >  1:913447/1‑139 (MQ=255)                              tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                         <  1:1119283/135‑1 (MQ=255)                              tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGc                                                                                                         >  2:1119283/1‑135 (MQ=255)                                      gcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtt                                                                                             >  1:1349470/1‑139 (MQ=255)                                          tACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                         >  1:1348501/1‑139 (MQ=255)                                               gcgcGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                    <  2:1314031/139‑1 (MQ=255)                                                         ccGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                          <  2:803431/139‑1 (MQ=255)                                                          cgcgcgcgCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCAt                                                                         >  1:406549/1‑139 (MQ=255)                                                            cgcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGcc                                                                                                        <  2:893484/106‑1 (MQ=255)                                                            cgcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGcc                                                                                                        >  1:893484/1‑106 (MQ=255)                                                             gcgcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                      >  2:69642/1‑139 (MQ=255)                                                                           tAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                        <  1:838404/139‑1 (MQ=255)                                                                           tAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                        <  1:57872/139‑1 (MQ=255)                                                                             aTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAg                                                      >  2:212280/1‑139 (MQ=255)                                                                              tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGa                                                     <  1:736216/139‑1 (MQ=255)                                                                                cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAtt                                                   <  1:964381/139‑1 (MQ=255)                                                                                  ccAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAattaatta                                                                                                                                  >  1:356826/1‑58 (MQ=255)                                                                                  ccAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAtttt                                                 >  2:238794/1‑139 (MQ=255)                                                                                       aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGt                                            <  1:212280/139‑1 (MQ=255)                                                                                               aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTc                                                                                                                            <  2:242374/51‑1 (MQ=255)                                                                                               aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTc                                                                                                                            >  1:242374/1‑51 (MQ=255)                                                                                               aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTa                                    >  1:160359/1‑139 (MQ=255)                                                                                                cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                        >  1:406086/1‑54 (MQ=255)                                                                                                cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGa                                                                                                                        <  2:406086/54‑1 (MQ=255)                                                                                                cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGGATACGCAACACTTGCCCCGGGTAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAg                                   >  1:1631002/1‑139 (MQ=255)                                                                                                cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAg                                   >  2:1079022/1‑139 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGGATACGCAACACTTGCCCCGGGTAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  2:1631002/139‑1 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  2:152836/139‑1 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  2:160359/139‑1 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  1:315023/139‑1 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  1:235877/139‑1 (MQ=255)                                                                                                             tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                      <  1:238794/139‑1 (MQ=255)                                                                                                              aTCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAcc                     >  2:619197/1‑139 (MQ=255)                                                                                                                 cagcagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGt                  >  2:1378097/1‑139 (MQ=255)                                                                                                                    cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg               >  2:48089/1‑139 (MQ=255)                                                                                                                       gCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaa            <  2:1349470/139‑1 (MQ=255)                                                                                                                         gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt          <  2:1348501/139‑1 (MQ=255)                                                                                                                          aTGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                        <  2:1412149/92‑1 (MQ=255)                                                                                                                          aTGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                        >  1:1412149/1‑92 (MQ=255)                                                                                                                                 gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGt                                <  1:253562/109‑1 (MQ=255)                                                                                                                                 gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGt                                >  2:253562/1‑109 (MQ=255)                                                                                                                                 gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt          >  2:67147/1‑131 (MQ=255)                                                                                                                                 gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAAt          <  1:67147/131‑1 (MQ=255)                                                                                                                                 gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCa  >  2:662556/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   GGAGGATAAACACGATCAGCACCAGACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCA  >  NC_000913/2796201‑2796467

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