Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	657,745	C→T	W65* (TGG→TAG)	crcB ←	F(‑) channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	657,745	0	C	T	100.0%	107.0 / NA	31	W65* (TGG→TAG)	crcB	F(‑) channel
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/17);  total (14/17)

AGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACT  >  NC_000913/657609‑657872                                                                                                                                         |                                                                                                                                aGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCtaca                                                                                                                              >  1:68847/1‑140 (MQ=255)  gACCCCAGAAGGTTGACGAAAACGTTAAGTAATGCCCAGCCAAAGCGCCCCTCTTGTAACAAAAACCCCACTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTacac                                                                                                                             <  1:299428/140‑1 (MQ=255)   aCCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACt                                                                                                                            >  1:167599/1‑140 (MQ=255)               tGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTc                                                                                                                >  1:66734/1‑140 (MQ=255)                   aaaaaCGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGAAAAAAAAACACCACTTCTGCCGAAAATGTTGAGATGTTTTTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCa                                                                                                            <  2:167599/139‑1 (MQ=255)                            cAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGGGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGa                                                                                                   >  1:186091/1‑140 (MQ=255)                                   gCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGc                                                                                            <  1:298148/140‑1 (MQ=255)                                   gCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGc                                                                                            >  1:404836/1‑140 (MQ=255)                                    cccAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCg                                                                                           >  1:209325/1‑140 (MQ=255)                                    cccAGCAAAGGCGCCCCTCTTGTAAAAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTTTAAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCg                                                                                           <  2:410254/140‑1 (MQ=255)                                       aGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                                                                        >  1:157348/1‑140 (MQ=255)                                       aGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                                                                        >  2:125416/1‑140 (MQ=255)                                           aaaGCGGCCCTCTTGTAACAAAAACCCCACTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtccta                                                                                    <  1:404294/140‑1 (MQ=255)                                                          tAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATg                                                                     <  2:66734/140‑1 (MQ=255)                                                           aaaaaaaaCGCCACTTCTGCCGAAATTTTTAAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGc                                                                    <  2:568380/137‑1 (MQ=255)                                                            aCAAAAACCCCACTTCTGCCAAAAATGTTGAGAAGTTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcc                                                                   <  2:404836/140‑1 (MQ=255)                                                                 aaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAAt                                                              >  2:533889/1‑140 (MQ=255)                                                                 aaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAAt                                                              >  2:300911/1‑140 (MQ=255)                                                                            tGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTg                                                   <  1:533889/140‑1 (MQ=255)                                                                               cGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCa                                                <  1:300911/140‑1 (MQ=255)                                                                                gAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAg                                               >  2:26505/1‑140 (MQ=255)                                                                                              ggTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGaa                                 <  1:125416/140‑1 (MQ=255)                                                                                                ttgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATc                               >  1:593694/1‑140 (MQ=255)                                                                                                   ttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGcc                            <  2:157348/140‑1 (MQ=255)                                                                                                         cGCCACAAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGc                      >  2:709020/1‑140 (MQ=255)                                                                                                            cacaAAATCCGGTGGTGATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGt                   >  1:482720/1‑140 (MQ=255)                                                                                                                        tggtgATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGc       <  1:686822/140‑1 (MQ=255)                                                                                                                         ggtgATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCa      <  1:26505/140‑1 (MQ=255)                                                                                                                          gtgATTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCAt     <  2:482720/140‑1 (MQ=255)                                                                                                                             aTTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGc       <  2:351406/135‑1 (MQ=255)                                                                                                                             aTTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGc       >  1:351406/1‑135 (MQ=255)                                                                                                                             aTTAATACTTTCTACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACt  <  2:68847/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                AGACCCCAGAAGGTTGACGAAAACGTTCAGTAATGCCCAGCCAAAGCGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACT  >  NC_000913/657609‑657872

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