Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,594,487	G→C	7.5%	R312R (CGC→CGG)	kgtP ←	alpha‑ketoglutarate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,594,487	0	G	C	7.5%	205.0 / 6.4	67	R312R (CGC→CGG)	kgtP	alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (2/3);  total (31/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.42e-01

GTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGT  >  NC_002947/1594340‑1594626                                                                                                                                                     |                                                                                                                                            gTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGAcggc                                                                                                                                            >  8:165042/1‑149 (MQ=255)    gATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACggcggc                                                                                                                                         >  5:76621/1‑149 (MQ=255)     aTCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGcc                                                                                                                                        >  6:254311/1‑149 (MQ=255)          ggTATAGCCGCTTACGATGATCAGTGCGGCCATGAGCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAAGATCGGCCGCCGTCCGATcc                                                                                                                                   <  7:248083/149‑2 (MQ=255)             ataGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGt                                                                                                                                >  1:211493/1‑149 (MQ=255)              taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGt                                                                                                                                <  5:254311/148‑1 (MQ=255)              taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               <  5:268379/149‑1 (MQ=255)              taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               <  5:263935/149‑1 (MQ=255)              taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               <  1:234767/149‑1 (MQ=255)                gCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACTTGGTCTGGATGGTGTGCAGGGTGCGGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCg                                                                                                                              <  1:271839/148‑1 (MQ=255)                gCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCg                                                                                                                              >  8:74756/1‑148 (MQ=255)                 ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               >  5:113914/1‑146 (MQ=255)                 ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               <  6:113914/146‑1 (MQ=255)                 ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGa                                                                                                                            >  4:275034/1‑149 (MQ=255)                   gCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACa                                                                                                                          <  2:211493/149‑1 (MQ=255)                         gatgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGaccac                                                                                                                    >  7:211814/1‑149 (MQ=255)                           tgatCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCg                                                                                                                  <  3:275034/149‑1 (MQ=255)                              acAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGAAGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTAGt                                                                                                                                >  5:311510/2‑132 (MQ=255)                                    cgaccATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCgg                                                                                                          <  7:288828/145‑1 (MQ=255)                                        cATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATc                                                                                                                                    >  3:286364/1‑118 (MQ=255)                                        cATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATc                                                                                                                                    <  4:286364/118‑1 (MQ=255)                                               aGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACAt                                                                                              <  6:76621/149‑1 (MQ=255)                                               aGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACAt                                                                                              >  3:16295/1‑149 (MQ=255)                                                 gaagaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGAGGGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACAt                                                                                              >  1:67388/1‑148 (MQ=255)                                                   agaaCGCACCCCACCATGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaac                                                                                          <  3:175100/149‑1 (MQ=255)                                                    gaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaa                                                                                           >  6:3377/1‑147 (MQ=255)                                                    gaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaa                                                                                           >  4:241631/1‑147 (MQ=255)                                                     aaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaacag                                                                                        <  7:74756/149‑1 (MQ=255)                                                       cGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaacagga                                                                                      <  7:165042/149‑1 (MQ=255)                                                        gCACCCCAGCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaacaggaa                                                                                     <  2:78410/149‑1 (MQ=255)                                                            cccACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                 >  7:169019/1‑149 (MQ=255)                                                            cccACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                 <  4:237030/149‑1 (MQ=255)                                                            cccACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTACAGGCACATGAACAGGAACAGc                                                                                 >  8:279737/1‑149 (MQ=255)                                                              caccaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                               <  2:231080/149‑1 (MQ=255)                                                               accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                     <  1:165844/126‑1 (MQ=255)                                                               accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                     >  2:165844/1‑126 (MQ=255)                                                               accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTg                                                                              <  3:241631/149‑1 (MQ=255)                                                               accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTg                                                                              >  2:152361/1‑149 (MQ=255)                                                               accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTg                                                                              >  2:157221/1‑149 (MQ=255)                                                                         gATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGcc                                                                        >  3:212277/1‑145 (MQ=255)                                                                         gATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCga                                                                      <  1:233414/147‑1 (MQ=255)                                                                           tGGTGTGTAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAAGATCGGCCGCCGTCCGATCCTgat                                                                                                                               <  5:202156/88‑3 (MQ=37)                                                                           tGGTGTGTAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAAGATCGGCCGCCGTCCGATCCTgat                                                                                                                               >  6:202156/1‑86 (MQ=37)                                                                              tgtgCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGAtggtg                                                               <  5:3377/149‑1 (MQ=255)                                                                                  cAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTc                                                            >  2:123933/1‑148 (MQ=255)                                                                                   aGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                          <  1:157221/149‑1 (MQ=255)                                                                                   aGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGa                                                          <  1:152361/149‑1 (MQ=255)                                                                                     ggTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGt                                                        <  1:123933/149‑1 (MQ=255)                                                                                      gTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGc                                                                                                  >  1:129375/1‑106 (MQ=255)                                                                                      gTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGc                                                                                                  <  2:129375/106‑1 (MQ=255)                                                                                        gCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCAc                                                     <  4:16295/149‑1 (MQ=255)                                                                                           gA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGa                                                  >  4:279390/1‑149 (MQ=255)                                                                                                 tGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTc                                             <  2:67388/149‑1 (MQ=255)                                                                                                        ggTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTACAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCa                                      <  7:279737/149‑1 (MQ=255)                                                                                                         gTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCAc                                     >  3:269195/1‑149 (MQ=255)                                                                                                                         gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACg                                    >  5:181598/1‑134 (MQ=255)                                                                                                                         gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACg                                    <  6:181598/134‑1 (MQ=255)                                                                                                                               cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTg                >  1:8932/1‑148 (MQ=255)                                                                                                                               cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTg                <  2:8932/148‑1 (MQ=255)                                                                                                                               cGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGc               <  4:212277/149‑1 (MQ=255)                                                                                                                                   ggCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGt           <  8:169019/149‑1 (MQ=255)                                                                                                                                    gCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATg            >  1:25144/1‑147 (MQ=255)                                                                                                                                       aTCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAggt       >  2:1530/1‑149 (MQ=255)                                                                                                                                       aTCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACAACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGAGCAGGTATTTCTGCATGTAggt       >  3:81502/1‑149 (MQ=255)                                                                                                                                           ggATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGgtggtg   >  5:219066/1‑149 (MQ=255)                                                                                                                                            gATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGACTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTc                  <  3:270094/133‑1 (MQ=255)                                                                                                                                            gATCGGCCGCCGTCCGATCTTGTCGGACAGACCACCGATCATCGACTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTc                  >  4:270094/1‑133 (MQ=255)                                                                                                                                            gATCGGACGGCGGCGGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAAAAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGgtggtg   >  6:270697/1‑148 (MQ=255)                                                                                                                                            gATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTggtgt  <  3:279390/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            GTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGT  >  NC_002947/1594340‑1594626

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: ‑