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%	181.8 / 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 (22/27);  total (22/27)

GTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAT  >  NC_000913/4430392‑4430663                                                                                                                                  |                                                                                                                                               gtcgtcGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTgg                                                                                                                           >  1:458731/1‑151 (MQ=255)     tcgGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGAt                                                                                                                       <  1:494211/151‑1 (MQ=255)        gCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTg                                                                                                                     <  1:313136/150‑1 (MQ=255)            ggTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAg                                                                                                                   >  1:164876/1‑148 (MQ=255)               cATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta                                                                                                             <  2:303425/151‑1 (MQ=255)               cATGGTGGCGATGCACTTTCAGTCACCGACTGGTCTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta                                                                                                             <  1:312223/151‑1 (MQ=255)                aTGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggtag                                                                                                            <  1:145125/151‑1 (MQ=255)                      gCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTacaa                                                                                                      >  2:115746/1‑151 (MQ=255)                       cGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGACTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTacaac                                                                                                     <  2:197733/151‑1 (MQ=255)                         aTGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAg                                                                                                   >  1:18271/1‑151 (MQ=255)                              cTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCg                                                                                              >  2:496103/1‑151 (MQ=255)                              cTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCg                                                                                              >  2:365221/1‑151 (MQ=255)                                  cAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGaaa                                                                                          >  1:521160/1‑151 (MQ=255)                                    gTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAAc                                                                                         >  1:253215/1‑150 (MQ=255)                                        ccGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAg                                                                                    <  1:184695/151‑1 (MQ=255)                                              ggTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCag                                                                              >  1:241945/1‑151 (MQ=255)                                               gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  2:369902/1‑151 (MQ=255)                                               gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  1:204104/1‑151 (MQ=255)                                               gtgtGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCaga                                                                             >  1:105589/1‑151 (MQ=255)                                                   ggAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGGAACCAAAGATCCAGAGaaa                                                                         >  2:57650/1‑151 (MQ=255)                                                      aGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCa                                                                      >  1:410936/1‑151 (MQ=255)                                                          tCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                  <  1:367441/151‑1 (MQ=255)                                                          tCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                  <  2:442829/151‑1 (MQ=255)                                                          tCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                  <  1:491432/151‑1 (MQ=255)                                                          tCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGc                                                                  <  2:100017/151‑1 (MQ=255)                                                                           gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTc                                                  <  1:225787/150‑1 (MQ=255)                                                                              tGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCAtt                                              >  1:223271/1‑151 (MQ=255)                                                                                aCGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTc                                             <  2:100748/150‑1 (MQ=255)                                                                                     ggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATc                                        >  2:414786/1‑150 (MQ=255)                                                                                         ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTAt                                   >  1:420205/1‑151 (MQ=255)                                                                                         ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTAt                                   >  2:323772/1‑151 (MQ=255)                                                                                              ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                              <  2:9966/151‑1 (MQ=255)                                                                                              ccGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTa                              <  1:500254/151‑1 (MQ=255)                                                                                                 aaaGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGt                           <  1:181082/151‑1 (MQ=255)                                                                                                 aaaGGTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGt                           <  1:96830/151‑1 (MQ=255)                                                                                                    ggTTTAAGTGGCTTCTTTGCCGTTTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGTAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  1:140738/1‑149 (MQ=255)                                                                                                    ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  1:274551/1‑149 (MQ=255)                                                                                                    ggTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGtt                          >  1:99405/1‑149 (MQ=255)                                                                                                     gTTTAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTc                         <  2:57076/149‑1 (MQ=255)                                                                                                        tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                    <  1:530806/151‑1 (MQ=255)                                                                                                        tAAGTGGCTTCTTTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGt                    <  2:179302/151‑1 (MQ=255)                                                                                                               cttcttTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACC‑AAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCgcgc            <  2:164369/151‑1 (MQ=255)                                                                                                                ttcttTGCCGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCgcgc            <  2:201969/151‑1 (MQ=255)                                                                                                                       ccGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATtg      <  1:53592/150‑1 (MQ=255)                                                                                                                       ccGGATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATtg      <  1:238434/150‑1 (MQ=255)                                                                                                                          gATTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  >  1:198854/1‑151 (MQ=255)                                                                                                                           aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  <  2:385270/150‑1 (MQ=255)                                                                                                                           aTTCCAGGTAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  <  2:480644/150‑1 (MQ=255)                                                                                                                           aTTCCAGGTAGCGGTTTTCGCGTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAt  <  1:31555/150‑1 (MQ=255)                                                                                                                                  |                                                                                                                                               GTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAAAGATCCAGAGAAATCACTGCCACGCGCGATTAACTCCATTCCGATCCGTATCATTATGTTCTACGTCTTCGCGCTGATTGTGAT  >  NC_000913/4430392‑4430663

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