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%	194.3 / NA	59	D22G (GAC→GGC)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/30);  total (29/30)

CACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTA  >  NC_000913/3379785‑3380071                                                                                                                                                     |                                                                                                                                           cACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGtcgccg                                                                                                                                          >  1:487823/1‑151 (MQ=255)  aCTGGCGTATGCTTTGTGTCTTCTCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGtcgccg                                                                                                                                          <  2:550017/150‑1 (MQ=255)   cTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCggg                                                                                                                                        >  2:576674/1‑151 (MQ=255)          aTGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATgg                                                                                                                                 >  2:363466/1‑151 (MQ=255)              tttGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGACTCTTTTTGAACATGGCCGGGTTGATGGCTGc                                                                                                                             >  1:520035/1‑151 (MQ=255)                tgtgTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCt                                                                                                                           >  2:50258/1‑151 (MQ=255)                    tCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCAt                                                                                                                       >  2:554609/1‑151 (MQ=255)                          cgTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAgcg                                                                                                                 <  1:71762/151‑1 (MQ=255)                          cgTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAgcg                                                                                                                 <  1:127792/151‑1 (MQ=255)                           gTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAgcgc                                                                                                                >  1:762011/1‑151 (MQ=255)                              tACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCt                                                                                                              <  1:180277/150‑1 (MQ=255)                              tACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGATGCCTTCATTAAGCGcc                                                                                                               >  1:601056/1‑149 (MQ=255)                                   cacaGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGt                                                                                                        >  2:303659/1‑151 (MQ=255)                                   cacaGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGt                                                                                                        <  1:827470/151‑1 (MQ=255)                                    acaGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGt                                                                                                        >  2:759798/1‑150 (MQ=255)                                     caGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAt                                                                                                      <  1:630751/151‑1 (MQ=255)                                      aGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAtt                                                                                                     <  2:730810/151‑1 (MQ=255)                                      aGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTAtt                                                                                                     <  2:285098/151‑1 (MQ=255)                                       gCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTg                                                                                                    >  1:123587/1‑151 (MQ=255)                                           aCCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGAt                                                                                                >  2:754744/1‑151 (MQ=255)                                            ccGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATg                                                                                               <  2:389659/151‑1 (MQ=255)                                            ccGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATg                                                                                               <  1:982/151‑1 (MQ=255)                                            ccGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATg                                                                                               <  2:127453/151‑1 (MQ=255)                                              gACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGtt                                                                                             <  2:181689/151‑1 (MQ=255)                                               aCCCGCGCCATTAGCCCACTCGTCCTGGGGGGTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTg                                                                                            <  1:282337/151‑1 (MQ=255)                                               aCCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTg                                                                                            <  1:331/151‑1 (MQ=255)                                                  cgcgCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTgggg                                                                                         >  1:248749/1‑151 (MQ=255)                                                  cgcgCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTgggg                                                                                         <  1:89794/151‑1 (MQ=255)                                                      ccATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAAc                                                                                     <  2:364849/151‑1 (MQ=255)                                                           aGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCtt                                                                                 <  2:157438/150‑1 (MQ=255)                                                                aCTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATg                                                                           >  2:362034/1‑151 (MQ=255)                                                                        tGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTAtt                                                                   >  1:147059/1‑151 (MQ=255)                                                                         gggggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTAt                                                                    <  1:46517/149‑1 (MQ=255)                                                                         gggggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTg                                                                  <  2:294095/151‑1 (MQ=255)                                                                            ggCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTa                                                               <  2:444077/151‑1 (MQ=255)                                                                             gCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTAt                                                              >  2:52212/1‑151 (MQ=255)                                                                                   ggCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTg                                                        <  2:639093/151‑1 (MQ=255)                                                                                     cGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaa                                                      <  2:359323/151‑1 (MQ=255)                                                                                      gTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaat                                                     >  2:247694/1‑151 (MQ=255)                                                                                        gCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaataa                                                   <  1:561622/151‑1 (MQ=255)                                                                                         cTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGaataat                                                  >  1:422061/1‑151 (MQ=255)                                                                                          tATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATc                                                 >  1:448341/1‑151 (MQ=255)                                                                                             tgatgaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCgg                                               <  2:367867/150‑1 (MQ=255)                                                                                               atgaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTgc                                            >  2:145749/1‑151 (MQ=255)                                                                                                 gaGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTgcg                                           <  2:463078/150‑1 (MQ=255)                                                                                                      ccTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCgttg                                     >  1:242726/1‑151 (MQ=255)                                                                                                        tGATAAATAATTTCAAGGCGGCTGACGGCATCTTTTTGAACGTAGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCCGGTATTGCGATGATGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGACgttgtt                                   >  1:700950/1‑151 (MQ=255)                                                                                                         gATAAATAATTTCCAGGCGGCTGACGGCCTCTTGTTGAACGGCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTc                                  <  1:420965/151‑1 (MQ=255)                                                                                                            aaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCAc                                <  2:193791/150‑1 (MQ=255)                                                                                                            aaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACg                               <  1:52190/151‑1 (MQ=255)                                                                                                            aaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACg                               >  1:843652/1‑151 (MQ=255)                                                                                                            aaataaTTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACg                               >  2:97696/1‑151 (MQ=255)                                                                                                             aataatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACgg                              >  1:502054/1‑151 (MQ=255)                                                                                                              ataatTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGt                             >  1:340620/1‑151 (MQ=255)                                                                                                                    tCCAGGCGGCTGACGGCCTCTTTTTGAACGGCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGATTGTCg                        <  1:69232/150‑1 (MQ=255)                                                                                                                          cggcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGaaaaaa                  >  1:713304/1‑150 (MQ=255)                                                                                                                            gcTGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAgg               <  2:344258/151‑1 (MQ=255)                                                                                                                              tGACGGCCTCTTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGcc             >  2:294559/1‑151 (MQ=255)                                                                                                                                     ctctTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTac      >  2:689347/1‑151 (MQ=255)                                                                                                                                      tctTTTTGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTac      >  1:419278/1‑150 (MQ=255)                                                                                                                                          tttGAACGTCGCCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTa  <  2:63638/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           CACTGGCGTATGCTTTGTGTCTTCGCGTTTACCCCACAGCTTACCGACCCGCGCCATTAGCCCACTCGTCCTGGGGGCTGGTGGCGTGCTATTGATGAGTTCCTGATAAATAATTTCCAGGCGGCTGACGGCCTCTTTTTGAACGTCGTCGGGTTGATGGCTGCCTTCATTAAGCGCCTTCAGGTATTGCGATGTTGGGGTAACGCTTTGCATGATCTTATTGTTATTCCTTGAATAATCGGTGCGCCGTTGTTCACGGTTGACGAAAAAAAGGCCGTTCTACACTA  >  NC_000913/3379785‑3380071

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