Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,042,803	T→C	Y146H (TAC→CAC)	yodB →	cytochrome b561 homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,803	0	T	C	100.0%	136.2 / NA	38	Y146H (TAC→CAC)	yodB	cytochrome b561 homolog
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/22);  total (16/22)

CTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAGGAAATGTGAACCATGGAC  >  NC_000913/2042668‑2042933                                                                                                                                        |                                                                                                                                   ctgGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACt                                                                                                                                 >  1:19503/1‑139 (MQ=255)   ggggATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACttt                                                                                                                               >  1:362171/1‑139 (MQ=255)                   gCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACg                                                                                                               <  2:369389/139‑1 (MQ=255)                             gAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCgcagc                                                                                                     <  1:375015/139‑1 (MQ=255)                             gAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCgcagc                                                                                                     <  1:366368/139‑1 (MQ=255)                              aaaaTCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCgcagca                                                                                                    >  2:237028/1‑139 (MQ=255)                                aaTCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCAct                                                                                                  >  1:538970/1‑139 (MQ=255)                                  tCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCActct                                                                                                <  2:526294/139‑1 (MQ=255)                                   cGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTTACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCtt                                                                                               >  2:356489/1‑139 (MQ=255)                                        aGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTtcatc                                                                                          <  2:437867/139‑1 (MQ=255)                                              cTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACa                                                                                    <  1:491267/139‑1 (MQ=255)                                                tGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACAtt                                                                                  <  1:441081/139‑1 (MQ=255)                                                   tttCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAg                                                                               <  2:538970/139‑1 (MQ=255)                                                    ttCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGa                                                                              <  2:788988/139‑1 (MQ=255)                                                     tCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGaa                                                                             <  2:786817/139‑1 (MQ=255)                                                     tCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGaa                                                                             <  1:237028/139‑1 (MQ=255)                                                             tctcCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGataata                                                                     >  1:294106/1‑139 (MQ=255)                                                               tcCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATAct                                                                   >  1:400326/1‑139 (MQ=255)                                                                 ccTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATActct                                                                 >  2:195813/1‑139 (MQ=255)                                                                 ccTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATActct                                                                 <  2:294106/139‑1 (MQ=255)                                                                     tGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTa                                                             <  2:454890/139‑1 (MQ=255)                                                                         aCCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCgca                                                                                                       <  1:253841/93‑1 (MQ=255)                                                                         aCCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCgca                                                                                                       >  2:253841/1‑93 (MQ=255)                                                                                aCAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGcc                                                  <  2:362171/139‑1 (MQ=255)                                                                                 cAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCt                                                 <  1:356489/139‑1 (MQ=255)                                                                                    cGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTcgc                                              >  1:319946/1‑139 (MQ=255)                                                                                    cGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTcgc                                              >  2:490514/1‑139 (MQ=255)                                                                                    cGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTcgc                                              >  2:117912/1‑139 (MQ=255)                                                                                                 cTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCg                                 >  2:66354/1‑139 (MQ=255)                                                                                                  tGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTACATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGa                                >  1:573720/1‑139 (MQ=255)                                                                                                   gATAAAAAATATTCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGAt                               <  1:296093/139‑1 (MQ=255)                                                                                                              ttCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATc                              <  2:803301/129‑1 (MQ=255)                                                                                                              ttCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATc                              >  1:803301/1‑129 (MQ=255)                                                                                                              ttCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAg                    <  1:195813/139‑1 (MQ=255)                                                                                                               tCACGAAACCTGGGCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAgg                   <  2:400326/139‑1 (MQ=255)                                                                                                                          gggCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAGGAAATGTGAAcc        <  1:632384/139‑1 (MQ=255)                                                                                                                           ggCAAATATAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAGGAAATGTGAACCa       >  1:373522/1‑139 (MQ=255)                                                                                                                                atatAGGCCACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAGGAAATGTGAACCATGGAc  <  2:72142/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   CTGGGGATTGCTTTGATGGCTTACAGTGGAAAATCGTGGAGTTTCCTTGGTTTCAATGTGTCTCCCTTTGTTACCCCAAACAGCGAAATTAAAGCACTGATAAAAAATATTCACGAAACCTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACATTCAGAAAGATAATACTCTGTTACGAATGATGCCTCGCCGCAAATAAGTCGATCCCCTTACAAGGAAATGTGAACCATGGAC  >  NC_000913/2042668‑2042933

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