Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,069,292	Δ1 bp	6.5%	coding (316/897 nt)	PP_5323 ←	M23/M37 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,069,292	0	G	.	6.5%	243.6 / 10.0	77	coding (316/897 nt)	PP_5323	M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base . (3/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

ACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTT  >  NC_002947/6069151‑6069438                                                                                                                                               |                                                                                                                                                     aCGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCGCCGGCTGCTGCCGCC‑GCGCAcc                                                                                                                                               <  3:94235/149‑2 (MQ=255)     gAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATAc                                                                                                                                           >  6:103861/1‑149 (MQ=255)         gCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCgg                                                                                                                                       >  2:389399/1‑149 (MQ=255)          cATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGGCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGGGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCggg                                                                                                                                      <  1:303894/149‑1 (MQ=255)                       cTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCtt                                                                                                                         >  3:122639/1‑149 (MQ=255)                       cTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCtt                                                                                                                         <  5:103861/149‑1 (MQ=255)                            gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                                    >  1:30787/1‑149 (MQ=255)                            gCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCg                                                                                                                    >  7:19526/1‑149 (MQ=255)                                tttGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAACAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTtgt                                                                                                                       >  3:394874/1‑142 (MQ=255)                                tttGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAACAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTtgt                                                                                                                       <  4:394874/142‑1 (MQ=255)                                   gCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGAc                                                                                                               >  8:61516/1‑147 (MQ=255)                                    cccGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCACCGCCTTGTGCGTTGACCAg                                                                                                            <  4:298332/149‑1 (MQ=255)                                     ccGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAgg                                                                                                           <  5:252742/149‑1 (MQ=255)                                        ggTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCGCCTGAGCACAATGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGACGCCTTGTGGGTTGACAAGGTTg                                                                                                        >  4:319614/1‑149 (MQ=255)                                        cgTCACTCAGGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTg                                                                                                        <  5:101969/148‑1 (MQ=255)                                            aCCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCa                                                                                                    <  2:353379/149‑1 (MQ=255)                                             cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGc                                                                                                                                    <  4:292888/116‑1 (MQ=255)                                             cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGc                                                                                                                                    >  3:292888/1‑116 (MQ=255)                                             cccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAg                                                                                                   >  4:50786/1‑149 (MQ=255)                                              ccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                                  >  7:249499/1‑149 (MQ=255)                                              ccATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                                  >  8:115229/1‑149 (MQ=255)                                               cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTg                                                                                                                 <  2:142591/133‑1 (MQ=255)                                               cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTg                                                                                                                 >  1:142591/1‑133 (MQ=255)                                               cATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGcc                                                                                                 <  3:50786/149‑1 (MQ=255)                                                agggCCTGCAGCAGGGTCGAACGGTGGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGTATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGc                                                                                                  <  4:345115/145‑1 (MQ=255)                                                aTGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCg                                                                                                >  3:95442/1‑149 (MQ=255)                                                   gCCTGCAGGAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCaa                                                                                             <  1:352739/149‑1 (MQ=255)                                                    ccTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCaa                                                                                             >  4:278573/1‑148 (MQ=255)                                                        cagcagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAActcta                                                                                        <  4:341575/149‑1 (MQ=255)                                                        cagcagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAActct                                                                                         <  3:2934/148‑1 (MQ=255)                                                           cagGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCTTTTTGTCCAGCCGCAActctacct                                                                                     <  7:61516/149‑1 (MQ=255)                                                                    aCGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                                          >  4:91658/1‑149 (MQ=255)                                                                     cGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                                          >  8:46674/1‑148 (MQ=255)                                                                      ggTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTTCGTGGTGCG‑CGGCGGAAGCAGCCGGCGCAAGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCAAGc                                                                                                  >  7:223526/1‑125 (MQ=255)                                                                      ggTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGACATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGcaca                                                                        <  1:389399/149‑1 (MQ=255)                                                                          gTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCa                                                                    <  4:122639/149‑1 (MQ=255)                                                                                     gcCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGt                                                         >  1:207790/1‑149 (MQ=255)                                                                                     gcCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGt                                                         <  3:278573/149‑1 (MQ=255)                                                                                      cccGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTCGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTg                                                        >  2:122257/1‑149 (MQ=255)                                                                                             cgccCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                                 <  2:30787/149‑1 (MQ=255)                                                                                             cgccCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                                 <  2:207790/149‑1 (MQ=255)                                                                                               cccccTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGc                                               >  5:2403/1‑149 (MQ=255)                                                                                                  ccTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCa                                            >  7:351396/1‑149 (MQ=255)                                                                                                     gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                                 <  3:394815/141‑1 (MQ=255)                                                                                                     gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGaaaa                                                 >  4:394815/1‑141 (MQ=255)                                                                                                     gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGc                                         >  4:43514/1‑149 (MQ=255)                                                                                                     gAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTACACCAGGTGATCCCCCCGCAACTCTACCTC‑T‑ACCCGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGc                                         >  1:132255/1‑149 (MQ=255)                                                                                                       gCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCa                                       >  5:300022/1‑149 (MQ=255)                                                                                                        cacaCTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCg                                                     >  4:266953/1‑134 (MQ=255)                                                                                                        cacaCTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCg                                                     <  3:266953/134‑1 (MQ=255)                                                                                                         acacTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcc                                     >  2:84313/1‑149 (MQ=255)                                                                                                          cacTGAGCATCTGCGTGGTTGG‑GGGGGGCAGCAGGCGGCGGACGTTACGCGGGGCATGGCCGGCTTGTTGGTTGTCCAGGTTTTCCATCCGCAACTCTAACTC‑TCACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcc                                     <  1:84313/149‑1 (MQ=255)                                                                                                          cacTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGAAGCTCGTTTCGTGCCGAAAAGCTCACGCCAcag                                    >  6:270878/1‑147 (MQ=255)                                                                                                           acTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGt                                   >  4:336235/1‑149 (MQ=255)                                                                                                               aGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCg                                 >  5:111607/1‑147 (MQ=255)                                                                                                               aGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGaa                               <  7:115229/149‑1 (MQ=255)                                                                                                                gCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACTCCACCGTCGAAg                              <  6:2403/149‑1 (MQ=255)                                                                                                                     tGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAc                                      <  1:364946/136‑1 (MQ=255)                                                                                                                     tGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCAc                                      >  2:364946/1‑136 (MQ=255)                                                                                                                      gtatcgtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAg                              >  8:334783/6‑143 (MQ=255)                                                                                                                      gtatcgtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAg                              <  7:334783/138‑1 (MQ=255)                                                                                                                      gCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        >  2:382976/1‑149 (MQ=255)                                                                                                                      gCGTGGTGCG‑CGGCGGCAACAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        >  1:392607/1‑149 (MQ=255)                                                                                                                       cGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        <  4:95442/148‑1 (MQ=255)                                                                                                                       cGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        <  3:91658/148‑1 (MQ=255)                                                                                                                        gtggtgCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACACCACCGTCGAAGGTTTTcgc                      >  1:215042/1‑149 (MQ=255)                                                                                                                           gtgCGCCGGCTGCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGtt                   >  6:28729/1‑149 (MQ=255)                                                                                                                            tgcgccggctgCTGCCGCC‑GCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTc                  >  1:394801/6‑149 (MQ=255)                                                                                                                            tgCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGGCGTCCAGCCGCAACTCTACCTCAC‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGtt                   <  1:122257/149‑1 (MQ=255)                                                                                                                                        agcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                                          <  4:159160/82‑1 (MQ=255)                                                                                                                                        agcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGca                                                                          >  3:159160/1‑82 (MQ=255)                                                                                                                                         gcagcCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTgg      <  6:300022/149‑1 (MQ=255)                                                                                                                                           tggcGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        >  4:87258/4‑129 (MQ=255)                                                                                                                                           tggcGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTc                        <  3:87258/126‑1 (MQ=255)                                                                                                                                             ccGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTg       >  7:300691/1‑144 (MQ=255)                                                                                                                                             ccGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTg       <  8:300691/144‑1 (MQ=255)                                                                                                                                             ccGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCtt  >  5:307557/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     ACGGGAACGCATAGCGGAAGCCCTGGGGCCGTTTGCCCGGGTCACCCATGGCCTGCAGCAGGGTCGAACGGTAGTTCAGCAGCAGCCCGGGCCGCCCCCTGAGCACACTGAGCATCTGCGTGGTGCG‑CGGCGGCAGCAGCCGGCGCACGATACGCGGGGCATCGCCGCCTTGTGCGTTGACCAGGTTGTCCAGCCGCAACTCTACCTC‑T‑ACCGGCACATGCAGCTCGTTTCGTGCCGAAAAGCTCACGCCACCGTCGAAGGTTTTCGCGTTCAGCAGTACCTGGCCTT  >  NC_002947/6069151‑6069438

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