Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,933	T→C	D22G (GAC→GGC)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,933	0	T	C	100.0%	245.0 / NA	70	D22G (GAC→GGC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/43);  total (27/43)

CACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCG  >  NC_000913/3379785‑3380075                                                                                                                                                     |                                                                                                                                               cACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGtcgccg                                                                                                                                              >  2:670912/1‑151 (MQ=255)      gCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTg                                                                                                                                         <  1:137795/151‑1 (MQ=255)       cGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTg                                                                                                                                         >  2:244311/1‑150 (MQ=255)       cGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTg                                                                                                                                         >  2:700703/1‑150 (MQ=255)        gTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGAt                                                                                                                                       <  1:906826/151‑1 (MQ=255)              tttGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGc                                                                                                                                 >  2:681594/1‑151 (MQ=255)                 gtgtCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCtt                                                                                                                              <  1:365929/151‑1 (MQ=255)                 gtgtCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCt                                                                                                                               <  1:144236/150‑1 (MQ=255)                        cgcgTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTa                                                                                                                         >  2:528376/1‑149 (MQ=255)                         gcgTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAgc                                                                                                                      <  1:97340/151‑1 (MQ=255)                            tttACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGCGGCCTGGCGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGCCGGCCTCTTTTTGAACGCCGCCGGGTTGATGGCTGCCTTCATTAAgcgc                                                                                                                    >  2:645372/1‑150 (MQ=255)                               aCCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTc                                                                                                                >  1:164752/1‑151 (MQ=255)                                   cacaGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAgg                                                                                                             >  2:203661/1‑150 (MQ=255)                                      aGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAtt                                                                                                         >  1:868369/1‑151 (MQ=255)                                       gCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTg                                                                                                        >  1:736852/1‑151 (MQ=255)                                        cTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGc                                                                                                       >  1:836515/1‑151 (MQ=255)                                           aCCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTGCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGc                                                                                                       <  2:673353/148‑1 (MQ=255)                                            ccGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGAt                                                                                                    <  1:252841/150‑1 (MQ=255)                                              gACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGtt                                                                                                 >  2:706281/1‑151 (MQ=255)                                              gACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGtt                                                                                                 <  2:575676/151‑1 (MQ=255)                                                     gcCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTaa                                                                                          <  1:868091/151‑1 (MQ=255)                                                      ccATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAAc                                                                                         <  1:259771/151‑1 (MQ=255)                                                         ttAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCt                                                                                      >  1:559412/1‑151 (MQ=255)                                                          tAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCt                                                                                      <  1:366717/150‑1 (MQ=255)                                                            gCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTg                                                                                   <  1:80405/151‑1 (MQ=255)                                                               cACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAGTGCGATGTTGGGGTAACGCTTTGCAt                                                                                <  1:493175/151‑1 (MQ=255)                                                                 cTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATg                                                                               >  1:487971/1‑150 (MQ=255)                                                                   cGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGTCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATc                                                                            <  1:650495/151‑1 (MQ=255)                                                                   cGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGa                                                                              <  2:754857/149‑1 (MQ=255)                                                                    gTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCt                                                                           <  2:258324/151‑1 (MQ=255)                                                                      ccTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTa                                                                         <  1:618161/151‑1 (MQ=255)                                                                         gggggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTag                                                                        <  2:246447/149‑2 (MQ=255)                                                                         gggggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTAtt                                                                       <  2:559569/150‑1 (MQ=255)                                                                            ggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTa                                                                   <  1:869509/151‑1 (MQ=255)                                                                             gCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTAt                                                                  <  2:495511/151‑1 (MQ=255)                                                                                ggtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTAtt                                                                 <  1:648015/149‑1 (MQ=255)                                                                                ggtggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTc                                                                <  1:493392/150‑1 (MQ=255)                                                                                   ggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCtt                                                             <  2:429746/150‑1 (MQ=255)                                                                                    gCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGa                                                           <  2:903402/151‑1 (MQ=255)                                                                                      gTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaat                                                         <  1:285549/151‑1 (MQ=255)                                                                                        gCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaataa                                                       <  1:698968/151‑1 (MQ=255)                                                                                        gCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaata                                                        >  1:740054/1‑150 (MQ=255)                                                                                        gCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaata                                                        <  2:827364/150‑1 (MQ=255)                                                                                           aTTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCg                                                    <  2:564112/151‑1 (MQ=255)                                                                                               atgaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTg                                                 >  2:791052/1‑150 (MQ=255)                                                                                                   gTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCg                                            >  2:434583/1‑151 (MQ=255)                                                                                                     tCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCgtt                                          >  2:694703/1‑151 (MQ=255)                                                                                                     tCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCgtt                                          >  1:890466/1‑151 (MQ=255)                                                                                                      ccTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCgttg                                         <  1:164322/151‑1 (MQ=255)                                                                                                         gATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCgttgtt                                       <  1:901739/150‑1 (MQ=255)                                                                                                          ataaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCa                                     <  1:439646/151‑1 (MQ=255)                                                                                                            aaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCAc                                    <  2:15077/150‑1 (MQ=255)                                                                                                             aataatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACgg                                  >  2:260060/1‑151 (MQ=255)                                                                                                                aatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGtt                                <  1:57747/150‑1 (MQ=255)                                                                                                                aatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGtt                                <  1:52565/150‑1 (MQ=255)                                                                                                                 atTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGa                              <  2:648880/151‑1 (MQ=255)                                                                                                                    tCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACg                            >  1:365445/1‑150 (MQ=255)                                                                                                                        ggcggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaa                       <  1:160358/151‑1 (MQ=255)                                                                                                                           ggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaaa                      >  1:456128/1‑149 (MQ=255)                                                                                                                             cTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGc                  >  2:782026/1‑151 (MQ=255)                                                                                                                               gACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGcc                 <  2:733532/150‑1 (MQ=255)                                                                                                                               gACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGcc                 <  2:81445/150‑1 (MQ=255)                                                                                                                                aCGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGt               <  1:318881/151‑1 (MQ=255)                                                                                                                                 cGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAATGCCGtt              <  2:577255/151‑1 (MQ=255)                                                                                                                                 cGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGCTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGtt              <  1:680306/151‑1 (MQ=255)                                                                                                                                  ggCCTCTTTTAGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGtt              >  2:618981/1‑150 (MQ=255)                                                                                                                                   gCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCt            >  2:433415/1‑151 (MQ=255)                                                                                                                                   gCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCt            >  1:835349/1‑151 (MQ=255)                                                                                                                                     ctctTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTac          >  1:636763/1‑151 (MQ=255)                                                                                                                                             gAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTAcgc   <  1:140195/150‑1 (MQ=255)                                                                                                                                              aaCGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTAcgcg  >  2:10405/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                               CACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTACGCG  >  NC_000913/3379785‑3380075

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