Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,521	A→G	I220V (ATA→GTA)	cycA →	D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,430,521	0	A	G	100.0%	94.3 / NA	28	I220V (ATA→GTA)	cycA	D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

GTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATG  >  NC_000913/4430380‑4430667                                                                                                                                              |                                                                                                                                                   gtgtCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGttt                                                                                                                                           >  2:122739/1‑151 (MQ=255)               cgtcgGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTgggg                                                                                                                             >  1:416414/1‑151 (MQ=255)                 tcgGCCTAGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGAt                                                                                                                           <  1:436769/151‑1 (MQ=255)                        ggTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAg                                                                                                                       >  1:94187/1‑148 (MQ=255)                           cATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta                                                                                                                 <  2:122338/151‑1 (MQ=255)                                tggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTac                                                                                                            <  1:286457/151‑1 (MQ=255)                                  gCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTacaa                                                                                                          >  2:15596/1‑151 (MQ=255)                                   cGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTacaac                                                                                                         <  1:16859/151‑1 (MQ=255)                                            ttCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGa                                                                                                >  2:200584/1‑151 (MQ=255)                                                           gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                                 >  2:178757/1‑151 (MQ=255)                                                           gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                                 >  1:405897/1‑151 (MQ=255)                                                           gagtgGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                                 >  1:30360/3‑151 (MQ=255)                                                                 aaGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATc                                                                           <  2:62876/151‑1 (MQ=255)                                                                      tCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                      <  2:26155/151‑1 (MQ=255)                                                                       cATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                      >  2:406105/1‑150 (MQ=255)                                                                                            aCGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTc                                                 <  2:26804/150‑1 (MQ=255)                                                                                                          ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                                  <  1:194773/151‑1 (MQ=255)                                                                                                             aaaGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGt                               <  1:282024/151‑1 (MQ=255)                                                                                                                ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                              >  2:213976/1‑149 (MQ=255)                                                                                                                 gTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                             <  2:284198/149‑1 (MQ=255)                                                                                                                 gTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                             <  2:37568/149‑1 (MQ=255)                                                                                                                    tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                        <  1:26421/151‑1 (MQ=255)                                                                                                                    tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                        <  1:163980/151‑1 (MQ=255)                                                                                                                     aaGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTc                       <  2:244961/151‑1 (MQ=255)                                                                                                                                       aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAtt     >  1:55208/1‑151 (MQ=255)                                                                                                                                       aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt      <  1:28307/150‑1 (MQ=255)                                                                                                                                       aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt      <  2:152142/150‑1 (MQ=255)                                                                                                                                          ccAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATg  >  1:415707/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                   GTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGATTATG  >  NC_000913/4430380‑4430667

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