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%	163.1 / NA	51	*150* (TAG→TAA)	ygaU	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/26);  total (25/26)

GACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTC  >  NC_000913/2796225‑2796466                                                                                                                 |                                                                                                                                  gACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTg                                                                                                         >  2:104451/1‑139 (MQ=255) gACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACGCTTCCGGAATACGCAACACTTg                                                                                                         >  1:156287/1‑139 (MQ=255)            cggcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATaaa                                                                                              >  1:1471789/1‑139 (MQ=255)                    cTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATcc                                                                                      <  2:1299063/139‑1 (MQ=255)                         gcgAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                        >  2:138442/1‑132 (MQ=255)                         gcgAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAt                                                                                        <  1:138442/132‑1 (MQ=255)                            agaTACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCtttt                                                                              <  2:487020/139‑1 (MQ=255)                                         gcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTa                                                                 >  1:66909/1‑139 (MQ=255)                                                   tAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCg                                                             <  1:23622/133‑1 (MQ=255)                                                   tAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCg                                                             >  2:23622/1‑133 (MQ=255)                                                   tAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                       >  1:142243/1‑139 (MQ=255)                                                     aTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                   >  2:620513/1‑109 (MQ=255)                                                     aTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                   <  1:620513/109‑1 (MQ=255)                                                       gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAt                                                   >  1:1241332/1‑139 (MQ=255)                                                       gcgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAt                                                   >  1:355408/1‑139 (MQ=255)                                                               aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGt                                           <  2:842285/139‑1 (MQ=255)                                                                   aCTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAg                                       <  2:1241332/139‑1 (MQ=255)                                                                       aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTa                                   >  1:583806/1‑139 (MQ=255)                                                                         gCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAattac                                                                                                                                >  1:370104/1‑44 (MQ=255)                                                                         gCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAattac                                                                                                                                <  2:370104/44‑1 (MQ=255)                                                                         gCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                          >  2:410308/1‑130 (MQ=255)                                                                         gCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                          <  1:410308/130‑1 (MQ=255)                                                                          ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCg                                <  2:66909/139‑1 (MQ=255)                                                                                cccATTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCg                          <  1:104451/139‑1 (MQ=255)                                                                                   aTTATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                       >  2:196864/1‑139 (MQ=255)                                                                                     tATCCAGCTGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGTCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                     <  1:894753/139‑1 (MQ=255)                                                                                     tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACg                                                                                                                  <  2:358000/46‑1 (MQ=255)                                                                                     tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACg                                                                                                                  >  1:358000/1‑46 (MQ=255)                                                                                     tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                     <  2:1471789/139‑1 (MQ=255)                                                                                     tATCCAGCAGGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGAc                     <  2:583806/139‑1 (MQ=255)                                                                                           gcagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                     >  1:1143492/1‑117 (MQ=255)                                                                                           gcagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                     <  2:1143492/117‑1 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGAt                                                                                                 >  2:1263108/1‑56 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGAt                                                                                                 <  1:1263108/56‑1 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                     <  2:820364/116‑1 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                     >  1:820364/1‑116 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg              >  1:1287203/1‑139 (MQ=255)                                                                                            cagGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTAACCGGGCTTTTTAGCATCGGTTTATTAGCCTCGAAGATTTTATTGTAAAGATTAGCGTTACCGTAGACCTGTTTg              >  2:506877/1‑139 (MQ=255)                                                                                             agGCGATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg             <  1:196864/139‑1 (MQ=255)                                                                                                 gATGGCTTGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTa                                               <  1:686168/101‑1 (MQ=255)                                                                                                 gATGGCTGGCAATTAATTACTCTTCCGGCATACGCAACACTTTCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTa                                               >  2:686168/1‑101 (MQ=255)                                                                                                 gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaa           <  2:228504/137‑1 (MQ=255)                                                                                                 gATGGCTGGCAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaa           >  1:228504/1‑137 (MQ=255)                                                                                                      ctggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                          >  2:837962/1‑101 (MQ=255)                                                                                                      ctggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                          <  1:837962/101‑1 (MQ=255)                                                                                                       tggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACt   <  2:156287/139‑1 (MQ=255)                                                                                                        ggcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTc  <  2:142243/139‑1 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     >  1:57294/1‑55 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCg                                                                                     <  2:57294/55‑1 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTa                                               >  2:88320/1‑93 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTa                                               <  1:88320/93‑1 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg             <  1:1154806/127‑1 (MQ=255)                                                                                                         gcAATTAATTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTgg             >  2:1154806/1‑127 (MQ=255)                                                                                                                 |                                                                                                                                  GACCCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTC  >  NC_000913/2796225‑2796466

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: ‑