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%	179.4 / NA	49	I220V (ATA→GTA)	cycA	D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/25);  total (24/25)

CCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAT  >  NC_000913/4430375‑4430663                                                                                                                                                   |                                                                                                                                               ccATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAg                                                                                                                                             >  1:125485/1‑150 (MQ=255)        gtCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGtttt                                                                                                                                      >  1:372632/1‑150 (MQ=255)                         gCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGa                                                                                                                    >  1:347210/1‑151 (MQ=255)                           cTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGc                                                                                                                  <  1:68283/151‑1 (MQ=255)                           cTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGc                                                                                                                  <  1:44305/151‑1 (MQ=255)                                 aTGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggtag                                                                                                            <  1:460307/151‑1 (MQ=255)                                    gtggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggtaggt                                                                                                          <  1:234094/150‑1 (MQ=255)                                     tggCGATGCACTTTCAGTCACCGTCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTac                                                                                                        <  1:343221/151‑1 (MQ=255)                                     tggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTac                                                                                                        <  1:170026/151‑1 (MQ=255)                                     tggCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTac                                                                                                        <  2:332525/151‑1 (MQ=255)                                        cGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTacaac                                                                                                     <  1:341478/151‑1 (MQ=255)                                           tGCACTTTCATTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGTTTTCGTGGGGATTGAGCTGGTAGGTACAACAGc                                                                                                  <  1:586114/151‑1 (MQ=255)                                           tGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGc                                                                                                  <  1:171224/151‑1 (MQ=255)                                               cTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCg                                                                                              >  1:383367/1‑151 (MQ=255)                                                 ttCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGa                                                                                            >  1:368190/1‑151 (MQ=255)                                                     gTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAAc                                                                                         >  2:189211/1‑150 (MQ=255)                                                          cggcTGGTGGGGATGCGTCCTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCTCGAAAGGTTTAAGTGTCTTCTTTGCCGGAGTCCAGGTATCGGGTTTCGCTTTCGTGGGGATTGAGAAGGTAGCTACAGCAGGTGCGGAGACCaaaac                                                                                   <  1:526718/148‑3 (MQ=255)                                                          cGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGa                                                                                   <  2:418330/151‑1 (MQ=255)                                                               ggTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCag                                                                              >  2:566988/1‑151 (MQ=255)                                                               ggTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCag                                                                              >  2:617158/1‑151 (MQ=255)                                                               ggTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCag                                                                              >  1:83123/1‑151 (MQ=255)                                                               ggTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGCTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCag                                                                              >  1:77146/1‑151 (MQ=255)                                                                gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  2:292941/1‑151 (MQ=255)                                                                gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  2:443315/1‑151 (MQ=255)                                                                gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  2:502309/1‑151 (MQ=255)                                                                gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  2:687735/1‑151 (MQ=255)                                                                    ggAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGaaa                                                                         >  2:519264/1‑151 (MQ=255)                                                                                      tttGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGAtt                                                       >  2:99809/1‑151 (MQ=255)                                                                                         gtgGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTaa                                                     >  1:458672/1‑150 (MQ=255)                                                                                         gtgGAATGACGGCGGCTGGCTCCCGAAAGGTTTAAGTGGCTTCTTTGACGGATTCCAGGTAGCGGATTTAGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCGGCGGTAACCAAAGATCCAGAGAAATCACTCCCACGCGCGATTaag                                                    >  1:360670/1‑150 (MQ=255)                                                                                             aaTGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCa                                                >  2:372080/1‑151 (MQ=255)                                                                                                 aCGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTc                                             >  2:435624/1‑150 (MQ=255)                                                                                                     cggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATc                                        <  2:99057/151‑1 (MQ=255)                                                                                                     cggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATc                                        <  2:375511/151‑1 (MQ=255)                                                                                                           gTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTAt                                   <  2:233628/150‑1 (MQ=255)                                                                                                              cccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCAtt                               >  2:3408/1‑151 (MQ=255)                                                                                                               ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                              <  1:229232/151‑1 (MQ=255)                                                                                                               ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                              <  2:381204/151‑1 (MQ=255)                                                                                                                     ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  2:463854/1‑149 (MQ=255)                                                                                                                     ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  2:508734/1‑149 (MQ=255)                                                                                                                     ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  2:308420/1‑149 (MQ=255)                                                                                                                      gTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                         <  1:297203/149‑1 (MQ=255)                                                                                                                         tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                    <  2:125483/151‑1 (MQ=255)                                                                                                                         tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                    <  2:330568/151‑1 (MQ=255)                                                                                                                         tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                    <  1:308419/151‑1 (MQ=255)                                                                                                                                 ttcttTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCgcgc            <  1:426194/151‑1 (MQ=255)                                                                                                                                        ccGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATtg      <  2:299574/150‑1 (MQ=255)                                                                                                                                            aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  <  2:396394/150‑1 (MQ=255)                                                                                                                                            aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  <  2:114035/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAT  >  NC_000913/4430375‑4430663

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