Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,957,269	A→G	G288G (GGT→GGC)	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,269	0	A	G	100.0%	128.2 / NA	36	G288G (GGT→GGC)	torY	cytochrome c quinol dehydrogenase TorY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/18);  total (18/18)

AATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTG  >  NC_000913/1957134‑1957404                                                                                                                                        |                                                                                                                                        aaTGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCt                                                                                                                                     <  2:128223/140‑1 (MQ=255)  aTGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTg                                                                                                                                    >  1:387127/1‑140 (MQ=255)              gCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTg                                                                                                                        <  1:134865/140‑1 (MQ=255)               cATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGa                                                                                                                       <  2:76821/140‑1 (MQ=255)                        tgtAATGGTTGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCattga                                                                                                              >  1:411691/1‑140 (MQ=255)                        tgtAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCattga                                                                                                              >  2:537700/1‑140 (MQ=255)                              ggTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCg                                                                                                        <  1:537700/140‑1 (MQ=255)                                  ggcggcAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTAtt                                                                                                    >  1:940910/1‑140 (MQ=255)                                  ggcggcAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTAtt                                                                                                    >  1:773652/1‑140 (MQ=255)                                  ggcggcAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTAtt                                                                                                    >  1:575152/1‑140 (MQ=255)                                    cggcAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCg                                                                                                  >  2:365210/1‑140 (MQ=255)                                         aTAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGc                                                                                             <  2:907055/140‑1 (MQ=255)                                         aTAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGc                                                                                             <  1:365210/140‑1 (MQ=255)                                                cATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCAt                                                                                      <  2:411691/140‑1 (MQ=255)                                                cATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCAt                                                                                      >  2:650171/1‑140 (MQ=255)                                                 aTGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCAtt                                                                                     >  1:940323/1‑140 (MQ=255)                                                  tGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTc                                                                                    <  1:292376/140‑1 (MQ=255)                                                  tGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTc                                                                                    <  1:816707/140‑1 (MQ=255)                                                            cgcAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACa                                                                          >  2:473260/1‑140 (MQ=255)                                                                aTAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGa                                                                      >  2:870451/1‑140 (MQ=255)                                                                 tAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAg                                                                     <  1:651340/140‑1 (MQ=255)                                                                         gTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGa                                                             >  1:116397/1‑140 (MQ=255)                                                                            aaGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATAc                                                          <  1:473260/140‑1 (MQ=255)                                                                            aaGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATAc                                                          <  1:217988/140‑1 (MQ=255)                                                                               tttttCAGCGTATTGCCATAGCGGGTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACctt                                                       >  2:248298/1‑140 (MQ=255)                                                                                     aGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCAc                                                 <  1:870451/140‑1 (MQ=255)                                                                                        gTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACtt                                               <  1:155385/139‑1 (MQ=255)                                                                                        gTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACtt                                               >  2:155385/1‑139 (MQ=255)                                                                                                 tAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTa                                     <  2:387127/140‑1 (MQ=255)                                                                                                        ttACGGTCCGCTAATGCGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAg                              >  2:312236/1‑140 (MQ=255)                                                                                                                        cGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGCCGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGtt                         >  2:482849/1‑129 (MQ=255)                                                                                                                        cGGGATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGtt                         <  1:482849/129‑1 (MQ=255)                                                                                                                          ggATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAg            >  1:490042/1‑140 (MQ=255)                                                                                                                          ggATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCATCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAg            >  2:789998/1‑140 (MQ=255)                                                                                                                           gATCGGTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAgg           <  1:446978/140‑1 (MQ=255)                                                                                                                                gTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTagag      >  2:1021827/1‑140 (MQ=255)                                                                                                                                gTAAGCGCGCCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTagag      <  2:80524/140‑1 (MQ=255)                                                                                                                                    gcgcgcCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCgtg  >  1:498061/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        AATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTATTGCCATAGCGGTTTACGGTCCGCTAATGCGGGATCGGTAAGCGCACCCTGTAAAGACGCCTGACGCCATTGATTGCCGTATTCGTCAGCCTGCCATTCCCCGTTTACAGTGAGCGCCTTGATACCTTCTTCACTTAATGTGGCTAGCGCCAGCCGTTGACCTGCCGCCAGGTAGAGCGTG  >  NC_000913/1957134‑1957404

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