Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,594,484	A→C	7.3%	R313R (CGT→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,484	0	A	C	7.3%	208.0 / 7.4	69	R313R (CGT→CGG)	kgtP	alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base C (2/3);  total (33/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 = 8.65e-01

CGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGT  >  NC_002947/1594336‑1594626                                                                                                                                                      |                                                                                                                                               cGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGa                                                                                                                                                <  8:39996/149‑1 (MQ=255)     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:263935/149‑1 (MQ=255)                  taGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               <  5:268379/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                                                                                                                               <  6:113914/146‑1 (MQ=255)                     ccGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTc                                                                                                                               >  5:113914/1‑146 (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                                                                                           >  4:241631/1‑147 (MQ=255)                                                        gaaCGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATgaa                                                                                           >  6:3377/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                                                                                 <  4:237030/149‑1 (MQ=255)                                                                cccACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGc                                                                                 >  7:169019/1‑149 (MQ=255)                                                                cccACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTACAGGCACATGAACAGGAACAGc                                                                                 >  8:279737/1‑149 (MQ=255)                                                                  caccaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGt                                                                               <  2:231080/149‑1 (MQ=255)                                                                   accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                     >  2:165844/1‑126 (MQ=255)                                                                   accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCa                                                                                                     <  1:165844/126‑1 (MQ=255)                                                                   accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTg                                                                              >  2:152361/1‑149 (MQ=255)                                                                   accaGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTg                                                                              <  3:241631/149‑1 (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                                                                                                  <  2:129375/106‑1 (MQ=255)                                                                                          gTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGc                                                                                                  >  1:129375/1‑106 (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                                    <  6:181598/134‑1 (MQ=255)                                                                                                                             gCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACg                                    >  5:181598/1‑134 (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)                                                                                                                                                      |                                                                                                                                               CGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCACCCCACCAGGTCTGGATGGTGTGCAGGGTGCTGA‑GGATGGGGACGGTGAACAGGGTGCCGAGCACACCGAAGGCGATCAGGATCGGACGGCGGCCGATCTTGTCGGACAGACCACCGATCATCGGCTGCAGGCACATGAACAGGAACAGCGTGGCGGCCGAGATGGTGGTCGAGTCACTGATGCTCATGCCCACGGTGTTGACCAGGTATTTCTGCATGTAGGTGGTGT  >  NC_002947/1594336‑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: ‑