Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,401,103	G→A	T151T (ACC→ACT)	nuoF ←	NADH:ubiquinone oxidoreductase, chain F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,401,103	0	G	A	100.0%	108.5 / NA	32	T151T (ACC→ACT)	nuoF	NADH:ubiquinone oxidoreductase, chain F
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

TTTCTTCCCCGCAGATGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCGGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAACGCGGAGATGAGCATACC  >  NC_000913/2401006‑2401213                                                                                                  |                                                                                                               tttCTTCCCCGCAGATGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTc                                                                       >  2:431039/1‑139 (MQ=255)      tcCCCGCAGATGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGatat                                                                  >  1:441506/1‑139 (MQ=255)              gATGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCa                                                          <  2:441506/139‑1 (MQ=255)                tGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACg                                                        <  1:431039/139‑1 (MQ=255)                    gcgcCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAgg                                                    >  1:861916/1‑139 (MQ=255)                         cTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTACCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCACATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGAt                                               >  2:555904/1‑139 (MQ=255)                          tGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCCTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATg                                              <  1:315891/139‑1 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGgcg                                                                                         >  1:235525/1‑84 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGgcg                                                                                         <  2:235525/84‑1 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCt                                                                         >  1:312101/1‑100 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCt                                                                         <  2:312101/100‑1 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAgg                                                    <  1:332433/121‑1 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAgg                                                    >  2:332433/1‑121 (MQ=255)                                      gACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTaa                                  >  2:233287/1‑139 (MQ=255)                                              gTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAg                          >  2:786130/1‑139 (MQ=255)                                                 cgaaatcgaaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAATTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAgcgc                       >  1:538987/1‑139 (MQ=255)                                                 cgaaatcgaaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAgcgc                       >  2:504521/1‑139 (MQ=255)                                                      tcgaaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAAcg                  <  1:504521/139‑1 (MQ=255)                                                      tcgaaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAAcg                  <  2:861916/139‑1 (MQ=255)                                                         aaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGa                                                                                     >  1:422022/1‑69 (MQ=255)                                                         aaaCCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGa                                                                                     <  2:422022/69‑1 (MQ=255)                                                                      tAATGTTTTTGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAACGCGGAGATGAGCATAcc  <  1:786130/139‑1 (MQ=255)                                                                             tttGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAAcg                  <  2:464700/116‑1 (MQ=255)                                                                             tttGCCAAGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAAcg                  >  1:464700/1‑116 (MQ=255)                                                                                   aaGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACgg                                                                                           >  1:181020/1‑37 (MQ=255)                                                                                   aaGCAGACCCGCTTCAGTGGCTTCGGCAATGGCACgg                                                                                           <  2:181020/37‑1 (MQ=255)                                                                                     gCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCg                                                                           >  2:297599/1‑51 (MQ=255)                                                                                     gCAGACCCGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCg                                                                           <  1:297599/51‑1 (MQ=255)                                                                                           ccGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAACGCGga              <  2:125198/106‑1 (MQ=255)                                                                                           ccGCTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAACGCGga              >  1:125198/1‑106 (MQ=255)                                                                                              cTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACg                                                        >  2:501667/1‑61 (MQ=255)                                                                                              cTTCAGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACg                                                        <  1:501667/61‑1 (MQ=255)                                                                                                  |                                                                                                               TTTCTTCCCCGCAGATGTAGCGCCCTGCCCCGGTATGGACGAACAGTTCGAAATCGAAACCTGTTCCCATAATGTTTTTGCCAAGCAGACCCGCTTCGGTGGCTTCGGCAATGGCACGGCGCAGATTAACTGCCGCTTCGATATATTCGCCACGCAGGAAGATGTAGCCACGGTAAGCTTTCAGCGCAAACGCGGAGATGAGCATACC  >  NC_000913/2401006‑2401213

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