Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,454	T→G	Y108D (TAT→GAT)	wecB →	UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,970,454	0	T	G	100.0%	275.6 / NA	73	Y108D (TAT→GAT)	wecB	UDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (38/35);  total (38/35)

CAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGC  >  NC_000913/3970306‑3970601                                                                                                                                                     |                                                                                                                                                    cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGAt                                                                                                                                                   <  1:451810/151‑1 (MQ=255) cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGAt                                                                                                                                                   >  1:356541/1‑151 (MQ=255)  aaCATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATc                                                                                                                                                  <  1:229708/151‑1 (MQ=255)   aCATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATc                                                                                                                                                  >  1:593011/1‑150 (MQ=255)    cATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGACTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATc                                                                                                                                                  >  1:675813/1‑149 (MQ=255)     aTAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCCGc                                                                                                                                               >  1:898059/1‑151 (MQ=255)        aTGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTa                                                                                                                                            <  2:476149/151‑1 (MQ=255)               cAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGt                                                                                                                                     >  2:97829/1‑151 (MQ=255)                aGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGtt                                                                                                                                    <  2:901131/151‑1 (MQ=255)                  gACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTgg                                                                                                                                  >  1:233569/1‑151 (MQ=255)                   aCAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGt                                                                                                                                 >  2:296322/1‑151 (MQ=255)                   aCAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACAGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGt                                                                                                                                 >  2:679756/1‑151 (MQ=255)                    cAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTc                                                                                                                                >  1:837001/1‑151 (MQ=255)                       ggCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACg                                                                                                                             <  1:705004/151‑1 (MQ=255)                        gCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGt                                                                                                                            <  2:772577/151‑1 (MQ=255)                         ccTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGtt                                                                                                                           >  1:387578/1‑151 (MQ=255)                         ccTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGt                                                                                                                            >  1:498628/1‑150 (MQ=255)                             aCAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAgg                                                                                                                       >  2:59111/1‑151 (MQ=255)                             aCAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAgg                                                                                                                       >  2:249979/1‑151 (MQ=255)                             aCAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAgg                                                                                                                       <  2:761233/151‑1 (MQ=255)                                gagaTAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTg                                                                                                                    <  2:545390/151‑1 (MQ=255)                                  gaTAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGt                                                                                                                  <  1:104988/151‑1 (MQ=255)                                      aCCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTgc                                                                                                              <  2:73868/151‑1 (MQ=255)                                        cTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTgcgc                                                                                                            >  1:679117/1‑151 (MQ=255)                                            cGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACg                                                                                                         <  1:146764/150‑1 (MQ=255)                                              gATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                                                                      <  1:551928/151‑1 (MQ=255)                                              gATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                                                                      <  1:303702/151‑1 (MQ=255)                                                 tCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGa                                                                                                    <  1:576848/150‑1 (MQ=255)                                                   tGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGAtct                                                                                                 >  1:243361/1‑151 (MQ=255)                                                    ggAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGAtctc                                                                                                <  1:561562/151‑1 (MQ=255)                                                        gggCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTAtt                                                                                            <  2:756094/151‑1 (MQ=255)                                                           cTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGc                                                                                         >  2:474476/1‑151 (MQ=255)                                                           cTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGc                                                                                         >  1:338540/1‑151 (MQ=255)                                                             aaaaCCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCg                                                                                       >  2:508880/1‑151 (MQ=255)                                                                aCCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTgg                                                                                    >  1:885569/1‑151 (MQ=255)                                                                    aTTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCgg                                                                                >  2:295808/1‑151 (MQ=255)                                                                     ttCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCgg                                                                                >  2:508259/1‑150 (MQ=255)                                                                           ccGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGc                                                                         >  1:541930/1‑151 (MQ=255)                                                                            cGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCt                                                                        <  2:110880/151‑1 (MQ=255)                                                                             gAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTa                                                                       <  2:80073/151‑1 (MQ=255)                                                                                   aaaCCAGACGTCGTGCTGGTTCACGGCGATACGACCACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGATGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCg                                                                   >  1:368614/1‑149 (MQ=255)                                                                                     aCCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACa                                                               <  1:828216/151‑1 (MQ=255)                                                                                          aCGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGAc                                                          >  2:924548/1‑151 (MQ=255)                                                                                           cgtcgtGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGAcc                                                         >  2:612146/1‑151 (MQ=255)                                                                                           cgtcgtGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGAcc                                                         >  2:327076/1‑151 (MQ=255)                                                                                           cgtcgtGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGAc                                                          <  2:754698/150‑1 (MQ=255)                                                                                              cgtGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCggg                                                      <  2:837517/151‑1 (MQ=255)                                                                                                     gTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTgg                                               >  1:905881/1‑151 (MQ=255)                                                                                                        cgcGGCGATACGACGACGACGCTGGCTACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGGTGGTCGCGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGGCGGAAGATGCTAACCGTACATTGACCGGGCATCTGGc                                              <  2:276345/147‑1 (MQ=255)                                                                                                         aCGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGAt                                           >  1:441196/1‑151 (MQ=255)                                                                                                           ggCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATg                                          >  2:98333/1‑150 (MQ=255)                                                                                                             cGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTAt                                       <  1:116903/151‑1 (MQ=255)                                                                                                              gATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATc                                      <  1:541952/151‑1 (MQ=255)                                                                                                                  cgacgacgacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACtt                                  <  1:669799/151‑1 (MQ=255)                                                                                                                       acgacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTtctctc                             <  2:120385/151‑1 (MQ=255)                                                                                                                       acgacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTtctctc                             >  1:272035/1‑151 (MQ=255)                                                                                                                         gacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCa                           <  2:882761/151‑1 (MQ=255)                                                                                                                         gacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCa                           <  2:803636/151‑1 (MQ=255)                                                                                                                         gacgCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCa                           <  2:801927/151‑1 (MQ=255)                                                                                                                            gCTGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAAcc                        <  1:297758/151‑1 (MQ=255)                                                                                                                              tGGCAACCAGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCg                       >  1:735478/1‑150 (MQ=255)                                                                                                                                      aGCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTccc              >  1:174325/1‑151 (MQ=255)                                                                                                                                       gCCTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCg             <  2:97640/151‑1 (MQ=255)                                                                                                                                        ccTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCgg            >  2:731315/1‑151 (MQ=255)                                                                                                                                        ccTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCgg            >  2:71676/1‑151 (MQ=255)                                                                                                                                        ccTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCgg            <  2:279561/151‑1 (MQ=255)                                                                                                                                         cTGGCGGCGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGc           >  1:333097/1‑151 (MQ=255)                                                                                                                                           ggcggcGTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGc           <  2:61899/149‑1 (MQ=255)                                                                                                                                            gcggcgTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCaa         >  2:33541/1‑150 (MQ=255)                                                                                                                                              ggcgTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCaaaa       <  1:679773/150‑1 (MQ=255)                                                                                                                                              ggcgTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAAc      >  1:714980/1‑151 (MQ=255)                                                                                                                                               gcgTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACt     <  2:338532/151‑1 (MQ=255)                                                                                                                                                 gTTTGATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAACACTtg   >  2:402703/1‑151 (MQ=255)                                                                                                                                                  tttgATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTtgc  >  1:748688/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCGATGTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGC  >  NC_000913/3970306‑3970601

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