Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,103,211	T→C	8.9%	L15P (CTG→CCG)	PP_5353 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,103,211	0	T	C	8.9%	254.8 / 11.0	90	L15P (CTG→CCG)	PP_5353	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (45/37);  new base C (3/5);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCG  >  NC_002947/6103067‑6103357                                                                                                                                                 |                                                                                                                                                   aaCATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACg                                                                                                                                                >  2:302869/1‑149 (MQ=255)    aTTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCt                                                                                                                                              >  5:210568/1‑148 (MQ=255)    aTTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTg                                                                                                                                             <  1:137007/149‑1 (MQ=255)     ttCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCt                                                                                                                                              >  5:25504/1‑147 (MQ=255)      tCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAg                                                                                                                                           >  4:382907/1‑149 (MQ=255)        aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                          <  1:302869/148‑1 (MQ=255)        aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGcg                                                                                                                                         <  1:358822/149‑2 (MQ=255)         aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                          >  8:380887/1‑147 (MQ=255)         aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                          <  7:380887/147‑1 (MQ=255)         aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCAGGCCGACGCAGAGcgc                                                                                                                                        <  5:202318/149‑3 (MQ=255)         aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGcga                                                                                                                                        >  5:167655/1‑147 (MQ=255)         aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                          <  6:167655/147‑1 (MQ=255)          gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTg                                                                                                                                       >  8:197429/1‑149 (MQ=255)          gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                          >  1:131517/1‑146 (MQ=255)          gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                          <  2:131517/146‑1 (MQ=255)            aaaaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGa                                                                                                                                     >  8:220813/1‑149 (MQ=255)             aaaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGAc                                                                                                                                    <  3:382907/149‑1 (MQ=255)              aaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACg                                                                                                                                   >  8:194164/1‑149 (MQ=255)               aGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGc                                                                                                                                  <  4:359108/149‑1 (MQ=255)                    gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                             >  2:331750/1‑149 (MQ=255)                    gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                             >  4:388294/1‑149 (MQ=255)                    gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                             >  8:332051/1‑149 (MQ=255)                    gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCAGGACGCCAgcg                                                                                                                             >  3:248582/1‑149 (MQ=255)                     aaaCACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcgc                                                                                                                            >  3:376471/1‑149 (MQ=255)                     aaaCACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcgc                                                                                                                            >  3:241139/1‑149 (MQ=255)                       acacTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGcc                                                                                                                           >  6:285474/1‑148 (MQ=255)                           tATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGAGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCCGACGCTGa                                                                                                                                            <  8:284969/127‑1 (MQ=255)                           tATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGAGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCCGACGCTGa                                                                                                                                            >  7:284969/1‑127 (MQ=255)                             tGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGtt                                                                                                                    <  1:331750/149‑1 (MQ=255)                               tCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGttt                                                                                                                   <  4:241139/148‑1 (MQ=255)                                cGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                         >  7:119672/1‑141 (MQ=255)                                cGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                         <  8:119672/141‑1 (MQ=255)                                 gCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTcgg                                                                                                                >  3:63105/1‑149 (MQ=255)                                     gTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTcggcgg                                                                                                             <  5:285474/148‑1 (MQ=255)                                        aTGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                         <  3:142928/149‑1 (MQ=255)                                          gtgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                         >  8:342259/1‑147 (MQ=255)                                          gtgCCCGGGAACCTGTGCACCATGACGCGAGTGCAATCCACAACAGGGGATATCCAACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGCCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                         >  6:28280/1‑147 (MQ=255)                                           tgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                         >  5:331404/1‑146 (MQ=255)                                           tgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCaa                                                                                                       <  2:298648/148‑1 (MQ=255)                                             cccGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAaga                                                                                                     <  6:25504/148‑1 (MQ=255)                                               cGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCt                                                                                                  >  2:124486/1‑149 (MQ=255)                                                 ggAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCg                                                                                                <  4:376471/149‑1 (MQ=255)                                                   aaCCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCggg                                                                                              <  7:220813/149‑1 (MQ=255)                                                      cTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTc                                                                                           >  7:205069/1‑149 (MQ=255)                                                         tgCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTcc                                                                                          <  7:285452/147‑1 (MQ=255)                                                           cACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                      <  8:227134/149‑1 (MQ=255)                                                            aCCATGACGCGCGTGCAATCCACAACAGGGGATATCACACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGcc                                                                                     >  6:145132/1‑149 (MQ=255)                                                              cATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCTGCAAGAGCTCGGGCTCCGCGCCTa                                                                                   <  2:224378/149‑1 (MQ=255)                                                              cATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTa                                                                                   >  4:160227/1‑149 (MQ=255)                                                               aTGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                         <  8:270665/110‑1 (MQ=255)                                                               aTGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTAAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                         >  7:270665/1‑110 (MQ=255)                                                                   cgcgcgTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccac                                                                              <  8:205069/149‑1 (MQ=255)                                                                   cgcgcgTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccac                                                                              <  7:194164/149‑1 (MQ=255)                                                                          gCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGAccc                                                                       >  4:66201/1‑149 (MQ=255)                                                                            aaTCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGc                                                                     >  6:334853/1‑149 (MQ=255)                                                                               ccACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTTGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                      >  6:105607/1‑129 (MQ=255)                                                                               ccACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTTGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                      <  5:105607/129‑1 (MQ=255)                                                                                cacaACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCa                                                                   <  1:369244/147‑1 (MQ=255)                                                                                cacaACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCa                                                                   >  2:369244/1‑147 (MQ=255)                                                                                 acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccacca                                                                            >  3:314477/1‑137 (MQ=255)                                                                                 acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccacca                                                                            <  4:314477/137‑1 (MQ=255)                                                                                 acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGc                                                                >  6:170639/1‑149 (MQ=255)                                                                                   aacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGATTCGGCTCCGGCAAGAGCTCGGGCTCCGCGCATATCCACCAGACCCGCCAGGCTa                                                              >  5:401643/1‑149 (MQ=255)                                                                                        gggATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAgg                                                                 <  2:362031/141‑1 (MQ=255)                                                                                        gggATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAgg                                                                 >  1:362031/1‑141 (MQ=255)                                                                                          gATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAAc                                                       <  7:197429/149‑1 (MQ=255)                                                                                               cccACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGc                                                  >  8:15452/1‑149 (MQ=255)                                                                                                       tGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTTACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCtgctgc                                           <  3:160227/148‑1 (MQ=255)                                                                                                         cAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgc                                        >  5:93721/1‑149 (MQ=255)                                                                                                          aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGTGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                       >  4:255825/1‑149 (MQ=255)                                                                                                          aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                       >  8:236717/1‑149 (MQ=255)                                                                                                          aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                       >  5:50864/1‑149 (MQ=255)                                                                                                             gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                            >  6:306922/1‑125 (MQ=255)                                                                                                             gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                            <  5:306922/125‑1 (MQ=255)                                                                                                              ttttCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGa                                   <  6:331404/149‑1 (MQ=255)                                                                                                                         tCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGc                        >  6:412979/1‑149 (MQ=255)                                                                                                                           gCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGcc                       <  7:332051/148‑1 (MQ=255)                                                                                                                            cACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGt                     <  1:271026/149‑1 (MQ=255)                                                                                                                             aCTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTg                    <  1:124486/149‑1 (MQ=255)                                                                                                                             aCTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTg                    <  3:66201/149‑1 (MQ=255)                                                                                                                                 gcgcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCAcc                >  6:189910/1‑149 (MQ=255)                                                                                                                                   gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg              <  5:334853/149‑1 (MQ=255)                                                                                                                                   gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg              >  6:141142/1‑149 (MQ=255)                                                                                                                                         cGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgc          <  7:342259/147‑1 (MQ=255)                                                                                                                                          gTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgccgc       <  4:63105/149‑1 (MQ=255)                                                                                                                                               ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg  >  5:34991/1‑149 (MQ=255)                                                                                                                                               ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg  >  7:292244/1‑149 (MQ=255)                                                                                                                                               ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg  >  3:368988/1‑149 (MQ=255)                                                                                                                                               ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg  <  6:34991/149‑1 (MQ=255)                                                                                                                                               ccTGACGATGAGCCTGGACGCCATCGCCAAGCGTTTCGGCGGCGGCAAGATCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGTCCGTACCTGGCCGTGCACCGGCTGCCGACAGCg  <  6:401643/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   AACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCG  >  NC_002947/6103067‑6103357

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