Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,103,217	T→A	6.4%	L17Q (CTG→CAG)	PP_5353 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,103,217	0	T	A	6.4%	281.9 / 5.2	92	L17Q (CTG→CAG)	PP_5353	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (46/40);  new base A (2/4);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

ATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCTTCG  >  NC_002947/6103070‑6103365                                                                                                                                                    |                                                                                                                                                     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                                                                                      <  7:194164/149‑1 (MQ=255)                                                                cgcgcgTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccac                                                                                      <  8:205069/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                                               >  5:50864/1‑149 (MQ=255)                                                                                                       aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                               >  8:236717/1‑149 (MQ=255)                                                                                                          gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                                    <  5:306922/125‑1 (MQ=255)                                                                                                          gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                                    >  6:306922/1‑125 (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                            <  3:66201/149‑1 (MQ=255)                                                                                                                          aCTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTg                            <  1:124486/149‑1 (MQ=255)                                                                                                                              gcgcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCAcc                        >  6:189910/1‑149 (MQ=255)                                                                                                                                gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg                      >  6:141142/1‑149 (MQ=255)                                                                                                                                gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg                      <  5:334853/149‑1 (MQ=255)                                                                                                                                      cGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgc                  <  7:342259/147‑1 (MQ=255)                                                                                                                                       gTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgccgc               <  4:63105/149‑1 (MQ=255)                                                                                                                                            ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          >  3:368988/1‑149 (MQ=255)                                                                                                                                            ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          <  6:34991/149‑1 (MQ=255)                                                                                                                                            ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          >  5:34991/1‑149 (MQ=255)                                                                                                                                            ccTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          >  7:292244/1‑149 (MQ=255)                                                                                                                                            ccTGACGATGAGCCTGGACGCCATCGCCAAGCGTTTCGGCGGCGGCAAGATCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGTCCGTACCTGGCCGTGCACCGGCTGCCGACAGCg          <  6:401643/149‑1 (MQ=255)                                                                                                                                                aCGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGc      >  5:308303/1‑149 (MQ=255)                                                                                                                                                  gCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          >  4:148458/1‑143 (MQ=255)                                                                                                                                                  gCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCg          <  3:148458/143‑1 (MQ=255)                                                                                                                                                  gCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCtt    >  4:273217/1‑149 (MQ=255)                                                                                                                                                   cTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCTTc   <  6:93721/149‑1 (MQ=255)                                                                                                                                                    tGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCTTcg  >  2:20791/1‑149 (MQ=255)                                                                                                                                                    tGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCTTcg  >  6:225786/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGCCAGCGGTGCTTCG  >  NC_002947/6103070‑6103365

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