Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	94,278	G→A	A371A (GCG→GCA)	murE →	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	94,278	0	G	A	100.0%	148.7 / NA	46	A371A (GCG→GCA)	murE	UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

TGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGACTTGCGCCGTTATG  >  NC_000913/94062‑94524                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       ggAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGgtgtg                                                                                                                                                                                                                      <  1:45963/250‑1 (MQ=255)  gAAAGCCATCTGATGGGAGATTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTCTCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGCGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGACAGCTGTGgtgtgt                                                                                                                                                                                                                     >  2:314354/1‑251 (MQ=255)                  gcgcTTTTAACGTCAGCACCCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGc                                                                                                                                                                                                     <  2:131459/251‑1 (MQ=255)                   cgcTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTAGCCACGGGCTGACCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTAGGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGCGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCg                                                                                                                                                                                                    <  2:133344/251‑1 (MQ=255)                    gctttaaacgtcctcaaaatgctgcTCGCGCTGGCGCCACTGTTGGCACTCGGCTATCCCCTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCCACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGGGGTGTGTCTTTGGCCGTGGTGGCg                                                                                                                                                                                                    >  2:198471/20‑250 (MQ=255)                         tAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATcgcg                                                                                                                                                                                              >  2:143538/1‑251 (MQ=255)                            cGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGAt                                                                                                                                                                                            >  1:161626/1‑250 (MQ=255)                            cGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATa                                                                                                                                                                                           <  2:18657/251‑1 (MQ=255)                            cGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGACTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGAt                                                                                                                                                                                            >  2:190107/1‑250 (MQ=255)                             gTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTAGGCTATCCACTGGCTGATCTGCTGAAACCCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATa                                                                                                                                                                                           >  1:303712/1‑250 (MQ=255)                                     ccTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGc                                                                                                                                                                                  <  1:285700/251‑1 (MQ=255)                                        gctgctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATGAAGGTAAGCGTc                                                                                                                                                                               >  2:305964/1‑251 (MQ=255)                                        gctgctCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGGGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTc                                                                                                                                                                               >  2:8714/1‑251 (MQ=255)                                         ccgttCGCGCTGGCGACACTGTTGGCCCTCGGCTCTCCACTGGCTGACCTGCGGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTcc                                                                                                                                                                              <  2:303712/247‑1 (MQ=255)                                             tCGCGCTGGCGACACTGTTTGCACTCGGCTATCCACTGGCTGATCTGCGGACAACCGCCGCGCGTCTGCAACAGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATGACGCGCATACGCAGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCAATGTGCGGGGAAGGTGAGGTGTGGGTTTGGCTGTGGTGGCGCTCGAGATAAAGGTCAGCGTGCACTg                                                                                                                                                                          >  2:107774/1‑251 (MQ=255)                                                   gggCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAGCCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGgc                                                                                                                                                                    <  2:47479/250‑1 (MQ=255)                                                     gCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCTCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGgcgc                                                                                                                                                                  <  2:225636/251‑1 (MQ=255)                                                     gCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGgcgc                                                                                                                                                                  <  1:188796/251‑1 (MQ=255)                                                           atGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGc                                                                                                                                                            <  2:225828/250‑1 (MQ=255)                                                                  cACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAag                                                                                                                                                       >  1:99050/1‑249 (MQ=255)                                                                        gCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCt                                                                                                                                               >  1:79965/1‑251 (MQ=255)                                                                           agccACTGGCTGATCTGCTGAAACCCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGa                                                                                                                                             <  2:226339/248‑1 (MQ=255)                                                                                                    cgccgcGCGTATGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGGGGGGACGGACGATacac                                                                                                                   >  1:58501/1‑248 (MQ=255)                                                                                                      ccgcGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGGCGGACGATaaaccc                                                                                                                 >  2:111390/1‑247 (MQ=255)                                                                                                       cgcgcgTCTGCAACAGGTTTGCGGACGTATGGAAATGTTCACTGCGCCAGGCAAACCGACGGCGGTGGCGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGGTCGCGATAAAGGACAGCGTCCACTGATGGGCGCAATTGCCGACGCGCTTGCTGACGTGGCggcgggggcgggagaacaccagc                                                                                                                >  2:308371/1‑230 (MQ=255)                                                                                                             tCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGGAGGGCGATAACCCGCGTACCg                                                                                                          >  2:137893/1‑251 (MQ=255)                                                                                                               ggCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCgaa                                                                                                        <  2:5496/250‑1 (MQ=255)                                                                                                                gCAACCGGGTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCgaag                                                                                                       <  1:49335/251‑1 (MQ=255)                                                                                                                   aCCGGTTTGCGGACGTATGGCAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGGGACGGGCGATAACCCCCGTCCCCAAGAAc                                                                                                    >  2:50305/1‑251 (MQ=255)                                                                                                                      ggTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGGAGAACcgc                                                                                                 >  1:35142/1‑251 (MQ=255)                                                                                                                         ttGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGt                                                                                               <  1:264323/250‑1 (MQ=255)                                                                                                                             ggACGTATGGAAGTGTGCCCTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTTCcat                                                                                          <  2:9614/251‑1 (MQ=255)                                                                                                                               aCGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACTATAACCCGCGTACCGAAGAACCGCGTGCcatca                                                                                        <  1:180048/251‑1 (MQ=255)                                                                                                                                   aTGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGGGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGAGGTGGGGACGGGCGGTAACCCGCCAACCgaagaaaaccgggccctcacca                                                                                     >  2:234416/1‑234 (MQ=255)                                                                                                                                      gAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGa                                                                                 <  2:58501/251‑1 (MQ=255)                                                                                                                                           gtTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTc                                                                            <  2:101619/251‑1 (MQ=255)                                                                                                                                                           aaaCCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCACGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGGTAAAACGCGTCCCGCAGAACCGCGTGCCATCATCAACGCCATTCTGGCGGGAATgatgat                                                            >  2:137062/1‑246 (MQ=255)                                                                                                                                                            aaCCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGTAATGTTAGAt                                                           <  1:78467/251‑1 (MQ=255)                                                                                                                                                                            ggATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTTGCGGGAATGTTAGGTGCCGGACATGCCaaa                                            >  1:86438/1‑250 (MQ=255)                                                                                                                                                                                 aCGCGCAGACGCCGGAGGCACTGGAAAAATCCTTACCGGCAGCGCGTCTGCCCTGTGCGGGCAGGCGTTTGTGTGGCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTTACTGACGAGACCCCTCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCTGGAATGTTAGATGCCCGACATGCCAAAGTGATg                                      <  2:284414/251‑1 (MQ=255)                                                                                                                                                                                           gCCGGATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTg                            <  2:286274/251‑1 (MQ=255)                                                                                                                                                                                              ggATGCACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTg                         <  1:55622/251‑1 (MQ=255)                                                                                                                                                                                                   cACTGGAAAAAGCCTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGGATGTTTGAAGCCGGGCATGCCAAAGTGATGGAAGGGCGTGGTGAAGCg                    >  1:87584/1‑251 (MQ=255)                                                                                                                                                                                                    aCTGGAAAAAGCCTTACAGGCAGCGCGGCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCgg                   <  2:18504/251‑1 (MQ=255)                                                                                                                                                                                                                cTTACAGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGGGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGGATGTTAAATGCCGGACATGCCAAAGTGAAGGAAGGCCGGGCTGAAGCGGTGACTTGCGccc       >  1:117847/1‑250 (MQ=255)                                                                                                                                                                                                                     aGGCAGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATAGTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGGGATGGAAGGCCGTGCTGATGCGGTGACTTGCTCCGTTATg  <  2:182293/251‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       TGAAAGCCATCTGATGGGCGCTTTTAACGTCAGCAACCTGCTGCTCGCGCTGGCGACACTGTTGGCACTCGGCTATCCACTGGCTGATCTGCTGAAAACCGCCGCGCGTCTGCAACCGGTTTGCGGACGTATGGAAGTGTTCACTGCGCCAGGCAAACCGACGGTGGTGGTGGATTACGCGCATACGCCGGATGCACTGGAAAAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGACGATAACCCGCGTACCGAAGAACCGCGTGCCATCATCAACGATATTCTGGCGGGAATGTTAGATGCCGGACATGCCAAAGTGATGGAAGGCCGTGCTGAAGCGGTGACTTGCGCCGTTATG  >  NC_000913/94062‑94524

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