Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,245,600	G→A	*147* (TAG→TAA)	ycgY →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,600	0	G	A	100.0%	149.4 / NA	46	*147* (TAG→TAA)	ycgY	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/23);  total (23/23)

CCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTA  >  NC_000913/1245468‑1245730                                                                                                                                     |                                                                                                                                   ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaaaa                                                                                                                              >  2:1361311/1‑139 (MQ=255)    aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           <  1:1361311/139‑1 (MQ=255)     cGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          >  2:458040/1‑139 (MQ=255)        cGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGc                                                                                                                       >  1:186657/1‑139 (MQ=255)         ggCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      <  1:783826/139‑1 (MQ=255)         ggCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      <  1:534937/139‑1 (MQ=255)           cGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTAt                                                                                                                    >  1:1409227/1‑139 (MQ=255)            gTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATg                                                                                                                   <  2:1409227/139‑1 (MQ=255)                       ggTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGc                                                                                                        >  2:388723/1‑139 (MQ=255)                                         cTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGgag                                                                                      <  1:458040/139‑1 (MQ=255)                                             cccGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      >  2:1300339/1‑103 (MQ=255)                                             cccGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      <  1:1300339/103‑1 (MQ=255)                                                  aGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTg                                                                                                         <  1:376202/111‑1 (MQ=255)                                                  aGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTg                                                                                                         >  2:376202/1‑111 (MQ=255)                                                   gAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          <  1:900659/93‑1 (MQ=255)                                                   gAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          >  2:900659/1‑93 (MQ=255)                                                                    ggCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                           <  2:186657/139‑1 (MQ=255)                                                                         ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTAt                                                                                                            >  2:1433866/1‑85 (MQ=255)                                                                         ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTAt                                                                                                            <  1:1433866/85‑1 (MQ=255)                                                                         ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTc                                                                    >  2:1494864/1‑125 (MQ=255)                                                                         ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTc                                                                    <  1:1494864/125‑1 (MQ=255)                                                                         ccACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTgg                                                      >  1:31792/1‑139 (MQ=255)                                                                                    aaaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACa                                                                                                   >  1:432776/1‑83 (MQ=255)                                                                                    aaaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACa                                                                                                   <  2:432776/83‑1 (MQ=255)                                                                                    aaaCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTgg                                           >  1:1275168/1‑139 (MQ=255)                                                                                      aCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          >  2:716761/1‑58 (MQ=255)                                                                                      aCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCAt                                                                                                                          <  1:716761/58‑1 (MQ=255)                                                                                      aCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGgag                                                                                      <  2:595845/94‑1 (MQ=255)                                                                                      aCATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGgag                                                                                      >  1:595845/1‑94 (MQ=255)                                                                                         tttGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                      <  2:1275168/139‑1 (MQ=255)                                                                                         tttGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                      <  1:388723/139‑1 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAAc                                                                                                                             >  1:445314/1‑50 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAAc                                                                                                                             <  2:445314/50‑1 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      <  2:1493888/57‑1 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                      >  1:1493888/1‑57 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                           >  1:472613/1‑116 (MQ=255)                                                                                           tGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                           <  2:472613/116‑1 (MQ=255)                                                                                                          gCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAggg                     >  2:610281/1‑139 (MQ=255)                                                                                                            cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           >  1:664260/1‑35 (MQ=38)                                                                                                            cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                           <  2:664260/35‑1 (MQ=38)                                                                                                               gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACg                                                                         >  2:280567/1‑82 (MQ=255)                                                                                                               gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACg                                                                         <  1:280567/82‑1 (MQ=255)                                                                                                               gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCtt                               <  2:124364/124‑1 (MQ=255)                                                                                                               gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCtt                               >  1:124364/1‑124 (MQ=255)                                                                                                                 cTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGcgtc              >  2:326886/1‑139 (MQ=255)                                                                                                                       aaaCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCt        <  1:149016/139‑1 (MQ=255)                                                                                                                             tCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTa  >  1:851180/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   CCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTA  >  NC_000913/1245468‑1245730

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