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%	260.7 / NA	66	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 (33/33);  total (33/33)

GCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAA  >  NC_000913/3914942‑3915234                                                                                                                                                    |                                                                                                                                                  gcattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCt                                                                                                                                                 <  1:756166/150‑1 (MQ=255) gcattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTg                                                                                                                                                >  2:33887/1‑151 (MQ=255) gcattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTg                                                                                                                                                <  2:1005221/151‑1 (MQ=255)  cattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTg                                                                                                                                                >  1:360633/1‑150 (MQ=255)   attgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGa                                                                                                                                               >  2:661721/1‑150 (MQ=255)   attgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGa                                                                                                                                               >  1:578507/1‑150 (MQ=255)    ttgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAAc                                                                                                                                             >  2:681851/1‑151 (MQ=255)     tgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAAcc                                                                                                                                            >  1:943145/1‑151 (MQ=255)     tgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAAcc                                                                                                                                            >  1:212781/1‑151 (MQ=255)       catGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCAt                                                                                                                                          <  1:941240/151‑1 (MQ=255)          gACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGc                                                                                                                                       <  2:422718/151‑1 (MQ=255)                  cccAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCgg                                                                                                                               >  1:994930/1‑151 (MQ=255)                   ccAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGc                                                                                                                              <  1:781319/151‑1 (MQ=255)                    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCa                                                                                                                             <  1:213409/151‑1 (MQ=255)                         gctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGt                                                                                                                        <  1:847837/151‑1 (MQ=255)                           tctGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTAt                                                                                                                      >  1:466731/1‑151 (MQ=255)                              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATg                                                                                                                     <  1:305233/149‑1 (MQ=255)                              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATg                                                                                                                     <  1:925303/149‑1 (MQ=255)                                   aaGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACt                                                                                                              >  1:801286/1‑151 (MQ=255)                                     gCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACttt                                                                                                            <  1:266661/151‑1 (MQ=255)                                              tAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGAt                                                                                                    <  1:679113/150‑1 (MQ=255)                                              tAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCACCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGAt                                                                                                    <  1:363757/150‑1 (MQ=255)                                                     gtcgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATa                                                                                             <  2:678566/150‑1 (MQ=255)                                                       cgtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACa                                                                                           <  2:746852/150‑1 (MQ=255)                                                        gtcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATg                                                                                         <  2:6113/151‑1 (MQ=255)                                                         tcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                        >  2:7326/1‑151 (MQ=255)                                                          cTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgcg                                                                                       >  1:56463/1‑151 (MQ=255)                                                              cAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTg                                                                                   >  1:184981/1‑151 (MQ=255)                                                                  gcCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCtt                                                                               <  2:1001194/151‑1 (MQ=255)                                                                   ccTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCttt                                                                              <  1:939726/151‑1 (MQ=255)                                                                   ccTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCttt                                                                              >  2:567558/1‑151 (MQ=255)                                                                   ccTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCttt                                                                              <  1:688378/151‑1 (MQ=255)                                                                    cTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTg                                                                             >  2:16119/1‑151 (MQ=255)                                                                          tAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGc                                                                        <  2:581930/150‑1 (MQ=255)                                                                                  cgccgcCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGa                                                                <  1:789056/150‑1 (MQ=255)                                                                                     cgccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATcgc                                                            <  1:935073/151‑3 (MQ=255)                                                                                     cgccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCAc                                                            >  1:606639/1‑151 (MQ=255)                                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                         >  2:76311/1‑151 (MQ=255)                                                                                           gACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCt                                                      >  2:527340/1‑151 (MQ=255)                                                                                            aCCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTc                                                     >  1:151531/1‑151 (MQ=255)                                                                                             ccGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTc                                                     >  2:714464/1‑150 (MQ=255)                                                                                             ccGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCa                                                    <  2:167870/151‑1 (MQ=255)                                                                                             ccGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCa                                                    <  1:743123/151‑1 (MQ=255)                                                                                                 acacCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGc                                                <  2:987745/151‑1 (MQ=255)                                                                                                    ccAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCatt                                             <  2:725749/151‑1 (MQ=255)                                                                                                    ccAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCatt                                             >  2:990002/1‑151 (MQ=255)                                                                                                    ccAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCatt                                             <  2:345907/151‑1 (MQ=255)                                                                                                          gAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCa                                        >  1:511863/1‑150 (MQ=255)                                                                                                           aaCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCaa                                       >  1:153963/1‑150 (MQ=255)                                                                                                           aaCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCa                                        <  2:342821/149‑1 (MQ=255)                                                                                                            aCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAAc                                      >  1:15328/1‑150 (MQ=255)                                                                                                             cGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACAt                                    <  2:97156/151‑1 (MQ=255)                                                                                                             cGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACAt                                    <  2:885577/151‑1 (MQ=255)                                                                                                              gTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACAt                                    >  1:625159/1‑150 (MQ=255)                                                                                                                    tgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTc                             >  2:802255/1‑151 (MQ=255)                                                                                                                      cgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCt                           <  2:128846/151‑1 (MQ=255)                                                                                                                         cTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGAc                        >  2:488625/1‑151 (MQ=255)                                                                                                                         cTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCAGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGAc                        >  2:574965/1‑151 (MQ=255)                                                                                                                                attcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAtt                 <  1:107815/151‑1 (MQ=255)                                                                                                                                 ttcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAttt                >  2:846699/1‑151 (MQ=255)                                                                                                                                  tcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg               >  2:363699/1‑151 (MQ=255)                                                                                                                                    gcagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg               >  2:549287/1‑149 (MQ=255)                                                                                                                                     cagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg               >  2:940466/1‑148 (MQ=255)                                                                                                                                         aTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACg         >  1:1004736/1‑150 (MQ=255)                                                                                                                                           cAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGaa       <  1:669518/150‑1 (MQ=255)                                                                                                                                               gACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTaa  <  1:777542/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAA  >  NC_000913/3914942‑3915234

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