Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,556	T→G	6.1%	K213Q (AAG→CAG)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,556	0	T	G	6.1%	205.9 / 2.8	65	K213Q (AAG→CAG)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base G (2/2);  total (29/36)
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.92e-01

GTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAC  >  NC_002947/626411‑626703                                                                                                                                                  |                                                                                                                                                    gTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGc                                                                                                                                                  >  8:224518/1‑149 (MQ=255)     cTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTAAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTa                                                                                                                                              <  7:333786/149‑1 (MQ=255)         tGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGc                                                                                                                                          <  4:211429/149‑1 (MQ=255)            cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                       <  8:269098/149‑1 (MQ=255)            cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                       <  5:147885/149‑1 (MQ=255)             gTCCAGCCTTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCACCAGGCGCTGAATGACCACATGCCGTGCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCAAg                                                                                                                                      <  2:250486/149‑1 (MQ=255)                cAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAGgct                                                                                                                                   <  7:224518/149‑2 (MQ=255)                  gCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATGGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAGgctca                                                                                                                                 <  7:338497/149‑4 (MQ=255)                     ttGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCAc                                                                                                                              >  6:313307/1‑149 (MQ=255)                         cAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgcag                                                                                                                          <  8:356880/149‑1 (MQ=255)                           gCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGTCGTACGGCCCAGGCGCCCACgcagc                                                                                                                         >  2:125636/1‑148 (MQ=255)                            ccGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCaa                                                                                                                       <  8:6352/149‑1 (MQ=255)                                  gcgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCga                                                                                                                  <  5:306262/148‑1 (MQ=255)                                   cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCgaga                                                                                                                <  1:147117/149‑1 (MQ=255)                                   cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCgaga                                                                                                                <  1:300154/149‑1 (MQ=255)                                         ggcttccAATGCTGTGGCGTGATCGGGTAACGAGCCAGACTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGATCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                          <  6:48672/144‑1 (MQ=255)                                         ggcTGCCAATTCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                          <  6:53965/147‑1 (MQ=255)                                              ccATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                     >  2:293869/1‑149 (MQ=255)                                               cATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                     >  2:22374/1‑148 (MQ=255)                                               cATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                     >  2:28185/1‑148 (MQ=255)                                                aTTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCa                                                                                                   <  3:132672/149‑1 (MQ=255)                                                aTTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCa                                                                                                   >  3:334955/1‑149 (MQ=255)                                                             gATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgat                                                                                      >  4:75295/1‑149 (MQ=255)                                                                 gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCTTCGGCACCAGCGACAGGatgatgatg                                                                                  <  5:313307/149‑1 (MQ=255)                                                                  ggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGccc                                                                                                                                        >  4:17361/1‑94 (MQ=255)                                                                  ggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGccc                                                                                                                                        <  3:17361/94‑1 (MQ=255)                                                                       cgcgCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGTCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccac                                                                             >  1:56753/1‑148 (MQ=255)                                                                              gagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCAt                                                                       >  8:282533/1‑147 (MQ=255)                                                                                 tACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                  >  8:182932/1‑85 (MQ=255)                                                                                 tACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                  <  7:182932/85‑1 (MQ=255)                                                                                        gcgcCGCCAGGCGCTGAATGACCACATGTCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGt                                                            <  2:56753/148‑1 (MQ=255)                                                                                          gccgccCGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                              <  4:48985/144‑1 (MQ=255)                                                                                          gccgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                              >  3:48985/1‑144 (MQ=255)                                                                                          gccgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgct                                                         <  3:217819/149‑1 (MQ=255)                                                                                          gccgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgct                                                         >  6:239977/1‑149 (MQ=255)                                                                                           ccgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgctg                                                        >  5:87972/1‑149 (MQ=255)                                                                                              ccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                     >  5:277032/1‑149 (MQ=255)                                                                                               cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                     >  8:131768/1‑148 (MQ=255)                                                                                               cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtg                                                    <  7:112694/149‑1 (MQ=255)                                                                                                 ggCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                  <  3:75295/149‑1 (MQ=255)                                                                                                     cTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                              <  1:22374/149‑1 (MQ=255)                                                                                                     cTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAAt                                              <  1:28185/149‑1 (MQ=255)                                                                                                        aaTGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGTCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATgg                                            <  1:125636/148‑1 (MQ=255)                                                                                                        aaTGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccacca                                                                           <  5:37483/117‑1 (MQ=255)                                                                                                        aaTGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccacca                                                                           >  6:37483/1‑117 (MQ=255)                                                                                                        aaTGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATgg                                            <  7:131768/148‑1 (MQ=255)                                                                                                          tGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGCTGCCCACCACCCTCAGCGACAGGTGCTGCTTGATGAATGGCAc                                         >  2:302416/1‑149 (MQ=255)                                                                                                                   gCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                  >  2:315523/1‑131 (MQ=255)                                                                                                                   gCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                                  <  1:315523/131‑1 (MQ=255)                                                                                                                   gCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaa                                >  3:204362/1‑149 (MQ=255)                                                                                                                    ccGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCg                                                                  >  8:341996/1‑114 (MQ=255)                                                                                                                    ccGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCg                                                                  <  7:341996/114‑1 (MQ=255)                                                                                                                      gttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaagaa                             <  1:189152/149‑1 (MQ=255)                                                                                                                       ttcgttcTTAACCTGTTCAGTGCGCCTTGGCCGTACGGCCATGGCGCCAACGCAGCAAGCCTATAATCTTCGGCACCAGCTACAGCATGATTATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGaagaag                            <  1:293869/149‑1 (MQ=255)                                                                                                                        tcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGaagaag                            <  4:204362/148‑1 (MQ=255)                                                                                                                         cgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGt                           <  7:199870/148‑1 (MQ=255)                                                                                                                          gttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAg                         <  4:334955/149‑1 (MQ=255)                                                                                                                           ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                        <  3:93067/149‑1 (MQ=255)                                                                                                                                tACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAg                   >  6:287637/1‑149 (MQ=255)                                                                                                                                   cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCCCCACCATCAGCGACAGGTGCGGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTg                >  7:162178/1‑149 (MQ=255)                                                                                                                                    tGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa               >  8:56032/1‑149 (MQ=255)                                                                                                                                          gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCCTCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgc         >  2:172739/1‑149 (MQ=255)                                                                                                                                            gAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccg       >  7:24009/1‑149 (MQ=255)                                                                                                                                               ccTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGATGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAATCTGACAAGgccgcc     >  8:112045/1‑148 (MQ=255)                                                                                                                                                cTTGGCCGTACGGCCCAGGCGCCCACGCAGAAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCCCCAACACGGAACGGTGGTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAGGGTGACCAGGCCGCGGa   >  7:280553/1‑149 (MQ=255)                                                                                                                                                 ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                     >  3:343545/1‑98 (MQ=255)                                                                                                                                                 ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                     <  4:343545/98‑1 (MQ=255)                                                                                                                                                 ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc  >  4:102189/1‑149 (MQ=255)                                                                                                                                                 ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc  >  4:263338/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTCGCTGATGCCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAC  >  NC_002947/626411‑626703

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