Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	469,271	G→A	G134E (GGG→GAG)	mdlA →	ABC transporter family protein MdlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	469,271	0	G	A	100.0%	119.5 / NA	40	G134E (GGG→GAG)	mdlA	ABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/21);  total (19/21)

CAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCA  >  NC_000913/469033‑469501                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       cAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCt                                                                                                                                                                                                                            <  1:53694/251‑1 (MQ=255)                  cGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGa                                                                                                                                                                                                           <  1:78046/251‑1 (MQ=255)                                 gatggCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTGCCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTTCCGCCGGAGAAGAGGGGCGGACTCTGGTGGATTCACTGGTGATGTGCTGCGCtgt                                                                                                                                                                                              <  2:29265/245‑1 (MQ=255)                                     gCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGAt                                                                                                                                                                                        <  2:166432/251‑1 (MQ=255)                                        gTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTTAAGCTGGTGGATTCACTGGTGGTGGGCTGCGCTGGGGTGATTat                                                                                                                                                                                     >  1:288894/1‑251 (MQ=255)                                           gtgGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGGGAGGGGCTGCGCTGTGTTGATTatga                                                                                                                                                                                   <  2:191919/250‑1 (MQ=255)                                                ttATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTCGCTGTTGAACTGCGTGAAGCTTATTACCGTCAGCTAACCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGCTCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGGGGTGGGCTGCGGTGTGTTGATTATGATGTCt                                                                                                                                                                              >  1:31650/1‑250 (MQ=255)                                                 tATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGGCGATCGCGTCGTGGTTGCCGCCGGAGGAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAc                                                                                                                                                                            >  2:72965/1‑251 (MQ=255)                                                  aTCTCCTGCGTTACGTCTGGCGGGTGTTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTAc                                                                                                                                                                            <  2:15107/250‑1 (MQ=255)                                                    ctcCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGCGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGCCCAATGAAGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGGGCTGACGCCGGGGGGTTCACTGGGGATGGGCTGCGCTGTGTTGATTTTGATGTCTACgg                                                                                                                                                                          >  2:88238/1‑249 (MQ=255)                                                      ccTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTTCCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTTGTGATGGGCTGCGCTGTGTTGATTattttttctcctcaa                                                                                                                                                                        <  1:197204/250‑14 (MQ=255)                                                           gTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTTGATTCACTGGTGATTGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGc                                                                                                                                                                  <  1:126306/251‑1 (MQ=255)                                                                       gggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGCACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGCGCTGAACCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGAcc                                                                                                                                                      >  1:31357/1‑251 (MQ=255)                                                                        ggTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCt                                                                                                                                                     >  1:99815/1‑251 (MQ=255)                                                                                           tCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGGGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCCAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGAt                                                                                                                                  >  2:85792/1‑251 (MQ=255)                                                                                              tATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGAt                                                                                                                                  <  1:162459/248‑1 (MQ=255)                                                                                               aTCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGcc                                                                                                                               >  2:153959/1‑250 (MQ=255)                                                                                                 cAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGt                                                                                                                            <  1:74600/251‑1 (MQ=255)                                                                                                   aCTGGATGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTgg                                                                                                                          >  2:130560/1‑250 (MQ=255)                                                                                                      ggCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGAt                                                                                                                         <  2:27622/249‑1 (MQ=255)                                                                                                       gCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGc                                                                                                                      >  1:38966/1‑251 (MQ=255)                                                                                                              aaCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGGGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGGGTTGATTATGATGTCTACGAAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATg                                                                                                               >  2:239064/1‑251 (MQ=255)                                                                                                               aCTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGa                                                                                                              <  2:125004/251‑1 (MQ=255)                                                                                                                  gCGTGAAGATTATTACCGTCAGCTAAGCCGACAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTAGCTCGCGTCGTGTTTGCCGCCGGGGAAGAGGCGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATTGCGCACATGAta                                                                                                            >  2:289375/1‑249 (MQ=255)                                                                                                                   cGTGAAGATTAATACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCACGTGCGACAAATGACGTAGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGAAGCTGGTGGATTCACTGGTGATGGGCTGAGCTGTGTTGCTTATGATGTCTACGCAACTTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCTGTGATGGCGATCATGATTac                                                                                                          >  2:131578/1‑250 (MQ=255)                                                                                                                          attattACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCAGGATTAAGCGCAAc                                                                                                   <  1:236103/251‑1 (MQ=255)                                                                                                                             attaCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGc                                                                                                <  2:221289/251‑1 (MQ=255)                                                                                                                               taCCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGACGTGTTTGCCGCCGGAGAAGAGGCGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGGAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGCTCATGATTAAGCGCAACGGCg                                                                                               >  1:188528/1‑250 (MQ=255)                                                                                                                                  cGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGc                                                                                           >  1:178794/1‑251 (MQ=255)                                                                                                                                  cGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATGTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGc                                                                                           <  2:56681/251‑1 (MQ=255)                                                                                                                                           aGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATg                                                                                   >  1:232717/1‑250 (MQ=255)                                                                                                                                              cGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACg                                                                               <  2:259049/251‑1 (MQ=255)                                                                                                                                               ggCAGCATCCTGAGATTTACCTCCGTCATCGCACCGGTGACCTCATGGCTCGTGCGCCCAATGACGACGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCCGGTGGATTCACTGGAGATGGGCTGAGCTGTGTTGATTATTATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCCTGAACGc                                                                              >  2:136965/1‑251 (MQ=255)                                                                                                                                                 cagtcaccTGAGTTTTTCCTGCGACATCGCTCCGGTGACCTCAGGGCTCGGGCGACCAATGACGTCGATTGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTTACGCTGGTGGATTCTCTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGAAGGTGACCTTATTTTCCCTGTTGCCTATGCCAGTGATGGCGATCATGATTAAGCTCAACGTCTATGCTTTGCATGAACGCtt                                                                            <  2:300101/245‑1 (MQ=255)                                                                                                                                                       ccTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCt                                                                      >  1:154820/1‑251 (MQ=255)                                                                                                                                                                             gCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGt                                                <  1:99667/251‑1 (MQ=255)                                                                                                                                                                             gCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGt                                                <  1:85792/251‑1 (MQ=255)                                                                                                                                                                                   gTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                          <  2:132928/251‑1 (MQ=255)                                                                                                                                                                                                                       gcgTCGTGGTTGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCGGTGTTGATTATGATGGCTACGCAAATTAGCTTGCAGTTGACCTTATTTGCCCGGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCTAGGCTGTGCATGAACGCGTTAAGCTGGCACAGGCTGCGTTTTCCAGTTTGAATGACCGCACCCGGGAAAGCCTCATCAGGATCTGCATg      <  2:204981/251‑1 (MQ=255)                                                                                                                                                                                                                            gtttttGCCGCCGGAGAAGAGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCa  >  2:163386/4‑250 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       CAGATCCTGATGTGGATCGCCACCATGGTGCTGATTGCCGTTGTGGTTTATCTCCTGCGTTACGTCTGGCGGGTATTGCTGTTTGGTGCGTCTTATCAACTGGCTGTTGAACTGCGTGAAGATTATTACCGTCAGCTAAGCCGGCAGCATCCTGAGTTTTACCTGCGTCATCGCACCGGTGACCTCATGGCTCGTGCGACCAATGACGTCGATCGCGTCGTGTTTGCCGCCGGAGAAGGGGTGCTGACGCTGGTGGATTCACTGGTGATGGGCTGCGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCA  >  NC_000913/469033‑469501

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