Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,957,999	A→G	M45T (ATG→ACG)	torY ←	cytochrome c quinol dehydrogenase TorY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,957,999	0	A	G	100.0%	183.2 / NA	52	M45T (ATG→ACG)	torY	cytochrome c quinol dehydrogenase TorY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/18);  total (34/18)

TATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTCCCT  >  NC_000913/1957862‑1958127                                                                                                                                          |                                                                                                                                 tatCTTTAGATGCCTTTAATTTTGCAAATAATTAATCCAACCCTGACTATGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTc                                                                                                                                <  2:280273/140‑1 (MQ=255)   tCTTTAGAGGCCTTTATTTTAGCAAATAAATCAGCCATCTCTGACTTTAGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGa                                                                                                                              <  2:309374/140‑1 (MQ=255)     tttAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACGGTTCCCTGATATTCTTCAAAAGGTTTACTCGTCGAGt                                                                                                                            >  1:218507/1‑140 (MQ=255)            gCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGa                                                                                                                     >  1:497733/1‑140 (MQ=255)               tttAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAg                                                                                                                  >  2:519563/1‑140 (MQ=255)               tttAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAg                                                                                                                  >  2:111055/1‑140 (MQ=255)                 tAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGc                                                                                                                >  1:651808/1‑140 (MQ=255)                   aTTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCaa                                                                                                              <  1:406840/140‑1 (MQ=255)                       aGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCGAAAGGTTTACTGGTGGAGTGGAAGGAAAGGCAAAATg                                                                                                          >  1:203974/1‑140 (MQ=255)                       aGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATg                                                                                                          >  1:162654/1‑140 (MQ=255)                         cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTTGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCt                                                                                                        >  2:550501/1‑140 (MQ=255)                         cAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGACAAAATGCt                                                                                                        <  2:624917/140‑1 (MQ=255)                             taaataaTCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTAt                                                                                                    <  1:431028/140‑1 (MQ=255)                                    tccatccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGt                                                                                             >  2:390024/1‑140 (MQ=255)                                       atccCTGACTTTGGAATCTGGCAATCGGCACCTTCGGCAGGTATCCCTATCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGtttt                                                                                          <  2:185285/140‑1 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAATGCCAAATGCCGTATCCGACGTTTTa                                                                                         >  1:242017/1‑140 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGttttt                                                                                         >  2:102615/1‑139 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGtt                                                                                            <  2:133607/137‑1 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGtt                                                                                            >  1:133607/1‑137 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCCGTATCCGACGTTTTa                                                                                         >  1:104855/1‑140 (MQ=255)                                        tccCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGAAAAAGCAAAATGCTGTATCCGACGTTTTa                                                                                         >  1:210428/1‑140 (MQ=255)                                         cccTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTAt                                                                                        >  1:626659/1‑140 (MQ=255)                                         cccTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTAt                                                                                        >  2:112548/1‑140 (MQ=255)                                         cccTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTGTACTCGTCTAGTGGCATGAAAGGCAAAATGCTGTATCCGACGTTTTAt                                                                                        >  2:216627/1‑140 (MQ=255)                                          ccTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCAAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATg                                                                                       >  2:67661/1‑140 (MQ=255)                                                        atGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTgcg                                                                         <  1:112548/140‑1 (MQ=255)                                                          ggCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAAGAAAGGCAGAATACTGTATCCGACGATTTAAGTAAGACTTTTTGCGcc                                                                       >  2:250260/1‑140 (MQ=255)                                                            cAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAg                                                                     >  1:603972/1‑140 (MQ=255)                                                                            cACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGcca                                                     >  1:97966/1‑140 (MQ=255)                                                                                 aTCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCaccaac                                                <  1:43066/140‑1 (MQ=255)                                                                                 aTCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCaccaac                                                >  2:43066/1‑140 (MQ=255)                                                                                        tCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAAcc                                         >  1:622900/1‑140 (MQ=255)                                                                                        tCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAAcc                                         >  2:6724/1‑140 (MQ=255)                                                                                                    aGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGGGTGGCAGGAAAGGCAAAATGCTGTAACCGAAGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCCCCAACCACAACCGCTATCAGCaaa                             <  2:647819/140‑1 (MQ=255)                                                                                                    aGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCaaa                             <  1:67661/140‑1 (MQ=255)                                                                                                    aGTGGACAGTTCCCTGAGATTCCTCAAAAGGATTAATCGACGAGTGGCCGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCaaa                             <  2:210428/140‑1 (MQ=255)                                                                                                       ggACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCaaa                             <  2:164198/137‑1 (MQ=255)                                                                                                       ggACAGTTCCCTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCaaa                             >  1:164198/1‑137 (MQ=255)                                                                                                                ccTGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAAAAACCCGCCGCCACCAAACACAACCGATATCCGCAAAAACAATAACCCa                 >  2:215460/1‑140 (MQ=255)                                                                                                                  tGATATTCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAGGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACTACAACCGCTATCAGCAAAAACAATAACCCAAt               >  2:412972/1‑140 (MQ=255)                                                                                                                       ttCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAGGACTTTTTGCGCCAGCAATAACCCGGCGCCACCAGCCAAAACCGCTATCAGAAAAAGCAATAACCCAATGCGtt          >  2:172899/1‑140 (MQ=255)                                                                                                                       ttCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGtt          >  1:361556/1‑140 (MQ=255)                                                                                                                       ttCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGtt          >  2:527461/1‑140 (MQ=255)                                                                                                                       ttCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAAAGCTGTATCCGACGTTTTATGTAAGACTTTTTGGGCCAGCAATAACCCGCCGCCACCACCCAAAACAGCTATCAGGAAAAACAAAAATCCAATGCGtt          >  2:270773/1‑140 (MQ=255)                                                                                                                        tCCTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGttt         >  1:381426/1‑140 (MQ=255)                                                                                                                         ccTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCTCCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGtttt        <  2:97966/140‑1 (MQ=255)                                                                                                                          cTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACGACCGCTACCAGGAAGAACAATGACCCAATGCGttttt       >  1:177264/1‑140 (MQ=255)                                                                                                                          cTCAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGttttt       >  1:627606/1‑140 (MQ=255)                                                                                                                            cAAAAGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTc     <  2:651808/140‑1 (MQ=255)                                                                                                                             aaaaGGTTTACTCGTCGTGTGGGAGGAAAGGGAAACTGCGGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTcc    <  2:218507/140‑1 (MQ=255)                                                                                                                             aaaaGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTcc    <  1:111055/140‑1 (MQ=255)                                                                                                                             aaaaGGTTTACTCGTCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAAAAACCCAATGCGTTTTTTcc    <  2:644075/140‑1 (MQ=255)                                                                                                                             aaaaGGTTTAATCGTCGAATGGCAGGAAAGGCAAAATACTGTATACGGCGTTTTATGTAAGTTTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTcc    <  2:162654/140‑1 (MQ=255)                                                                                                                               aaGGTTTACTCGTCGCGTGGCACGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTCCCt  <  1:250260/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                 TATCTTTAGATGCCTTTAATTTAGCAAATAAATAATCCATCCCTGACTTTGGAATATGGCAATCGGCACATTCCGCACGTATCCCTTTCTGGTTCGAAAAGTGGACAGTTCCCTGATATTCCTCAAAAGGTTTACTCATCGAGTGGCAGGAAAGGCAAAATGCTGTATCCGACGTTTTATGTAAGACTTTTTGCGCCAGCAATAACCCGCCGCCACCAACCACAACCGCTATCAGCAAAAACAATAACCCAATGCGTTTTTTCCCT  >  NC_000913/1957862‑1958127

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