Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,267	A→G	D135G (GAC→GGC)	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,267	0	A	G	100.0%	171.0 / NA	46	D135G (GAC→GGC)	cycA	D‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/27);  total (19/27)

AGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTG  >  NC_000913/4430125‑4430412                                                                                                                                               |                                                                                                                                                  aGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTc                                                                                                                                            <  1:160717/150‑1 (MQ=255)   tGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGcc                                                                                                                                         >  2:181503/1‑151 (MQ=255)    gACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCt                                                                                                                                        <  1:295009/151‑1 (MQ=255)       ttCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGc                                                                                                                                     >  2:409222/1‑151 (MQ=255)          gCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTgg                                                                                                                                  <  2:614796/151‑1 (MQ=255)                   ttACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGAt                                                                                                                           >  1:196478/1‑149 (MQ=255)                    tACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAg                                                                                                                         <  1:276336/150‑1 (MQ=255)                     gcTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTgcaggtg                                                                                                                       >  2:254152/2‑145 (MQ=255)                      cTCGGGCCGTTGGCAGGATATTTCACCGGCTGGACTTACTCGTTCTGCTGGGTTGGAACCGGTATGGCAGACGGGGAGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGCCGCCTCGCTGGCggcgt                                                                                                                            >  2:270596/1‑142 (MQ=255)                        cGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGACGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGtgctg                                                                                                                    <  1:741283/151‑1 (MQ=255)                             cGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTgctgctgct                                                                                                               <  2:777966/151‑1 (MQ=255)                                gggCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGa                                                                                                             <  2:823936/150‑1 (MQ=255)                                  gCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGc                                                                                                          <  1:991673/151‑1 (MQ=255)                                   cAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCt                                                                                                         <  1:460703/151‑1 (MQ=255)                                           ttCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCg                                                                                                 >  2:587138/1‑151 (MQ=255)                                            tCACCGGGTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCGCTCCGGCTGGGTCGCCGCGCTGGCGGTGATAGTTCTACTGCTGACGCTCAAtct                                                                                                   <  2:620496/148‑1 (MQ=255)                                            tCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGc                                                                                                                                          >  1:336778/1‑109 (MQ=255)                                            tCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGc                                                                                                                                          <  2:336753/109‑1 (MQ=255)                                                       aCTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGaaaa                                                                                     <  2:143719/151‑1 (MQ=255)                                                          tACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGt                                                                                  >  1:743451/1‑151 (MQ=255)                                                               gTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGt                                                                             <  1:608946/151‑1 (MQ=255)                                                                  ctgctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCgg                                                                              >  2:880770/1‑147 (MQ=255)                                                                         gTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGt                                                                   <  1:181504/151‑1 (MQ=255)                                                                               aCCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTgg                                                              <  1:69094/150‑1 (MQ=255)                                                                                   gTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGGGGTCTGGttt                                                           >  2:1006252/1‑149 (MQ=255)                                                                                       ggCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTGTCCGGGTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGGTGACGCTGAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCgatg                                                     <  2:684143/151‑1 (MQ=255)                                                                                       ggCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCgatg                                                     <  1:484320/151‑1 (MQ=255)                                                                                          aGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCgatgat                                                   >  2:979974/1‑150 (MQ=255)                                                                                              gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTTGGGGAGATGGAGTTCTGGTTTGCGATGATCaaa                                               >  2:852860/1‑150 (MQ=255)                                                                                              gtggtgGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTATGGTTTGCGATGATCaaaa                                              <  2:574934/151‑1 (MQ=255)                                                                                                ggtggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAt                                             <  1:1029406/150‑1 (MQ=255)                                                                                                ggtggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAt                                             <  1:852908/150‑1 (MQ=255)                                                                                                 gtggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcg                                           <  1:384273/151‑1 (MQ=255)                                                                                                   ggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAtcgtc                                         <  2:281014/151‑1 (MQ=255)                                                                                                            ggCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATtgtg                                >  1:523333/1‑151 (MQ=255)                                                                                                              cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTgtgt                               >  1:617549/1‑150 (MQ=255)                                                                                                              cTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTc                              >  1:56275/1‑151 (MQ=255)                                                                                                                 aTGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCt                           <  1:682218/151‑1 (MQ=255)                                                                                                                     tCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTACGGTGAGATGGAGTTTTGGTTTGCGGTGATCAAAATCGTCGCCATTGTGTCGCTGAtt                       >  2:896608/1‑151 (MQ=255)                                                                                                                             ggTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                  <  2:860929/148‑1 (MQ=255)                                                                                                                             ggTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGtcgtcg                <  2:852654/150‑1 (MQ=255)                                                                                                                                 ccccGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCt            <  2:74838/150‑1 (MQ=255)                                                                                                                                   ccGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGGGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTgg          >  2:548423/1‑150 (MQ=255)                                                                                                                                        ctctcCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCAtgg    >  1:641358/1‑151 (MQ=255)                                                                                                                                          ctcCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTc        <  2:707095/145‑1 (MQ=255)                                                                                                                                          ctcCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCAtggtg  >  2:975207/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                  AGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTG  >  NC_000913/4430125‑4430412

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