Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,915,089	G→T	H38N (CAT→AAT)	glmU ←	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,915,089	0	G	T	100.0%	244.4 / NA	62	H38N (CAT→AAT)	glmU	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (28/34);  total (28/34)

GCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAC  >  NC_000913/3914942‑3915236                                                                                                                                                    |                                                                                                                                                    gcattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTg                                                                                                                                                  <  2:855229/151‑1 (MQ=255)  cattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGa                                                                                                                                                 <  1:246347/151‑1 (MQ=255)   attgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGaa                                                                                                                                                <  1:359003/151‑1 (MQ=255)     tgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAAcc                                                                                                                                              >  1:666585/1‑151 (MQ=255)        atGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATc                                                                                                                                           <  2:279152/151‑1 (MQ=255)           aCCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCt                                                                                                                                        <  1:313081/151‑1 (MQ=255)               gTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTcc                                                                                                                                    >  1:100237/1‑151 (MQ=255)                     aGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCa                                                                                                                               >  1:416571/1‑150 (MQ=255)                        tgctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAggg                                                                                                                           >  1:972424/1‑151 (MQ=255)                        tgctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATCCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAggg                                                                                                                           >  1:754106/1‑151 (MQ=255)                            ctGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTa                                                                                                                         <  2:794406/149‑1 (MQ=255)                              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATgc                                                                                                                      >  2:1010382/1‑150 (MQ=255)                               ccTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATgc                                                                                                                      >  2:51843/1‑149 (MQ=255)                                  gAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCAc                                                                                                                 <  1:194513/151‑1 (MQ=255)                                   aaGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACt                                                                                                                >  1:595848/1‑151 (MQ=255)                                    aGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACtt                                                                                                               <  1:317561/151‑1 (MQ=255)                                             ttAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGa                                                                                                          >  1:859301/1‑147 (MQ=255)                                                aGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCgg                                                                                                   <  2:331244/151‑1 (MQ=255)                                                 ggTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGa                                                                                                          >  1:355516/1‑143 (MQ=255)                                                     gtcgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATa                                                                                               <  1:562772/150‑1 (MQ=255)                                                     gtcgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATa                                                                                               <  2:803088/150‑1 (MQ=255)                                                       cgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACAt                                                                                            >  2:799362/1‑151 (MQ=255)                                                       cgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACAt                                                                                            <  2:959976/151‑1 (MQ=255)                                                         tcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                          >  1:66295/1‑151 (MQ=255)                                                         tcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                          >  2:375301/1‑151 (MQ=255)                                                          cTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGACGATCGGAATACATgc                                                                                          >  2:997374/1‑150 (MQ=255)                                                            ttCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                          <  2:473343/148‑1 (MQ=255)                                                                gcgcCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGt                                                                                      <  1:834396/148‑1 (MQ=255)                                                                  gcCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCtt                                                                                 >  1:200857/1‑151 (MQ=255)                                                                   ccTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCttt                                                                                <  2:840709/151‑1 (MQ=255)                                                                   ccTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCttt                                                                                <  2:972523/151‑1 (MQ=255)                                                                           aGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCgg                                                                        <  2:161423/151‑1 (MQ=255)                                                                             cAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCa                                                                      <  2:989212/151‑1 (MQ=255)                                                                               gATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCa                                                                      >  2:406779/1‑149 (MQ=255)                                                                                aTCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                                   <  2:759493/151‑1 (MQ=255)                                                                                 tCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                                   >  2:13538/1‑150 (MQ=255)                                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                           >  1:504132/1‑151 (MQ=255)                                                                                          tGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGc                                                         >  1:486976/1‑151 (MQ=255)                                                                                           gACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCt                                                        <  2:595903/151‑1 (MQ=255)                                                                                               gTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCAt                                                     >  1:517969/1‑150 (MQ=255)                                                                                                 acacCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAg                                                   <  1:22283/150‑1 (MQ=255)                                                                                                  cacCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCa                                                 <  1:212765/151‑1 (MQ=255)                                                                                                             cGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACAt                                      <  1:888652/151‑1 (MQ=255)                                                                                                              gTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATa                                     <  1:987940/151‑1 (MQ=255)                                                                                                                gcgctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATAcg                                   <  2:59347/151‑1 (MQ=255)                                                                                                                gcgctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATAc                                    <  2:390421/150‑1 (MQ=255)                                                                                                                 cgctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATAcgc                                  <  2:756061/151‑1 (MQ=255)                                                                                                                   ctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGt                                >  2:99119/1‑151 (MQ=255)                                                                                                                   ctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGt                                >  2:825217/1‑151 (MQ=255)                                                                                                                    tgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTc                               >  1:70835/1‑151 (MQ=255)                                                                                                                     gcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTcc                              >  1:771349/1‑151 (MQ=255)                                                                                                                      cgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCt                             <  2:878408/151‑1 (MQ=255)                                                                                                                       gcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTg                            <  2:307063/151‑1 (MQ=255)                                                                                                                       gcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTg                            <  1:722795/151‑1 (MQ=255)                                                                                                                           aaTTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACt                         <  2:200869/150‑1 (MQ=255)                                                                                                                            attcattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGt                       <  2:747897/151‑1 (MQ=255)                                                                                                                               cattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTaa                     <  2:238627/150‑1 (MQ=255)                                                                                                                               cattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAt                    >  1:273637/1‑151 (MQ=255)                                                                                                                                 ttcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAttt                  >  1:839964/1‑151 (MQ=255)                                                                                                                                   cGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTga                <  2:644434/151‑1 (MQ=255)                                                                                                                                                 cATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAc  >  2:957901/1‑151 (MQ=255)                                                                                                                                                 cATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCAGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAc  >  1:180586/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAC  >  NC_000913/3914942‑3915236

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