Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,401,795	G→A	R86R (CGC→CGT)	nuoE ←	NADH:ubiquinone oxidoreductase, chain E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,401,795	0	G	A	100.0%	170.6 / NA	45	R86R (CGC→CGT)	nuoE	NADH:ubiquinone oxidoreductase, chain E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/27);  total (18/27)

TCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGGCGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAGCGCTT  >  NC_000913/2401647‑2401940                                                                                                                                                     |                                                                                                                                                  tCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACg                                                                                                                                                 >  2:248310/1‑151 (MQ=255)  ccccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACg                                                                                                                                                 <  1:603671/150‑1 (MQ=255)   cccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACg                                                                                                                                                 <  1:611281/149‑1 (MQ=255)    ccAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGAcc                                                                                                                                              >  1:407780/1‑151 (MQ=255)       ggCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAAc                                                                                                                                           <  2:407589/151‑1 (MQ=255)         cagcaAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCg                                                                                                                                         >  1:344103/1‑151 (MQ=255)             aaGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTg                                                                                                                                     >  2:579286/1‑151 (MQ=255)               gTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGc                                                                                                                                   <  2:365402/151‑1 (MQ=255)                  ggCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCgg                                                                                                                                 <  2:15273/150‑1 (MQ=255)                  ggCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGa                                                                                                                                <  1:7552/151‑1 (MQ=255)                             aaaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACt                                                                                                                      <  1:617097/150‑1 (MQ=255)                              aaGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGt                                                                                                                    <  1:626740/151‑1 (MQ=255)                                 cGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGt                                                                                                                    >  1:150720/1‑148 (MQ=255)                                   gCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAAc                                                                                                               <  1:526431/151‑1 (MQ=255)                                            tGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCcac                                                                                                       >  1:130484/1‑150 (MQ=255)                                                gTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCtt                                                                                                  >  2:463980/1‑151 (MQ=255)                                                  ttGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCg                                                                                                >  2:344229/1‑151 (MQ=255)                                                    gCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGa                                                                                               <  1:414418/150‑1 (MQ=255)                                                     cccTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACg                                                                                             >  2:385344/1‑151 (MQ=255)                                                          gTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCt                                                                                        <  1:233082/151‑1 (MQ=255)                                                           tttGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCtt                                                                                       >  2:457034/1‑151 (MQ=255)                                                                 gTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCgg                                                                                  <  1:579229/150‑1 (MQ=255)                                                                     aGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATAc                                                                             <  1:694779/151‑1 (MQ=255)                                                                            tctcGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGcac                                                                      >  1:317577/1‑151 (MQ=255)                                                                             ctcGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGcac                                                                      >  1:156022/1‑150 (MQ=255)                                                                               cGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCt                                                                                                   >  1:629668/1‑119 (MQ=255)                                                                               cGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCt                                                                                                   <  2:629725/119‑1 (MQ=255)                                                                               cGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACAt                                                                    <  2:176397/150‑1 (MQ=255)                                                                                gAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCg                                                                  <  1:38175/151‑1 (MQ=255)                                                                                 aGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                                 >  2:524733/1‑151 (MQ=255)                                                                                 aGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                                 <  1:641167/151‑1 (MQ=255)                                                                                 aGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                                 <  1:530338/151‑1 (MQ=255)                                                                                 aGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCgg                                                                 >  1:291852/1‑151 (MQ=255)                                                                                       gccTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATcg                                                           <  1:101458/151‑1 (MQ=255)                                                                                        ccTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATcgc                                                          <  2:674970/151‑1 (MQ=255)                                                                                               aCCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGAt                                                    <  1:173920/150‑1 (MQ=255)                                                                                               aCCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATc                                                   <  2:134746/151‑1 (MQ=255)                                                                                                                 atGACAGACCACGCTGTCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCa                                  <  2:317610/150‑1 (MQ=255)                                                                                                                                 tCACAATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGa                  >  1:704683/1‑150 (MQ=255)                                                                                                                                  cacaATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAAc                <  2:101070/151‑1 (MQ=255)                                                                                                                                  cacaATAACGGATCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAAc                <  1:641767/151‑1 (MQ=255)                                                                                                                                             aTCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAgcg     >  2:406746/1‑151 (MQ=255)                                                                                                                                             aTCACATGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAgcg     <  1:535146/151‑1 (MQ=255)                                                                                                                                                acaTGACGACCAACCGGCTGGCGGAAGATCTGACTGTCGAACGTTGCCACACTTTCTACGTCGCCTGCCGGCATACCCCGCACATCGGCGATCGCGTGGATCGCACCATCCGGCCCCCAGCCACGCTGTTTCTGAACGATTTTCAgt      <  1:241840/147‑2 (MQ=255)                                                                                                                                                acaTGACGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATTGCGTGGATCGCACCATCCGGCACCCCGCCACGCTGCTTCTGAACGATTTTCAGCGCtt  >  1:769017/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCCCCAGGCAGCAAGTTGGCAGCAGCGTAAAGCGGCCATCAAATGTCGTTTGCCCTGGTTTGATGTTCAGCTTTTTCTCGAGCGCCGCCTGAATACCCTGATAACCGTTGATATGACAGACCACGCTGTCACAATAACGGATCACATGGCGACCAACCGGCTGGCGGAAGATCTGACTGTAGAACGTTGCCACACCTTCGACGTCGCTTGCCGGAATACCCAGCACATCGGCGATCGCGTGGATCGCACCATCCGGCACCCAGCCACGCTGCTTCTGAACGATTTTCAGCGCTT  >  NC_000913/2401647‑2401940

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