Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,299,087	G→C	6.8%	P339A (CCG→GCG)	recD ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,299,087	0	G	C	6.8%	184.4 / 3.6	59	P339A (CCG→GCG)	recD	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease
Reads supporting (aligned to +/- strand):  ref base G (33/22);  new base C (2/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCCACTGACCTGCTGAAGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGTGCT  >  NC_002947/5298941‑5299234                                                                                                                                                   |                                                                                                                                                     gCCACTGACCTGCTGAAGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCggcggc                                                                                                                                                     >  5:347301/1‑148 (MQ=255)   cACTGACCTGCTGAAGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAg                                                                                                                                                  >  5:166561/1‑149 (MQ=255)        aCCTGCTGAAGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATc                                                                                                                                             <  6:347301/149‑1 (MQ=255)                aaGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCATCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAg                                                                                                                                       <  3:303765/147‑1 (MQ=255)                aaGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAgg                                                                                                                                      <  1:4089/148‑1 (MQ=255)                 aGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGtt                                                                                                                                    >  1:29801/1‑149 (MQ=255)                   ccACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTgg                                                                                                                                  >  5:54664/1‑149 (MQ=255)                   ccACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTgg                                                                                                                                  >  1:204628/1‑149 (MQ=255)                     aCGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGcc                                                                                                                                <  2:204628/149‑1 (MQ=255)                     aCGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGCCAAGGCATCGAGCAGGTTGGc                                                                                                                                 >  1:105748/1‑148 (MQ=255)                     aCGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGCCAAGGCATCGAGCAGGTTGGc                                                                                                                                 >  7:62116/1‑148 (MQ=255)                         tGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCcatca                                                                                                                            >  5:119146/1‑149 (MQ=255)                            aGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCAttt                                                                                                                         >  7:146901/1‑149 (MQ=255)                            aGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCAttt                                                                                                                         >  3:164920/1‑149 (MQ=255)                             gTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGGACCGCGCCGGCCTCGACCGAGGCCAACTGGGCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTc                                                                                                                        <  1:308529/149‑1 (MQ=255)                                 ccGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAgg                                                                                                                    <  6:54664/149‑1 (MQ=255)                                      gAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGat                                                                                                               <  2:29801/149‑1 (MQ=255)                                      gAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGa                                                                                                                >  1:192704/1‑148 (MQ=255)                                       aGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGat                                                                                                               >  5:280978/1‑148 (MQ=255)                                        gTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGCCCGCGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGGCAGGATTTCCAGGTCGatca                                                                                                             >  2:291595/1‑149 (MQ=255)                                             aGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGCCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGAc                                                                                                        <  8:62116/149‑1 (MQ=255)                                             aGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGCCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGAc                                                                                                        <  2:105748/149‑1 (MQ=255)                                                    ctcAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGt                                                                                                                   >  2:178234/1‑131 (MQ=255)                                                    ctcAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑Ttcgt                                                                                                >  2:291267/1‑149 (MQ=255)                                                    ctcAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑Ttcgt                                                                                                >  3:40732/1‑149 (MQ=255)                                                           tCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGaccac                                                                                         >  3:223834/1‑149 (MQ=255)                                                             ccGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑GTCGTCGaccac                                                                                         >  1:124024/1‑147 (MQ=255)                                                                   caGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACT‑GTCCTTGTCACCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATc                                                                                >  6:257982/1‑149 (MQ=255)                                                                           ccAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCCTCGAGCAGGTTGGCCATCATTTCCAGCTCGATCATCGACGC‑TTCGTCGACCACCAGCCCATCcaccgg                                                                          >  5:239787/1‑148 (MQ=255)                                                                           ccAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCcaacgg                                                                          >  3:158679/1‑148 (MQ=255)                                                                              gCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCcaacggcaac                                                                      >  2:17119/1‑149 (MQ=255)                                                                               cACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTAGTCTCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGCTTTGCCATCATTTCCAGCTCGATAATCGAAGC‑TTCGTCGACCACCAGCACATCcaacggcaacg                                                                     >  4:202380/1‑149 (MQ=255)                                                                                aCCGCGCCGGCCTCGACCGAGGCCAACT‑GTCCTTGTCACCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACggg                                                                   <  5:257982/149‑1 (MQ=255)                                                                                 ccGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACggg                                                                   <  6:166561/149‑1 (MQ=255)                                                                                  cgcgCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGt                                                                  <  8:146901/149‑1 (MQ=255)                                                                                    cgcCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTg                                                                <  1:17119/149‑1 (MQ=255)                                                                                     gttggCCTCGGTCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGc                                                               >  7:232655/4‑149 (MQ=255)                                                                                     gcCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGc                                                               <  6:280978/149‑1 (MQ=255)                                                                                        ggCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAg                                                            >  2:51149/1‑149 (MQ=255)                                                                                           cTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGCCGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑GTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCat                                                         >  6:190807/1‑149 (MQ=255)                                                                                           cTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGc                                                           >  4:7985/1‑147 (MQ=255)                                                                                           cTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGc                                                           <  3:7985/147‑1 (MQ=255)                                                                                           cTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCat                                                         >  8:1384/1‑149 (MQ=255)                                                                                           cTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCat                                                         >  6:364262/1‑149 (MQ=255)                                                                                                 cGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGc                                                   >  4:250054/1‑149 (MQ=255)                                                                                                     gCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTCCCAGCATGATGGCGa                                                 >  4:265059/1‑147 (MQ=255)                                                                                                             gTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCcgcg                                       >  6:368788/1‑149 (MQ=255)                                                                                                             gTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCcgc                                        <  4:164920/148‑1 (MQ=255)                                                                                                                       gccAGCAGCACCAGAGGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATGATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGggcg                             <  3:202380/148‑1 (MQ=255)                                                                                                                       cccAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGggcg                             >  3:103099/1‑149 (MQ=255)                                                                                                                       cccAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCAACGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCcgc                                        <  4:54434/138‑1 (MQ=255)                                                                                                                       cccAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCAACGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCcgc                                        >  3:54434/1‑138 (MQ=255)                                                                                                                               caCCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCACGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGa                     <  1:291267/149‑1 (MQ=255)                                                                                                                                aCCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGa                     <  3:244637/148‑1 (MQ=255)                                                                                                                                aCCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCCAg                    <  6:119146/149‑1 (MQ=255)                                                                                                                                   aGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAg                 <  5:364262/149‑1 (MQ=255)                                                                                                                                     aCGCGCCCGCGGCGGCAAGGCATCGCGCAGCGTGGCCATCATTTCCAGGTCGATCATCGACACGTTCGTCGACCACCAGCACATCCAGCGGCACCGGGTTGCCAGCATGATGGCGAAAGAGCCGCGCGCCCGGGCGGCTGCCGAGCAgg                >  1:103499/1‑149 (MQ=255)                                                                                                                                      cgcgCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCg              <  4:40732/149‑1 (MQ=255)                                                                                                                                         gcCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATg           >  1:198670/1‑149 (MQ=255)                                                                                                                                           ccGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCa         <  5:368788/149‑1 (MQ=255)                                                                                                                                               ggcggcAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGt     >  2:87125/1‑149 (MQ=255)                                                                                                                                                  ggcAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGTGCt  <  1:51149/149‑1 (MQ=255)                                                                                                                                                  ggcAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGTGCt  <  1:87125/149‑1 (MQ=255)                                                                                                                                                  ggcAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGTGCt  <  8:232655/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCACTGACCTGCTGAAGCCACGCTGCAGTCGCCGGCGAGTAAAAGCCCTCCTCAGCATCCCGGCACAGGTCACCCAGCACCGCGCCGGCCTCGACCGAGGCCAACTGGTCCTTGTCGCCCAGCAGCACCAGACGCGCCCGCGGCGGCAAGGCATCGAGCAGGTTGGCCATCATTTCCAGGTCGATCATCGACGC‑TTCGTCGACCACCAGCACATCCAACGGCAACGGGTTGCCAGCATGATGGCGAAAGTGCCGCGAGCCCGGGCGGCTGCCGAGCAGGCGATGCACGGTGCT  >  NC_002947/5298941‑5299234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑