Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,302,268	C→T	Q363* (CAG→TAG)	oppA →	oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,302,268	0	C	T	100.0%	72.2 / NA	24	Q363* (CAG→TAG)	oppA	oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/9);  total (15/9)

TCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGT‑ACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACT  >  NC_000913/1302051‑1302483                                                                                                                                                                                                                           |                                                                                                                                                                                                                        tCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGT‑ACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAActgc                                                                                                                                                                                        >  1:169444/1‑251 (MQ=255)                            caaTTACAAAATTAAAAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGT‑ACCGCACGGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTAAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGNc                                                                                                                                                              <  1:348129/248‑1 (MQ=255)                                 aaaaaaTAAACAACCAAAAACCCCCATTAAACGAGGTGCTGGTG‑GTGCCCGCATGGAAACTGGGTAGGGGCCGCGATATCATGGTTAAAAAAGTGAAAGCGCGGGGCAACAGGCCCGCCTATGGTTACACTCCACCGTATACTGAGGGCGCAAAAGTGACTCAGCCAGAAGGGTTTGGCTGGAGCTAGGAAAAACTTATCGAAGAAGCGAAAAAACTGCAGGCTGAAGCGGGTTATACCGCAGACAAACCg                                                                                                                                                        <  1:298198/248‑1 (MQ=255)                                                                                                   aTGGACCGCGATATCATTTTTAATAAAGTGAAAGCGCAGGGCAACATCCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAAAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTc                                                                                       <  1:443689/251‑1 (MQ=255)                                                                                                           ctatatCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTgctggc                                                                                                                                                                                    >  1:406829/3‑150 (MQ=255)                                                                                                                tCATTTTTAATAAAGGGAAATCGCAGGGCAAAATGCCCGCCTATGGTTACACTCCACCGTATACGGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGt                                                                                                                                                                          <  1:319276/155‑1 (MQ=255)                                                                                                                tCATTGTTNANNNANNGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTNTACTGATGGCGCAAANTTNACTCNGCCNGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGt                                                                                                                                                                          >  2:319276/1‑155 (MQ=255)                                                                                                                             aGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATAATGATGGCGCAAAATTGCCTAAGCCCGAATGGCTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAACAAACTGCTGGCTGAAGCGGGTTATACCGCa                                                                                                                                                                 >  1:66651/1‑151 (MQ=255)                                                                                                                                       cAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGCCTCCGCCAGACTGGTTTGGCTGGCGCTAGGAAAAACGTAACGCAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAAAACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACt                                                   >  1:415306/1‑251 (MQ=255)                                                                                                                                            naaCATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGc                                                                                                                                                                                                   >  1:93532/2‑102 (MQ=255)                                                                                                                                            cAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATc                                                                                                                      >  1:316665/1‑179 (MQ=255)                                                                                                                                                  gCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTCGGACAAAAGTAACGAAGAAGCGACAAACCTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGgaaga                                                                            >  1:334777/1‑215 (MQ=255)                                                                                                                                                    cnGCNNATGGNTACACTCCACCGTATACTGATGGCNCAAAATTGACTCAGCCAGAATGGTTGGGCTGNANCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCTGGNTATACCGCAGNCAANNCGTTNANCNTCAACCTGTTGTATAACACCTCCGATCTGNATAAAAAGCTGGCGATTGCTGCCNNTTNNNTGNNgaaga                                                                            <  2:334777/213‑1 (MQ=255)                                                                                                                                                                   cTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTACCGAAGAAGCGAAAAAAATGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGCGTGGAAAACGTTCCTc                       >  1:177629/1‑251 (MQ=255)                                                                                                                                                                               cTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAg                                                                                                                                                                               >  1:7642/1‑87 (MQ=255)                                                                                                                                                                               cTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAg                                                                                                                                                                               <  2:7642/87‑1 (MQ=255)                                                                                                                                                                                   tGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAActgctg                                                                                                                                                                                      >  2:19590/1‑76 (MQ=255)                                                                                                                                                                                   tGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAActgctg                                                                                                                                                                                      <  1:19590/76‑1 (MQ=255)                                                                                                                                                                                        cAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCTAGGAAAAACGTAAAGCAGAAGCGAACAAACTGCTGGCTGACGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGAAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACt  >  1:295110/1‑251 (MQ=255)                                                                                                                                                                                                        aGAATGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGTAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGCACATCAACATGTTGTATAACACCTCAGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAAAATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGaaaa                               >  1:218816/1‑206 (MQ=255)                                                                                                                                                                                                           aTGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGAccn                                                                                                                                                 >  2:31454/1‑88 (MQ=255)                                                                                                                                                                                                           aTGGTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCa                                                                                                                                                 <  1:31454/89‑1 (MQ=255)                                                                                                                                                                                                             ggTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATNNACCtgtnnt                                                                                                                                     >  2:22743/1‑99 (MQ=255)                                                                                                                                                                                                             ggTTTGGCTGGAGCTAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACg                             <  1:22743/203‑1 (MQ=255)                                                                                                                                                                                                              gTTTGNCTGGANCNNGGAAAAACGTCACGAAGAANCGAAAAAACTgctggn                                                                                                                                                                                    >  2:190952/1‑50 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                        TCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGT‑ACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACT  >  NC_000913/1302051‑1302483

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