Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,395,344	A→G	5.7%	E762E (GAA→GAG)	hsdR →	type I restriction modification system endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,395,344	0	A	G	5.7%	258.7 / 3.8	87	E762E (GAA→GAG)	hsdR	type I restriction modification system endonuclease
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base G (2/3);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAAT  >  NC_002947/5395196‑5395489                                                                                                                                                       |                                                                                                                                                   aGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGaa                                                                                                                                                    <  2:251047/149‑1 (MQ=255)  gCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAAc                                                                                                                                                   <  8:234406/149‑1 (MQ=255)     tAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTgg                                                                                                                                                >  5:60032/1‑149 (MQ=255)       actactTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGt                                                                                                                                               >  1:372412/1‑148 (MQ=255)          actTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTa                                                                                                                                             >  4:297281/1‑147 (MQ=255)           ctTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTaaa                                                                                                                                           <  4:303547/148‑1 (MQ=255)              cGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAggg                                                                                                                                        <  4:132305/148‑3 (MQ=38)                ccAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGCGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAgggatc                                                                                                                                     <  1:95366/149‑6 (MQ=37)                 cAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGTTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTg                                                                                                                                    >  8:167344/1‑149 (MQ=255)                   aaTACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGTACCCATCGAACTGGTTAAAGCCTTTggc                                                                                                                                  <  1:49625/149‑1 (MQ=255)                     tACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcgg                                                                                                                                <  1:14608/149‑1 (MQ=255)                       ctacGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACTCCGTCAACGGCTTTGAGCACATCGAGGAAATC‑‑TAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTTCACCCATCGAACTGGTTAAAGCCTTTGGCGTCa                                                                                                                              <  2:241219/146‑1 (MQ=255)                       cTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTggcggc                                                                                                                               <  6:29851/148‑1 (MQ=255)                        tCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCa                                                                                                                              >  7:53363/1‑148 (MQ=255)                             tGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGTTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTg                                                                                                                       <  7:167344/149‑1 (MQ=255)                               cAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTgg                                                                                                                      <  7:327895/148‑1 (MQ=255)                               cAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTgg                                                                                                                      <  1:293605/148‑1 (MQ=255)                                  cgcgcCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTa                                                                                                                   >  6:15291/1‑148 (MQ=255)                                        aGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACaaa                                                                                                            <  3:68933/149‑1 (MQ=255)                                              tGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATc                                                                                                      >  5:134903/1‑149 (MQ=255)                                              tGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATc                                                                                                      >  1:288496/1‑149 (MQ=255)                                                gAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCa                                                                                                    <  8:53363/149‑1 (MQ=255)                                                 aaGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCAt                                                                                                   <  2:203954/149‑1 (MQ=255)                                                         ctgGATAAATACGCCGACACCGGCATTGAGCACATCGATGTCCGCATAAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAAc                                                                                               <  1:329032/147‑1 (MQ=255)                                                         ctgGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTgg                                                                                            <  6:24938/148‑1 (MQ=255)                                                         ctgGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTg                                                                                             <  2:288496/147‑1 (MQ=255)                                                         ctgGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGa                                                                                           <  5:259073/149‑1 (MQ=255)                                                         ctgGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCATTAACCAGCTCAGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTgg                                                                                            <  4:229414/148‑1 (MQ=255)                                                          ctgATAAATACGCCGACACCGGCATTTAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCAACCATCGAACTGGTTAAAGCCTTTTGCGGCAAAGGCTGGCTACAACAAAGCTATCCATGAACTGGa                                                                                           <  1:47486/147‑1 (MQ=255)                                                            gATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAg                                                                                         <  7:180520/148‑1 (MQ=255)                                                                  tACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTTACCAGCTCGGTGAACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCt                                                                                  >  4:252493/1‑149 (MQ=255)                                                                  tACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCt                                                                                  >  5:223066/1‑149 (MQ=255)                                                                  tACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCt                                                                                  <  5:15291/149‑1 (MQ=255)                                                                  tACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCt                                                                                  >  5:13750/1‑149 (MQ=255)                                                                     gCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTa                                                                               >  8:337995/1‑149 (MQ=255)                                                                     gCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTa                                                                               >  1:94281/1‑149 (MQ=255)                                                                     gCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTa                                                                               >  4:267158/1‑149 (MQ=255)                                                                      ccGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACACATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTAc                                                                              >  6:314824/1‑149 (MQ=255)                                                                           acCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAg                                                                         >  5:67874/1‑149 (MQ=255)                                                                             cGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGAcgcg                                                                                    >  5:14610/1‑136 (MQ=255)                                                                             cGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGAcgcg                                                                                    <  6:14610/136‑1 (MQ=255)                                                                             cGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGt                                                                        >  4:18403/1‑148 (MQ=255)                                                                                      gCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGctc                                                              <  7:109907/149‑1 (MQ=255)                                                                                       cacaTCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGctc                                                              >  6:319334/1‑148 (MQ=255)                                                                                          aTCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTAcc                                                          >  7:203306/1‑149 (MQ=255)                                                                                                 aCATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtggg                                                   >  7:137295/1‑149 (MQ=255)                                                                                                   aTC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtg                                                 <  2:372412/149‑1 (MQ=255)                                                                                                   aTC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGAtgggtg                                                 >  3:161967/1‑149 (MQ=255)                                                                                                   aTC‑‑AAAATCCTCCAGCACGATCCCTTTAACCAGCTCGGTGCACCCATAGAACTGGTTAAAGCATTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAACACGCGCTGTACGCCAGTTAATCCTCTCTACCCGTtggctg                                                 >  1:91089/1‑149 (MQ=255)                                                                                                     c‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtggg                                               >  5:157744/1‑149 (MQ=255)                                                                                                          aaTCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCaa                                            <  5:319334/149‑1 (MQ=255)                                                                                                           aTCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAAc                                           <  3:245497/149‑1 (MQ=255)                                                                                                             cctccAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAAcgg                                         >  5:292867/1‑149 (MQ=255)                                                                                                                 cAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACCACAAGGCTATCCATGAACTGGAAGCCGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAAcggcg                                       >  4:330223/1‑147 (MQ=255)                                                                                                                  aGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAg                                    <  3:18403/149‑1 (MQ=255)                                                                                                                   gCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtggg                                               >  4:31398/1‑137 (MQ=255)                                                                                                                   gCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTtgggtggg                                               <  3:31398/137‑1 (MQ=255)                                                                                                                    cTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCa                                                                          >  4:4346/1‑109 (MQ=255)                                                                                                                    cTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCa                                                                          <  3:4346/109‑1 (MQ=255)                                                                                                                       gATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAGAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTg                               >  2:340818/1‑149 (MQ=255)                                                                                                                       gATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTg                               >  8:292834/1‑149 (MQ=255)                                                                                                                       gATCCCTTTAACCAGATCGGTGCACCCAACGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAAAAAAGATATCCATGAACTGCAAGACGCGCTGTACGCCAGTTAAACCGCTCTACCCGTTGGGTGGCCAACGCCGGAGCATTg                               >  2:388807/1‑149 (MQ=255)                                                                                                                           ccTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTGCCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCa                           >  1:351289/1‑149 (MQ=255)                                                                                                                           ccTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCa                           >  2:359659/1‑149 (MQ=255)                                                                                                                            cTTTAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCAc                          >  4:160836/1‑149 (MQ=38)                                                                                                                            cTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCAc                          <  6:67874/149‑1 (MQ=255)                                                                                                                             tttAACCAGTTCGATGGGTGCACCGAGCTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCAcc                         >  8:197669/1‑149 (MQ=38)                                                                                                                             tttAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCAcc                         >  5:68243/1‑149 (MQ=255)                                                                                                                                aaCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCt                        >  8:186272/1‑147 (MQ=255)                                                                                                                                 aCCAGCTCGGTGCACCCATCGAACTGTTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTaa                      <  3:297281/148‑1 (MQ=255)                                                                                                                                 aCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTaa                      <  2:94281/148‑1 (MQ=255)                                                                                                                                 aCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAag                     >  5:26703/1‑149 (MQ=255)                                                                                                                                 aCCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAag                     >  3:59862/1‑149 (MQ=255)                                                                                                                                   cAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAagca                   <  7:337995/149‑1 (MQ=255)                                                                                                                                      cTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGTAACGGCGGAGCCTTGTCCACCTAAGCAGcc                <  6:60032/149‑1 (MQ=255)                                                                                                                                              aCCCTTCTAATTTGTTAATGCTTTTGGCGGCAAA‑GCTTGCTACAACAAAGCTATCCATGAAGTGGAAGACGCGCTGTACGCCAGTTAATGCGCTGTACGCGTTGGGTGGGCATCGGCGTAGCCTTGTCCACCTAAGCAGCCAAGTCGAt        <  2:91089/149‑1 (MQ=255)                                                                                                                                              aCCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCCACTCCAt        >  6:43348/1‑149 (MQ=255)                                                                                                                                              aCCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCAt        >  1:346825/1‑149 (MQ=255)                                                                                                                                              aCCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCAt        >  1:134410/1‑149 (MQ=255)                                                                                                                                                ccATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTa      <  2:115465/149‑1 (MQ=255)                                                                                                                                                 cATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTAc     <  1:359659/149‑1 (MQ=255)                                                                                                                                                  aTCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACa    <  4:59862/149‑1 (MQ=255)                                                                                                                                                   tCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACaa   <  2:346825/149‑1 (MQ=255)                                                                                                                                                   tCGAACTGGTTAAAGCCTTTGGAGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACaa   <  4:104561/149‑1 (MQ=255)                                                                                                                                                   ccGAGCTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACaa   <  3:160836/148‑1 (MQ=255)                                                                                                                                                    cGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAAt  >  5:133583/1‑149 (MQ=255)                                                                                                                                                     gAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTCGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAagca                   <  6:116927/131‑1 (MQ=255)                                                                                                                                                     gAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTAGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAagca                   >  5:116927/1‑131 (MQ=255)                                                                                                                                                       |                                                                                                                                                   AGCGTAACTACTTCGCCAAATACTCCGGTGCAGCGCGCCAGGTGCTGGAAGCGCTGCTGGATAAATACGCCGACACCGGCATTGAGCACATCGAGGACATC‑‑AAAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAA‑GCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCCTTGTCCACCTAAGCAGCCAACTCCATTACAAT  >  NC_002947/5395196‑5395489

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