Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,805,580	A→G	5.9%	T119A (ACG→GCG)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,805,580	0	A	G	5.9%	194.5 / 2.2	68	T119A (ACG→GCG)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base A (38/26);  new base G (2/2);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TGCACCGGGCCTTGTATGTCGAGGGGATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAGCACCCGC  >  NC_002947/2805437‑2805728                                                                                                                                                  |                                                                                                                                                     tGCACCGGGCCTTGTATGTCGAGGGGATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GAcgccg                                                                                                                                                 >  6:320656/1‑149 (MQ=255) tGCACCGGGCCTTGTATGTCGAGGGGATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GAcgccg                                                                                                                                                 >  3:61020/1‑149 (MQ=255)      cGGGCCTTGTATGTCGAGGGGATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTc                                                                                                                                            >  7:99565/1‑149 (MQ=255)                          gATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GGCGCCGCAGTCGGTGCGGCGCTTCGATGCGt                                                                                                                        >  8:52948/1‑149 (MQ=255)                          gATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GGCGCCGCAGTCGGTGCGGCGCTTCGATGCGt                                                                                                                        >  2:294094/1‑149 (MQ=255)                            tCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTtg                                                                                                                      <  3:266049/149‑1 (MQ=255)                              gACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTtgct                                                                                                                    <  4:61020/149‑1 (MQ=255)                               aCTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTtgctg                                                                                                                   >  1:106389/1‑149 (MQ=255)                                           tctgCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATGACGGGGGATGGTGCCA‑GCGAGAACTGTTTTGTCTTCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAg                                                                                                       <  3:236071/146‑1 (MQ=255)                                             aGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAgg                                                                                                      >  5:172391/1‑148 (MQ=255)                                               cACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTa                                                                                                    <  8:239626/148‑1 (MQ=255)                                               cACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTa                                                                                                    >  7:239626/1‑148 (MQ=255)                                                 ccTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAat                                                                                                  >  8:98233/1‑148 (MQ=255)                                                  cTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAatca                                                                                                >  6:149279/1‑149 (MQ=255)                                                   tGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAatcat                                                                                               >  7:57980/1‑149 (MQ=255)                                                    gTCCAGCCATGTCAAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCACGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAatcatc                                                                                              >  5:11351/1‑149 (MQ=255)                                                     tCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCg                                                                                                         >  8:77581/1‑137 (MQ=255)                                                     tCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCg                                                                                                         <  7:77581/137‑1 (MQ=255)                                                       cAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGc                                                                                           <  6:172391/149‑1 (MQ=255)                                                       cAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGGCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTg                                                                                            <  8:26081/148‑1 (MQ=255)                                                        aGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGc                                                                                           >  3:42182/1‑148 (MQ=255)                                                          ccATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAg                                                                                        >  2:69530/1‑149 (MQ=255)                                                             tGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTg                                                                                     <  1:69530/149‑1 (MQ=255)                                                               tCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGc                                                                                           >  8:254910/1‑141 (MQ=255)                                                               tCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGc                                                                                           <  7:254910/141‑1 (MQ=255)                                                                 cAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGtt                                                                                      >  4:168709/1‑144 (MQ=255)                                                                 cAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGtt                                                                                      <  3:168709/144‑1 (MQ=255)                                                                 cAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGaa                                                                                  <  2:106389/148‑1 (MQ=255)                                                                 cAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAg                                                                                 >  3:248512/1‑149 (MQ=255)                                                                  aGCGGTGTGGCGCTGATCACGGTGGATGCTGCCAGGCCACAACTGAATTGTCATCATACCCGGCGCCAACGGCCTGCT‑CACGCCGCAGTCGGTgcgg                                                                                                                                     >  1:54535/1‑97 (MQ=255)                                                                  aGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTTTTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTgcgg                                                                                                                                     <  2:54535/96‑1 (MQ=255)                                                                           gcgcTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGccagccag                                                                       >  7:256040/1‑149 (MQ=255)                                                                                gATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCg                                                                  >  8:110925/1‑149 (MQ=255)                                                                                 aTCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCAcc                                                                   <  8:256040/147‑1 (MQ=255)                                                                                     cGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTgg                                                                                    <  8:110122/126‑1 (MQ=255)                                                                                     cGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTgg                                                                                    >  7:110122/1‑126 (MQ=255)                                                                                       gTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTg                                                           <  4:248512/149‑1 (MQ=255)                                                                                          gATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GGCGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGac                                                        <  7:52948/149‑1 (MQ=255)                                                                                          gATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑CGCGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGac                                                        <  1:294094/149‑1 (MQ=255)                                                                                              ctgcCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTg                                                    <  4:42182/149‑1 (MQ=255)                                                                                               tgcCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTgg                                                   >  6:76448/1‑149 (MQ=255)                                                                                                  ca‑gccagAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGc                                                  >  8:128570/1‑147 (MQ=255)                                                                                                      ccagAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCagag                                             >  3:188181/1‑149 (MQ=255)                                                                                                      ccagAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCagag                                             >  3:304878/1‑149 (MQ=255)                                                                                                       cagAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAgg                                            <  8:57980/149‑1 (MQ=255)                                                                                                       cagAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAgg                                            >  1:339630/1‑149 (MQ=255)                                                                                                        agAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GAcgccg                                                                                                                                                 >  3:190701/1‑47 (MQ=255)                                                                                                        agAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GAcgccg                                                                                                                                                 <  4:190701/47‑1 (MQ=255)                                                                                                          aaCTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCa                                         <  4:188181/149‑1 (MQ=255)                                                                                                                 ttGTCATCATTCCCGGCGCCAACGGCCTGCTAGACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCt                                   >  1:220663/1‑149 (MQ=255)                                                                                                                   gTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTggg                                >  4:184559/1‑149 (MQ=255)                                                                                                                   gTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTgg                                 <  7:110925/148‑1 (MQ=255)                                                                                                                    tcatcaTTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGc                               <  5:76448/149‑1 (MQ=255)                                                                                                                              cGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTa                                                                                                    >  4:77298/1‑70 (MQ=255)                                                                                                                              cGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTa                                                                                                    <  3:77298/70‑1 (MQ=255)                                                                                                                              cGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGa                     >  5:259252/1‑149 (MQ=255)                                                                                                                                gcgcCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATc                   >  4:322769/1‑149 (MQ=255)                                                                                                                                     aaCGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTtg                                                                                                                      >  4:29981/1‑45 (MQ=255)                                                                                                                                     aaCGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTtg                                                                                                                      <  3:29981/45‑1 (MQ=255)                                                                                                                                     aaCGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAgg                                            <  8:183895/119‑1 (MQ=255)                                                                                                                                     aaCGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAgg                                            >  7:183895/1‑119 (MQ=255)                                                                                                                                          ccTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAg         >  2:230979/1‑149 (MQ=255)                                                                                                                                            tgct‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGc                                    >  6:272745/1‑120 (MQ=255)                                                                                                                                            tgct‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGc                                    <  5:272745/120‑1 (MQ=255)                                                                                                                                              ct‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAGCAcc     >  8:13861/1‑149 (MQ=255)                                                                                                                                               t‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATc                   >  8:148614/1‑134 (MQ=255)                                                                                                                                               t‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATc                   <  7:148614/134‑1 (MQ=255)                                                                                                                                               t‑gACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAGCAccc    >  5:114658/1‑149 (MQ=255)                                                                                                                                                  aCGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAGCACCCGc  >  7:101035/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGCACCGGGCCTTGTATGTCGAGGGGATCGACTGTCAGGGCGTCAGCACCTGTCCAGCCATGTCCAGCGGTGTGGCGCTGATCACGGTGGATGCTGCCA‑GCCAGAACTGTATTGTCATCATTCCCGGCGCCAACGGCCTGCT‑GACGCCGCAGTCGGTGCGGCGCTTCGATGCGTTGCTGCAGGCTGCCGAGGTAATCATCTGCCAGTTGGAAGTGCCAGCCAGCACCGTGGCCTGGACACTGGCCAGAGGGCACGAGCTGGGCAAGCAGGTGATCCTGAACCCAGCACCCGC  >  NC_002947/2805437‑2805728

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