Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,546,804	T→A	D44E (GAT→GAA)	murP →	N‑acetylmuramic acid‑specific PTS enzyme IICB component/anhydro‑N‑acetylmuramic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,546,804	0	T	A	100.0%	134.4 / NA	40	D44E (GAT→GAA)	murP	N‑acetylmuramic acid‑specific PTS enzyme IICB component/anhydro‑N‑acetylmuramic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (20/20);  total (20/20)

AAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGAT  >  NC_000913/2546566‑2547046                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   aaaCGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGAGCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCaaa                                                                                                                                                                                                                                        <  1:316070/251‑1 (MQ=255)           tgtgATGGTGCTGAACAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCCGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGa                                                                                                                                                                                                                                >  2:234639/1‑249 (MQ=255)            gtgATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAgg                                                                                                                                                                                                                             >  1:76428/1‑251 (MQ=255)                       gAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTAGTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTGCACTGGTTGAACCCAATATCAAAACTCTGGAAGGCGTGAAGGGCg                                                                                                                                                                                                                  <  1:299472/251‑1 (MQ=255)                                 gATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACCGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGCCGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGCCGTCATTTTGAc                                                                                                                                                                                                        >  2:256993/1‑251 (MQ=255)                                        cAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGAc                                                                                                                                                                                                 >  1:271517/1‑251 (MQ=255)                                            aGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAgg                                                                                                                                                                                             >  2:136424/1‑251 (MQ=255)                                               tAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGc                                                                                                                                                                                          >  1:413731/1‑251 (MQ=255)                                                  aaaaCGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGCCAAGGAATAACCCATGGCCAAAGAGATCAGAAGTGAACTTCTGAACACCATTCTTACCAGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTCACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACCGGTTGAACCAAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGGGCAgg                                                                                                                                                                                       >  1:335155/1‑251 (MQ=255)                                                     aCGCCTGGATCAACTCGGCGGCTTTATTCTTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGGCTGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGGTCACTGGGTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCg                                                                                                                                                                                    <  1:525026/251‑1 (MQ=255)                                                           ggATCAACACGGCGGCTTTATGCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGttttt                                                                                                                                                                               <  2:271517/250‑1 (MQ=255)                                                                aaCACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCt                                                                                                                                                                         <  2:141405/251‑1 (MQ=255)                                                                                                   gAATAACCCATTGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcga                                                                                                                                      <  1:119670/251‑1 (MQ=255)                                                                                                      tAACCCATGGCCAACGAGATCAGCCGTGAACTTCTGCACACCATTTTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAACCGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACAGGGAGCAGGTCGCTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCCAgct                                                                                                                                   >  2:290751/1‑251 (MQ=255)                                                                                                           cATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTgg                                                                                                                              >  1:383578/1‑251 (MQ=255)                                                                                                                   aaGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGcc                                                                                                                      >  2:208718/1‑250 (MQ=255)                                                                                                                    agagagCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCTGGAAATATCGCCAGGTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCa                                                                                                                      <  2:428310/250‑1 (MQ=255)                                                                                                                    agagATCAGCAGTGGACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAgaa                                                                                                                      >  1:428310/1‑248 (MQ=255)                                                                                                                     gagaTCAGCAGAGAACTTCTGAAACCCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCAATGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCACAAGCGATGCGCGAGCTGCTGGGAGAAGCAcc                                                                                                                    >  2:57289/1‑251 (MQ=255)                                                                                                                      agaTCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCg                                                                                                                   <  2:6610/251‑1 (MQ=255)                                                                                                                                         gAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGaa                                                                                                 >  1:204534/1‑250 (MQ=255)                                                                                                                                         gAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGaa                                                                                                 >  1:480937/1‑250 (MQ=255)                                                                                                                                            cacCATACTTACCCGGGTCGGCCGACCGGGAAATACCGCCTGTTGGGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGCTGCACCCAATATCAAAACGCTGGGAGGGGTGAAGGGAGTCATTGTGACCAGCGACCAGGTGCAGGTCGTATTTGGGGCTGGTAAAGCTCATCGGGCCGCAAAAGCGATGAGCGAGCTGCTGATAGAAGCACCGGTACAGGATGCCGCAGAAATcg                                                                                             <  2:203660/251‑1 (MQ=255)                                                                                                                                                    ttACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGccc                                                                                       >  1:264544/1‑249 (MQ=255)                                                                                                                                                                  ggACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGtt                                                                       >  2:152219/1‑251 (MQ=255)                                                                                                                                                                            atatCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCcaaa                                                             <  2:220585/251‑1 (MQ=255)                                                                                                                                                                                         tgGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCCGTTAAAAGCCAAACAAACCTCCGGCg                                                >  2:503866/1‑251 (MQ=255)                                                                                                                                                                                             aaCTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGc                                             <  1:348541/250‑1 (MQ=255)                                                                                                                                                                                             aaCTGTATGACGCCCCTGCGTCTGTGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCAGTTTGACCAGCGTCCAGGTGCAGGTCGTCTTTTGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAACCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGca                                            <  2:335155/251‑1 (MQ=255)                                                                                                                                                                                                 gTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGcaac                                          <  1:234639/249‑1 (MQ=255)                                                                                                                                                                                                        gcgcCTGCGTCTGGGTGTACATGACAGTTCACTGTTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGAATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTATCg                                 >  2:290960/1‑251 (MQ=255)                                                                                                                                                                                                        gcgcCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCg                                 <  2:413731/251‑1 (MQ=255)                                                                                                                                                                                                                    gggTGTACATGACAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCGCCGGCGTGCAACAATTTCTCGCCAAATGCGcca                     <  2:256087/251‑1 (MQ=255)                                                                                                                                                                                                                            aTGACAGTTCACTGGTTGAACCCAATAGCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGTCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGGACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTc             <  1:537803/251‑1 (MQ=255)                                                                                                                                                                                                                              gACAGTTCGCTGGTTGAACCCAATATCCAAACGCTGGGAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCCAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCAc           <  1:436720/251‑1 (MQ=255)                                                                                                                                                                                                                               aCAGTTCACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCAcg          >  1:270517/1‑251 (MQ=255)                                                                                                                                                                                                                                  gTTCACTGGTTGCACCCAATATCAAAACGCTGGAAGGCGTGAAGGCCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGAGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTCAAAGCCAAACAAACCTCCGGCGGGCAACAATTTCTCGCCAAATTCGCCACCATCTTCAcgccc       >  1:179867/1‑250 (MQ=255)                                                                                                                                                                                                                                     cACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAACCCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCt     <  1:158008/250‑1 (MQ=255)                                                                                                                                                                                                                                     cACTGGTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTg    <  1:383651/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggTTGAACCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGAt  <  1:136424/249‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   AAACGGCCATTGTGATGGTGCTGAAAAATCTCGATGCCGCAGAAGCTAAAAAACGCCTGGATCAACACGGCGGCTTTATTCGTCAGGTTTTAGACAAGGAATAACCCATGGCCAAAGAGATCAGCAGTGAACTTCTGAACACCATTCTTACCCGTGTCGGCGGACCGGGAAATATCGCCAGTTGTGGTAACTGTATGACGCGCCTGCGTCTGGGTGTACATGACAGTTCACTGGTTGATCCCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAGGATGCCGCAGAAATCGCCGCCCAGAACAAACGTCAGTTAAAAGCCAAACAAACCTCCGGCGTGCAACAATTTCTCGCCAAATTCGCCACCATCTTCACGCCGCTGAT  >  NC_000913/2546566‑2547046

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