Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,512,984	T→A	6.3%	I670F (ATC→TTC)	ftsK ←	DNA translocase FtsK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,512,984	0	T	A	6.3%	192.2 / 2.5	63	I670F (ATC→TTC)	ftsK	DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGCCGTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGCGAAT  >  NC_002947/4512840‑4513132                                                                                                                                                 |                                                                                                                                                     cgcCTTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCAGGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGtt                                                                                                                                                   <  4:316919/148‑1 (MQ=255) cgcCGTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTg                                                                                                                                                  <  7:285206/149‑1 (MQ=255) cgcCGTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTg                                                                                                                                                  <  5:93077/149‑1 (MQ=255)  gcCGTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTgg                                                                                                                                                 >  4:194156/1‑149 (MQ=255)        cACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGa                                                                                                                                           <  2:285346/149‑1 (MQ=255)          cccGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGAt                                                                                                                                          >  3:142058/1‑148 (MQ=255)          cccGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGAt                                                                                                                                          >  3:212572/1‑148 (MQ=255)             ggATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAgg                                                                                                                                      <  3:194156/149‑1 (MQ=255)                  ggCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCgg                                                                                                                                  <  8:288494/148‑1 (MQ=255)                   gCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGt                                                                                                                                 >  3:245145/1‑148 (MQ=255)                         tGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCgg                                                                                                                                  <  4:73707/141‑1 (MQ=255)                         tGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCgg                                                                                                                                  >  3:73707/1‑141 (MQ=255)                           gTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTc                                                                                                                        >  2:213895/1‑149 (MQ=255)                              cccGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCAc                                                                                                                     >  5:90999/1‑149 (MQ=255)                               ccGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCAcc                                                                                                                    >  8:243692/1‑149 (MQ=255)                                  gcggcATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGTCTTGATCAGGCCGGTGATCACGTCCACCGAc                                                                                                                 >  8:538/1‑149 (MQ=255)                                      cATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGgc                                                                                                             <  7:299066/149‑1 (MQ=255)                                         gctgtacATGCCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGACCGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCAATTCAACCGACGGGCGCt                                                                                                          >  4:79501/7‑149 (MQ=255)                                            cAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTggg                                                                                                       <  4:245145/149‑1 (MQ=255)                                             aGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGt                                                                                                      <  5:206024/149‑1 (MQ=255)                                                   tCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAg                                                                                                <  8:14625/149‑1 (MQ=255)                                                   ccGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTTGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAg                                                                                                <  3:79501/148‑1 (MQ=255)                                                     gCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCCCCGACGGGCGCTGGGTGGCGAGGa                                                                                              >  6:158497/1‑149 (MQ=255)                                                     gCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTACTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGa                                                                                              >  3:330006/1‑149 (MQ=255)                                                      ccATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGAt                                                                                             <  1:313570/149‑1 (MQ=255)                                                        aTGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGTATCa                                                                                           <  7:243692/149‑1 (MQ=255)                                                             ccAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTg                                                                                       >  4:235402/1‑148 (MQ=255)                                                                      gTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGTCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGgcc                                                                             <  7:538/149‑1 (MQ=255)                                                                       ttCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGgc                                                                              <  6:90999/147‑1 (MQ=255)                                                                       ttCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGAAGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACg                                                                                                                >  5:154536/1‑113 (MQ=255)                                                                       ttCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGAAGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACg                                                                                                                >  5:160167/1‑113 (MQ=255)                                                                       ttCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGAAGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACg                                                                                                                <  6:160167/113‑1 (MQ=255)                                                                       ttCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGAAGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACg                                                                                                                <  6:154536/113‑1 (MQ=255)                                                                                  ccTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTc                                                                 >  2:224522/1‑149 (MQ=255)                                                                                      gTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTgc                                                              <  1:351773/148‑1 (MQ=255)                                                                                       tCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTgcgc                                                            <  3:235402/149‑1 (MQ=255)                                                                                         gatgatGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGa                                                          <  1:213895/149‑1 (MQ=255)                                                                                           tgatgGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGa                                                                                              <  7:349582/111‑1 (MQ=255)                                                                                           tgatgGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGa                                                                                              >  8:349582/1‑111 (MQ=255)                                                                                                gTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACggg                                                    <  4:212572/148‑1 (MQ=255)                                                                                                  aCGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGa                                                 >  4:251758/1‑149 (MQ=255)                                                                                                    gcgAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAgg                                                                                               >  3:10953/1‑101 (MQ=255)                                                                                                    gcgAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAgg                                                                                               <  4:10953/101‑1 (MQ=255)                                                                                                      gAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGTCTTCTGCGCGATACGGGCGATCAg                                             >  8:209248/1‑149 (MQ=255)                                                                                                       aGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGc                                            <  4:142058/149‑1 (MQ=255)                                                                                                            aTCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTc                                       >  3:32101/1‑149 (MQ=255)                                                                                                            aTCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTc                                       >  3:361345/1‑149 (MQ=255)                                                                                                             tCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTc                                       >  2:326209/1‑148 (MQ=255)                                                                                                                 gCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCg                                      <  1:224522/145‑1 (MQ=255)                                                                                                                 gCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACct                                  >  4:271513/1‑149 (MQ=255)                                                                                                                   tCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACcttc                                >  4:267779/1‑149 (MQ=255)                                                                                                                   tCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACcttc                                >  1:225867/1‑149 (MQ=255)                                                                                                                     gACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACcttc                                <  5:286279/147‑1 (MQ=255)                                                                                                                          cTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGAc                                    >  5:60815/1‑138 (MQ=255)                                                                                                                          cTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGAc                                    <  6:60815/138‑1 (MQ=255)                                                                                                                          cTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCcga                         >  5:206629/1‑149 (MQ=255)                                                                                                                             gAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCcgac                        <  3:251758/147‑1 (MQ=255)                                                                                                                                    aTACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGatcatcat               <  4:330006/149‑1 (MQ=255)                                                                                                                                      aCGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGAtcatcatca             >  7:278876/1‑149 (MQ=255)                                                                                                                                             gggATGTTGGCCTTGATCAGGCAGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGc      >  4:110081/1‑149 (MQ=255)                                                                                                                                              ggATGTTGGCCTTGATCAGGCCGGTGATCCCGTCCCCCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCCGCTCTTCGCCCTTCTTGCCGACGATCATCATCATGTCGGc      >  5:331473/1‑148 (MQ=255)                                                                                                                                              ggATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCTGCTCTTCGACCTTCTTGCAGACGATCATCATCATGTCGGc      >  5:325424/1‑148 (MQ=255)                                                                                                                                                aTGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGCGaa   <  4:32101/149‑1 (MQ=255)                                                                                                                                                 tGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGCGAAt  <  8:278876/149‑1 (MQ=255)                                                                                                                                                 tGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGCGAAt  <  1:296708/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGCCGTGCACCCGGATCGGCAGACTGGTGCCCGGCGGCATGTACAGCATGTCGCCATGGCCCAGCAGCTGTTCGGCACCGCCCTGGTCGATGATGGTACGCGAGTCGATCTTGCTCGACACCTGGAACGCCATACGAGTCGGGATGTTGGCCTTGATCAGGCCGGTGATCACGTCCACCGACGGGCGCTGGGTGGCGAGGATCAGGTGAATACCGGCCGCCCGTGCCTTCTGCGCGATACGGGCGATCAGCTCTTCGACCTTCTTGCCGACGATCATCATCATGTCGGCGAAT  >  NC_002947/4512840‑4513132

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