Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,848,343	C→G	13.9%	A36A (GCC→GCG)	PP_3397 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,848,343	0	C	G	13.9%	195.6 / 22.3	72	A36A (GCC→GCG)	PP_3397	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base G (4/6);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

CAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGCGCCCGCC  >  NC_002947/3848197‑3848490                                                                                                                                                   |                                                                                                                                                    cAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTcgccgc                                                                                                                                                   >  6:41686/1‑149 (MQ=255)     gtttgtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAACGTCTCAGCCTTCAAGGATGCAGCATGAACACTTGGTCGCGCCGTGCGCTGGTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               <  3:252847/145‑1 (MQ=255)     ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               <  8:5063/149‑1 (MQ=255)     ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               <  1:336000/149‑1 (MQ=255)     ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACATGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               <  3:247923/149‑1 (MQ=255)      gtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGt                                                                                                                                              >  1:234193/1‑149 (MQ=255)              gACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTc                                                                                                                                      >  2:192185/1‑149 (MQ=255)               aCTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTc                                                                                                                                      >  1:319815/1‑148 (MQ=255)                cTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGa                                                                                                                                    <  2:234193/149‑1 (MQ=255)                   cGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATc                                                                                                                                  <  7:279957/148‑1 (MQ=255)                   cGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTGTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATcc                                                                                                                                 >  4:285654/1‑149 (MQ=255)                       gACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCgggg                                                                                                                             <  3:162001/149‑1 (MQ=255)                           tggtgctggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCt                                                                                                                         <  2:319815/149‑1 (MQ=255)                                  ggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               >  3:90329/1‑120 (MQ=255)                                  ggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                               <  4:90329/120‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCtt                                                                                                                  <  2:211278/148‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                                 >  1:62036/1‑149 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                                 <  5:90262/149‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACATGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                                 >  1:270527/1‑149 (MQ=255)                                        aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                                             >  4:152420/1‑148 (MQ=255)                                         tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                                       >  5:174255/1‑137 (MQ=255)                                         tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                                       <  6:174255/137‑1 (MQ=255)                                          gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAg                                                                                                          >  1:248659/1‑149 (MQ=255)                                          gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAg                                                                                                          >  6:102705/1‑149 (MQ=255)                                           aaCAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAgg                                                                                                         >  1:318968/1‑149 (MQ=255)                                              aGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGCTCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                                      >  8:199004/1‑149 (MQ=255)                                              aGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                                      <  1:192185/149‑1 (MQ=255)                                                aCAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGc                                                                                                     >  5:67359/1‑148 (MQ=255)                                                          aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgtg                                                                                          <  4:299439/149‑1 (MQ=255)                                                          aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgtg                                                                                          <  5:142809/149‑1 (MQ=255)                                                           gCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTgg                                                                                         >  2:362600/1‑149 (MQ=255)                                                           gCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTgg                                                                                         <  2:318968/149‑1 (MQ=255)                                                                 aaGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAg                                                                                   <  2:248659/149‑1 (MQ=255)                                                                   ggaggCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGa                                                                                 <  5:41686/149‑1 (MQ=255)                                                                         agcaTGAACGCTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                                           >  4:336091/1‑149 (MQ=255)                                                                         agcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                                           >  1:171250/1‑149 (MQ=255)                                                                         agcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                                           >  5:19853/1‑149 (MQ=255)                                                                          gcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGCTGTGGTACAAGGACAGCCTc                                                                          <  6:67359/149‑1 (MQ=255)                                                                           caTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCg                                                                         <  2:62036/149‑1 (MQ=255)                                                                                  acTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATc                                                                  >  8:46321/1‑149 (MQ=255)                                                                                      ggTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATccc                                                                <  8:254738/147‑1 (MQ=255)                                                                                       gTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTgg                                                             <  3:152420/149‑1 (MQ=255)                                                                                        tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                                        >  3:336743/1‑137 (MQ=255)                                                                                        tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                                        <  4:336743/137‑1 (MQ=255)                                                                                        tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGa                                                            <  1:137679/149‑1 (MQ=255)                                                                                          gcgcCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGc                                                          >  3:169435/1‑149 (MQ=255)                                                                                           cgcCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCt                                                         >  3:351136/1‑149 (MQ=255)                                                                                                  cgcTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAg                                                  >  8:341517/1‑149 (MQ=255)                                                                                                         ggCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTg                                           >  4:257915/1‑149 (MQ=255)                                                                                                           ccgccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGc                                          >  8:223650/1‑148 (MQ=255)                                                                                                             gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                                   <  3:68786/137‑1 (MQ=255)                                                                                                             gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                                       >  4:4849/1‑149 (MQ=255)                                                                                                             gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                                       >  3:308331/1‑149 (MQ=255)                                                                                                             gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACATTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                                   >  4:68786/1‑137 (MQ=255)                                                                                                               cgGTTTCTCTCTTACCGTTACCAGCGCCGGGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccacca                                     <  2:270527/149‑1 (MQ=255)                                                                                                               cgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccacca                                     <  1:362600/149‑1 (MQ=255)                                                                                                                      tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                                   <  4:36738/128‑1 (MQ=255)                                                                                                                      tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                                   >  3:36738/1‑128 (MQ=255)                                                                                                                        tctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTg                            >  7:345519/1‑149 (MQ=255)                                                                                                                         ctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgc                           >  8:286993/1‑149 (MQ=255)                                                                                                                          ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                            >  7:207261/1‑131 (MQ=255)                                                                                                                          ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                            <  8:207261/131‑1 (MQ=255)                                                                                                                             ccGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgcgc                         >  7:29122/1‑147 (MQ=255)                                                                                                                              atcagacAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAg                      >  1:137333/7‑149 (MQ=255)                                                                                                                               gTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAgg                     <  8:345519/149‑1 (MQ=255)                                                                                                                                   acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                                       <  6:22684/126‑1 (MQ=255)                                                                                                                                   acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                                       >  5:22684/2‑127 (MQ=255)                                                                                                                                   acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCa                 >  2:120691/2‑149 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCaa                                                                    <  4:77112/97‑1 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCaa                                                                    >  3:77112/1‑97 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGTCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc                >  7:276039/1‑149 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGCGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc                <  3:285654/149‑1 (MQ=255)                                                                                                                                             gTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGgcgc       <  7:286993/149‑1 (MQ=255)                                                                                                                                                cgccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCCCCGGGTGcgcct     <  3:4849/148‑2 (MQ=255)                                                                                                                                                cgccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGcgccc     <  3:257915/148‑1 (MQ=255)                                                                                                                                                 gccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGAGGTACAAGGACAGCCTCGACCAAGCCCTGGAGCTGTGCCAGTCGCGTGCCCCCAGCCGCCGTGCGCCCGGCGCCCCCCGGGcgcccgc   >  4:277987/1‑149 (MQ=255)                                                                                                                                                  ccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCGGGCGCCCCGGGTGcgcctggc  <  3:336091/149‑5 (MQ=255)                                                                                                                                                  ccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGcgcccgcc  <  2:171250/149‑1 (MQ=255)                                                                                                                                                  ccgcGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGcgcccgcc  <  6:19853/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGCGCCCGCC  >  NC_002947/3848197‑3848490

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