Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,774,812	G→A	8.9%	P234S (CCA→TCA) ‡	PP_5521 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,774,812	0	G	A	8.9%	220.9 / 8.3	78	P234S (CCA→TCA) ‡	PP_5521	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (4/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

TCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCC  >  NC_002947/2774664‑2774960                                                                                                                                                          |                                                                                                                                                     tCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTgg                                                                                                                                                      <  8:171904/149‑1 (MQ=255) tCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTgg                                                                                                                                                      <  1:433092/149‑1 (MQ=255)  ccAGAACAAGATGGCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGgc                                                                                                                                                     <  7:333906/149‑1 (MQ=255)   cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGgca                                                                                                                                                    <  4:172808/149‑1 (MQ=255)   cAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGgca                                                                                                                                                    <  2:28565/149‑1 (MQ=255)   cAGAAAAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGgca                                                                                                                                                    <  6:300007/149‑1 (MQ=255)     gAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGgcagc                                                                                                                                                  <  2:187457/149‑1 (MQ=255)        cAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGTCCGGTGCTCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGcagcag                                                                                                                                                <  1:102078/149‑1 (MQ=255)        cAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAgg                                                                                                                                               <  8:410686/149‑1 (MQ=255)        cAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGCACATCCGCCGACACCGCACCATGCTCCATGCTCAGCACCAGGCTGATGTCGTCAGCGT‑CGGTGCCCTCGGCCAGCTCTTAGCGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAgg                                                                                                                                               >  8:332320/1‑149 (MQ=255)         aaGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGc                                                                                                                                              <  4:268771/149‑1 (MQ=255)           gATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGc                                                                                                                                              >  1:393460/1‑147 (MQ=255)           gATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGCCACAGCAGGcac                                                                                                                                            <  5:393611/149‑3 (MQ=255)           gATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGCGC‑‑‑CTGCTGCCACAGCAGGcac                                                                                                                                            <  3:157238/149‑3 (MQ=255)            aTGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCTGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAg                                                                                                                                           <  1:196425/149‑1 (MQ=255)             tGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGt                                                                                                                                          <  6:36819/149‑1 (MQ=255)              gTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTc                                                                                                                                         >  6:202374/1‑149 (MQ=255)                     gCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCtt                                                                                                                                  <  1:178751/149‑1 (MQ=255)                       cTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGa                                                                                                                                >  8:147805/1‑149 (MQ=255)                          gCTCACCCACGACGATCAGCCGTCCATCCGCCGACACCGCACAATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGt                                                                                                                              <  8:107958/148‑1 (MQ=255)                           cTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTg                                                                                                                             <  3:137882/148‑1 (MQ=255)                            tCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTggg                                                                                                                           <  2:381405/149‑2 (MQ=255)                            tCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGc                                                                                                                           >  3:358078/1‑149 (MQ=255)                             cACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGc                                                                                                                           >  8:321380/1‑148 (MQ=255)                              aCCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGcc                                                                                                                          <  4:210356/148‑1 (MQ=255)                               cccACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTc                                                                                                                                         <  3:171029/132‑1 (MQ=255)                               cccACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTc                                                                                                                                         >  4:171029/1‑132 (MQ=255)                                ccACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGcc                                                                                                                          >  6:100333/1‑146 (MQ=255)                                ccACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGcc                                                                                                                          <  5:100333/146‑1 (MQ=255)                                  acgacgATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTAc                                                                                                                     >  2:146797/1‑149 (MQ=255)                                    gacgaTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACg                                                                                                                    >  7:286751/1‑148 (MQ=255)                                     acgaTCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGaa                                                                                                                  <  1:173481/149‑1 (MQ=255)                                               ggCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCg                                                                                                        <  6:15946/149‑1 (MQ=255)                                                gCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGt                                                                                                                              <  5:77222/126‑1 (MQ=255)                                                gCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGt                                                                                                                              >  6:77222/1‑126 (MQ=255)                                                      cgccgACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCa                                                                                                 <  1:3364/149‑1 (MQ=255)                                                         cgACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCAc                                                                                                >  4:425568/1‑147 (MQ=255)                                                             acCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcagcag                                                                                          >  6:216755/1‑149 (MQ=255)                                                             acCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACcagcag                                                                                          >  3:329220/1‑149 (MQ=255)                                                               cGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGGAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTa                                                                                        <  1:124437/149‑1 (MQ=255)                                                                    catgctccatgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGc                                                                                   >  3:186485/1‑149 (MQ=255)                                                                      tgctccatgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTg                                                                                 <  5:202374/149‑1 (MQ=255)                                                                       gctccatgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTgg                                                                                >  4:9851/1‑149 (MQ=255)                                                                        ctccatgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGc                                                                               >  6:126804/1‑149 (MQ=255)                                                                          ccatgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCa                                                                             >  1:230192/1‑149 (MQ=255)                                                                           catgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCAt                                                                            >  3:406153/1‑149 (MQ=255)                                                                             tgctcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCg                                                                          <  5:170947/149‑1 (MQ=255)                                                                                tcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGa                                                                                                                   <  6:91961/105‑1 (MQ=255)                                                                                tcAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGa                                                                                                                   >  5:91961/1‑105 (MQ=255)                                                                                  aGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCgggg                                                                       >  6:14293/1‑147 (MQ=255)                                                                                    ccccAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAgg                                                                   >  2:183095/1‑149 (MQ=255)                                                                                     cccAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAggg                                                                  <  4:358078/149‑1 (MQ=255)                                                                                      ccAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAggg                                                                  >  1:158923/1‑148 (MQ=255)                                                                                          gCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGAt                                                               >  8:299948/1‑147 (MQ=255)                                                                                          gCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCa                                                             <  7:199555/149‑1 (MQ=255)                                                                                              aTGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAgcg                                                          <  8:286751/148‑1 (MQ=255)                                                                                               tttggttcGGGT‑GGGTGCCCTCGGGCATCTCTTTGAGCACGGCGCCCACACTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCa                                                             <  1:146797/141‑1 (MQ=255)                                                                                                  cgtcCGGGT‑CGGTGCCCTCGGACAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGcggtc                                                     >  8:225961/1‑149 (MQ=255)                                                                                                   gtcCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGcggtcg                                                    >  1:178432/1‑149 (MQ=255)                                                                                                   gtcCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGcggtcg                                                    <  3:9851/149‑1 (MQ=255)                                                                                                      cGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGcggtcggt                                                  <  7:321380/148‑1 (MQ=255)                                                                                                      cGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGgtg                                                 <  7:147805/149‑1 (MQ=255)                                                                                                             ggTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGggc                                            <  3:425568/148‑1 (MQ=255)                                                                                                                 cccTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAg                                       >  4:131215/1‑149 (MQ=255)                                                                                                                  ccTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAgg                                      <  2:393460/149‑1 (MQ=255)                                                                                                                         cAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACg                               >  6:409610/1‑149 (MQ=255)                                                                                                                           gCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGc                             <  1:183095/149‑1 (MQ=255)                                                                                                                               ttCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCg                          >  3:20361/1‑148 (MQ=255)                                                                                                                                  gAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAg                                                                    >  6:219242/1‑103 (MQ=255)                                                                                                                                  gAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAg                                                                    <  5:219242/103‑1 (MQ=255)                                                                                                                                  gAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAAc                       <  4:329220/148‑1 (MQ=255)                                                                                                                                    gCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCGCCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGGCGGTGTGGGGGGCAGGGCGAACGGCTCGCACggg                    >  8:365111/1‑149 (MQ=255)                                                                                                                                     cAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACg                               >  6:82908/1‑137 (MQ=255)                                                                                                                                     cAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACg                               <  5:82908/137‑1 (MQ=255)                                                                                                                                           tcgC‑‑‑CTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGggc                                            <  4:237995/117‑1 (MQ=255)                                                                                                                                           tcgC‑‑‑CTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGggc                                            >  3:237995/3‑119 (MQ=255)                                                                                                                                           tcgC‑‑‑CTGCTGCCACAGCAGGCGAGTAGCACCTTGAGTGGCGAGTACGATGATACCTTCGGCGCTCACCAGCACTACTGGCAGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGACGAACGGCTCGAACTGGTAACCt              >  5:226860/3‑149 (MQ=255)                                                                                                                                            cgC‑‑‑CTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACgg                     >  1:173404/2‑141 (MQ=255)                                                                                                                                             gC‑‑‑CTGCTGCCACAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTgg            >  7:276965/1‑149 (MQ=255)                                                                                                                                              c‑‑‑cTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTgg            <  4:406153/148‑1 (MQ=255)                                                                                                                                                  ctgctgTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggcg          <  2:178432/149‑1 (MQ=255)                                                                                                                                                  ctgctgTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggcg          >  2:214630/1‑149 (MQ=255)                                                                                                                                                     ctgTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggcg          <  7:360370/146‑1 (MQ=255)                                                                                                                                                     ctgTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTggcg          >  8:360370/1‑146 (MQ=255)                                                                                                                                                     ctgTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCa       >  4:170749/1‑149 (MQ=255)                                                                                                                                                         ggCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAAccc   >  7:264769/1‑149 (MQ=255)                                                                                                                                                          gcagcaGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAAcccc  <  7:225961/149‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                     TCCAGAACAAGATGTCGGCTGCCTTGCTCACCCACGACGATCAGCCGGCCATCCGCCGACACCGCACCATGCTCCATGCTCAGCCCCAGGCTGATGTCGTCCGGGT‑CGGTGCCCTCGGCCAGCTCTTCGAGCAC‑GCGTGC‑‑‑CTGCTGTGGCAGCAGGCGAGTCGCACCTTGAGTGGCCAGTACGAAGATACCTTCGGCGCTCACCAGCAGTACTCGCTGGCCATCGGGGAAGGGGATCAGCGCGGTCGGTGTGGGCGGCAGGTCGAACGGCTCGAACGGGTAACCTGGCGGCAACCCC  >  NC_002947/2774664‑2774960

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