Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	345,649	G→A	*82* (TAG→TAA)	yahM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,649	0	G	A	100.0%	160.8 / NA	50	*82* (TAG→TAA)	yahM	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/24);  total (26/24)

AAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCAAA  >  NC_000913/345513‑345783                                                                                                                                         |                                                                                                                                       aaaaaGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTaaag                                                                                                                                      <  1:636686/139‑1 (MQ=255) aaaaaGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCCAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTaaag                                                                                                                                      <  1:265066/139‑1 (MQ=255)     aGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTa                                                                                                                                  >  1:673789/1‑139 (MQ=255)                  gagaGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCa                                                                                                                     >  1:493453/1‑139 (MQ=255)                    gagTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCAcc                                                                                                                   <  2:474229/139‑1 (MQ=255)                     agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                                  >  1:410379/1‑139 (MQ=255)                      gTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGc                                                                                                                 <  2:673789/139‑1 (MQ=255)                          gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            <  2:514764/124‑1 (MQ=255)                          gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                            >  1:514764/1‑124 (MQ=255)                             gAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCa                                                                                                          <  2:493453/139‑1 (MQ=255)                                aTGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACTACGCATTACGGGGGCTATCATCGCCGGTGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                                       <  2:920907/139‑1 (MQ=255)                                        cTAAAGACTCAAGTAAACTTATCGGAAATAAGTACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTaa                                                                                               >  2:327261/1‑139 (MQ=255)                                        cTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTaa                                                                                               >  1:834215/1‑139 (MQ=255)                                          aaaGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             <  2:410379/139‑1 (MQ=255)                                                 cAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGt                                                                                      <  1:935691/139‑1 (MQ=255)                                                  aaGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTa                                                                                     >  1:1026473/1‑139 (MQ=255)                                                        aCTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCg                                                                                                  >  1:993307/1‑120 (MQ=255)                                                        aCTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCg                                                                                                  <  2:993307/120‑1 (MQ=255)                                                           tATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATa                                                                            <  2:1232676/139‑1 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAgg                                                                                         >  1:674044/1‑125 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAgg                                                                                         <  2:674044/125‑1 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAg                                                                                    <  1:673358/130‑1 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAg                                                                                    >  2:673358/1‑130 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATac                                                                           >  1:756853/1‑139 (MQ=255)                                                            aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATac                                                                           >  2:756441/1‑139 (MQ=255)                                                              cGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                             >  1:148187/1‑135 (MQ=255)                                                              cGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                             <  2:148187/135‑1 (MQ=255)                                                                    tAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCaataa                                                                   <  2:1043519/139‑1 (MQ=255)                                                                       ggACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                                <  2:1026473/139‑1 (MQ=255)                                                                           cACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCAcc                                                            >  2:1097963/1‑139 (MQ=255)                                                                                        ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTg                                               >  1:1065899/1‑139 (MQ=255)                                                                                        ggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTg                                               >  2:1223778/1‑139 (MQ=255)                                                                                            atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                                       >  2:22613/1‑79 (MQ=255)                                                                                            atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                                       <  1:22613/79‑1 (MQ=255)                                                                                                tcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAAcgcg                                       <  1:217928/139‑1 (MQ=255)                                                                                                    cGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAg                                                                                    >  2:161202/1‑90 (MQ=255)                                                                                                    cGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAg                                                                                    <  1:161202/90‑1 (MQ=255)                                                                                                     gATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCa                                  <  2:527496/139‑1 (MQ=255)                                                                                                      aTGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCAt                                 <  1:423746/139‑1 (MQ=255)                                                                                                       tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                                >  1:462467/1‑139 (MQ=255)                                                                                                        ggCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             >  2:1192570/1‑77 (MQ=255)                                                                                                        ggCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                             <  1:1192570/77‑1 (MQ=255)                                                                                                             tGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCa                          >  2:125096/1‑139 (MQ=255)                                                                                                               aTTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTc                                                                                                                      >  2:1073192/1‑45 (MQ=255)                                                                                                               aTTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTc                                                                                                                      <  1:1073192/45‑1 (MQ=255)                                                                                                                        aTCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCacgacg               <  2:756853/139‑1 (MQ=255)                                                                                                                          cAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGca             >  1:608838/1‑139 (MQ=255)                                                                                                                             gTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAg          >  1:821300/1‑139 (MQ=255)                                                                                                                                   tACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                                <  1:758067/79‑1 (MQ=255)                                                                                                                                   tACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCa    <  1:1223778/139‑1 (MQ=255)                                                                                                                                   tACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGAAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                                >  2:758067/1‑79 (MQ=255)                                                                                                                                     cGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCaaa  >  1:595902/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       AAAAAGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGGCAAA  >  NC_000913/345513‑345783

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