Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,305,657	A→G	6.7%	L532S (TTG→TCG)	recC ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,305,657	0	A	G	6.7%	212.6 / 4.4	74	L532S (TTG→TCG)	recC	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (3/2);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCGCTGAC  >  NC_002947/5305518‑5305805                                                                                                                                             |                                                                                                                                                       gCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATAGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCCACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCCGGc                                                                                                                                               >  4:160671/1‑149 (MQ=255)    gATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCg                                                                                                                                            >  6:52125/1‑149 (MQ=255)            ccGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCGAGCCTGGGCTTGCCCGccgg                                                                                                                                    <  8:35034/149‑4 (MQ=255)             cgcgTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTgc‑g                                                                                                                                  >  1:270176/1‑149 (MQ=255)               cgTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGc‑gcg                                                                                                                                >  5:283289/1‑149 (MQ=255)                 tCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑gcggc                                                                                                                              <  7:80847/149‑1 (MQ=255)                    aGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTc                                                                                                                           <  3:205977/149‑1 (MQ=255)                            gATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGcccc                                                                                                                   <  6:283289/149‑1 (MQ=255)                                tCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACACCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACg                                                                                                               >  5:227908/1‑149 (MQ=255)                                  aTCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAACAACCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGCCACGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACg                                                                                                               >  1:98444/1‑148 (MQ=255)                                   tCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGAt                                                                                                            >  5:195906/1‑149 (MQ=255)                                   tCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACACCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGAt                                                                                                            <  6:227908/149‑1 (MQ=255)                                       aGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGc                                                                                                          >  6:185537/1‑147 (MQ=255)                                           ttCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAgcgc                                                                                                    >  1:218274/1‑149 (MQ=255)                                           ttCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAgcgc                                                                                                    >  1:213001/1‑149 (MQ=255)                                              aaTGCGATGACAGGCCTCACCCACCCCCACGTCATAACCCATC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCTATCCCGGCTTGCAAGCCCAGTCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCt                                                                                                 <  7:260938/149‑1 (MQ=255)                                              aaTGCCATCACAGGCCTCACCCATCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCt                                                                                                 >  5:35080/1‑149 (MQ=255)                                                 gCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGa                                                                                             >  7:169966/1‑149 (MQ=255)                                                 gCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCCCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGa                                                                                             >  1:309878/1‑149 (MQ=255)                                                   cATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATc                                                                                           >  8:52099/1‑149 (MQ=255)                                                    aTCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑Cg                                                                                              >  5:205364/1‑145 (MQ=255)                                                    aTCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑Cg                                                                                              <  6:205364/145‑1 (MQ=255)                                                        caGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCg                                                                                      <  5:185537/149‑1 (MQ=255)                                                                acccacccCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCACCg                                                                              >  7:116513/1‑149 (MQ=255)                                                                 cccacccCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGt                                                                             <  8:136509/149‑1 (MQ=255)                                                                 cccacccCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGt                                                                             >  8:29896/1‑149 (MQ=255)                                                                    acccCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCGCT‑CGATCCAGCGGTGCAGCGTggg                                                                          <  7:52099/149‑1 (MQ=255)                                                                      ccccACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCACGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCCGCGCCCT‑CGATCCAGCGGTCCAGCGgggg                                                                          >  8:264068/1‑147 (MQ=255)                                                                        ccACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAgg                                                                      >  2:352745/1‑149 (MQ=255)                                                                         cACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGt                                                                     <  6:195906/149‑1 (MQ=255)                                                                         cACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGt                                                                     >  6:255656/1‑149 (MQ=255)                                                                                  aCCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGCTCCAGCGGTGCAGCGTGGGCAGGTCGTTTTcgc                                                            >  2:230686/1‑149 (MQ=255)                                                                                  aCCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGt                                                                             >  4:61253/1‑132 (MQ=255)                                                                                  aCCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGt                                                                             <  3:61253/132‑1 (MQ=255)                                                                                   cccAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAg                                                                       <  4:295792/137‑1 (MQ=255)                                                                                   cccAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAg                                                                       >  3:295792/1‑137 (MQ=255)                                                                                   cccAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTcgc                                                            <  6:275792/148‑1 (MQ=255)                                                                                   cccAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTcgc                                                            >  5:275792/1‑148 (MQ=255)                                                                                    ccAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTcgcg                                                           >  1:177530/1‑148 (MQ=255)                                                                                              tGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAg                                                 <  1:352745/149‑1 (MQ=255)                                                                                               gcggcgCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGa                                                                                             <  4:96981/104‑1 (MQ=255)                                                                                               gcggcgCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGa                                                                                             >  3:96981/1‑104 (MQ=255)                                                                                                cggcgCAGGCCGAAACGCCAACTGGTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGc                                                <  8:67098/148‑1 (MQ=255)                                                                                                 ggcgCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCgg                                              <  5:52125/149‑1 (MQ=255)                                                                                                        gCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGa                                       <  2:218274/149‑1 (MQ=255)                                                                                                        gCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGa                                       <  2:213001/149‑1 (MQ=255)                                                                                                        gCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGa                                       >  4:56423/1‑149 (MQ=255)                                                                                                         ccGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCgg                                              >  5:345527/1‑141 (MQ=255)                                                                                                         ccGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCgg                                              <  6:345527/141‑1 (MQ=255)                                                                                                          cGAAACGCGAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTCCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACt                                     >  6:228292/1‑149 (MQ=255)                                                                                                          cGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACt                                     >  4:273759/1‑149 (MQ=255)                                                                                                               cGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCgt                                <  2:177530/149‑2 (MQ=255)                                                                                                               cGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCgg                                <  3:56423/149‑1 (MQ=255)                                                                                                               cGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCgg                                <  1:230686/149‑1 (MQ=255)                                                                                                               cGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCgg                                <  6:35080/149‑1 (MQ=255)                                                                                                                gCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGCGGCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCgg                                >  1:269028/1‑149 (MQ=255)                                                                                                                 ccAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCa                              >  4:164523/1‑149 (MQ=255)                                                                                                                   aaCTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACg                            <  5:228292/149‑1 (MQ=255)                                                                                                                     ctgttctgttCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTc                          <  2:270176/149‑1 (MQ=255)                                                                                                                        ttctgttcGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCcag                       >  1:266342/1‑149 (MQ=255)                                                                                                                        ttctgttcGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCcag                       >  5:35001/1‑149 (MQ=255)                                                                                                                          ctgttcGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCcag                       >  6:136330/1‑147 (MQ=255)                                                                                                                               cGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTc                 >  8:124496/1‑148 (MQ=255)                                                                                                                               cGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTc                 <  8:169966/148‑1 (MQ=255)                                                                                                                                gAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTc                 >  3:250284/1‑147 (MQ=255)                                                                                                                                 aGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAg              <  6:35001/149‑1 (MQ=255)                                                                                                                                  gcgcagcggGCGAGCCTGGGCTTGCCCGCCGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAg              >  7:227737/6‑148 (MQ=255)                                                                                                                                  gcgcagcggGCGAGCCTGGGCTTGCCCGCCGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAg              <  8:227737/143‑1 (MQ=255)                                                                                                                                  gCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGc             >  8:141150/1‑149 (MQ=255)                                                                                                                                   cccGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCa            >  8:177424/1‑149 (MQ=255)                                                                                                                                    ccGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCAc           <  7:124496/149‑1 (MQ=255)                                                                                                                                       gcggGCGAGCCTGGGCTTGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTc        >  4:84320/1‑149 (MQ=255)                                                                                                                                        cggGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCg       <  3:273759/149‑1 (MQ=255)                                                                                                                                         gggCGAGCCTGGGCTTGCCCGCCGG‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCGc      >  1:164330/1‑149 (MQ=255)                                                                                                                                             aaGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCTGCGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCGCTGa   >  1:277906/1‑149 (MQ=255)                                                                                                                                             aaGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCGCTGAc  <  2:98444/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                       GCCGATCAGTGCCGCGTCCAGGCCACCGATTTCATCGTAGGGTTCAATGCCATCACAGGCCTCACCCACCCCCACGGCATAACCCAAC‑AGCATGCGGCGCAGGCCGAAACGCCAACTGTTCTGTTCGAGCCCGGCGGGCAAGCCCAGGCTCGCCCGCTGC‑GCGGCGTCAAGGCCCCAACGGATGCCAGCGCCCT‑CGATCCAGCGGTGCAGCGTGGGCAGGTCGTTTTCGCGGATGCCGAAGCGGGCGCGGACTGCCGGCACGTCCAGCAGGTCGAGCACTTCGCTGAC  >  NC_002947/5305518‑5305805

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