Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,896,312	A→G	5.8%	intergenic (‑137/+48)	cysU ← / ← sbp‑II	sulfate ABC transporter permease CysT/sulfate ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,896,312	0	A	G	5.8%	249.4 / 3.1	85	intergenic (‑137/+48)	cysU/sbp‑II	sulfate ABC transporter permease CysT/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base G (2/3);  total (39/46)
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.80e-01

CGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGT  >  NC_002947/5896168‑5896460                                                                                                                                                  |                                                                                                                                                      cgcgACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGTCGCTTCCACAAGGAACCATGATCTCCTGTGGGAGCGGCCtt                                                                                                                                                   <  7:164869/149‑1 (MQ=21) cgcgACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCtt                                                                                                                                                   <  3:302772/149‑1 (MQ=255)  gcgACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTtg                                                                                                                                                  >  1:196349/1‑149 (MQ=255)     aCATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGc                                                                                                                                               <  5:334683/149‑1 (MQ=255)       aTAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCcg                                                                                                                                             >  5:138748/1‑149 (MQ=255)          aaTATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCcgcg                                                                                                                                           >  1:139473/1‑148 (MQ=255)              tCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATgg                                                                                                                                      <  8:232723/149‑1 (MQ=255)                ccTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGTGGGCTGCTTGGCAGCCCATTCGCGACACAAGCCCCCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATggg                                                                                                                                     >  1:259434/1‑149 (MQ=255)                ccTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGc                                                                                                                                    >  8:100124/1‑149 (MQ=255)                 cTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGcc                                                                                                                                   >  2:195424/1‑149 (MQ=255)                 cTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGAACATGGTTCCTTGTGGAAGCGAc                                                                                                                                                      >  5:279120/1‑130 (MQ=255)                 cTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGAACATGGTTCCTTGTGGAAGCGAc                                                                                                                                                      <  6:279120/130‑1 (MQ=255)                    gTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑c                                                                                                                                <  1:250412/148‑1 (MQ=255)                       gAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATggg                                                                                                                                     <  6:266860/141‑1 (MQ=255)                       gAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATggg                                                                                                                                     >  5:266860/1‑141 (MQ=255)                        aaCGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CAcgc                                                                                                                            <  1:195424/148‑1 (MQ=255)                         aCGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCgg                                                                                                                          <  6:138748/149‑1 (MQ=255)                           gCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGcc                                                                                                                        >  4:131305/1‑149 (MQ=255)                                aaGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCt                                                                                                                   >  6:81615/1‑149 (MQ=255)                                   cccGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTc                                                                                                                >  1:294843/1‑149 (MQ=255)                                   cccGGAAACAATCCCCGGCCGGTGAAAGCTTTTGAATCCTG‑TAGCTGCTTGGCAGCCCATTCGCGACACAAAGCCCCACCCACATGACAACATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCCGCCCCGCATTc                                                                                                                >  2:121246/1‑149 (MQ=255)                                      ggAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCtt                                                                                                             <  7:55487/149‑1 (MQ=255)                                       gAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCttt                                                                                                            >  4:171746/1‑149 (MQ=255)                                          acacTCCCGGGCCGGGGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGgc                                                                                                         <  7:229077/149‑1 (MQ=255)                                            aaTCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATAATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGgcgc                                                                                                       >  2:316024/1‑149 (MQ=255)                                             aTCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCt                                                                                                      <  7:174459/149‑1 (MQ=255)                                                ccGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGgcgc                                                                                                       >  3:78662/1‑145 (MQ=255)                                                ccGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGgcgc                                                                                                       <  4:78662/145‑1 (MQ=255)                                                 cGGGCCGGTGCAAGCTTTTGAATCCTC‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGAGCCGCGAATGGGCCG‑CACGCGGCCCCGCATTCCTTTGGCGCTTACt                                                                                                  >  5:145460/1‑149 (MQ=255)                                                  gggCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACtg                                                                                                 >  2:101349/1‑149 (MQ=255)                                                   ggCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACtgt                                                                                                <  1:308399/149‑1 (MQ=255)                                                   ggCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACtgt                                                                                                <  3:168122/149‑1 (MQ=255)                                                        gTGCAAGCTTTTGAAACCTG‑GGGCTGCTTGCCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCCTGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTg                                                                                           <  6:104202/149‑1 (MQ=255)                                                           cAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGGTCCCACAGGAGATCATGGTTCCTTGTTGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGAt                                                                                         <  1:121246/148‑1 (MQ=255)                                                           cAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGAt                                                                                         <  2:139473/148‑1 (MQ=255)                                                               cTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATc                                                                                    <  7:100124/149‑1 (MQ=255)                                                                ttttGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCt                                                                                   >  4:250329/1‑149 (MQ=255)                                                                ttttGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCt                                                                                   >  1:101851/1‑149 (MQ=255)                                                                  ttGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTg                                                                                  <  3:149020/148‑1 (MQ=255)                                                                     aaTCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTc                                                                               >  1:236148/1‑148 (MQ=255)                                                                       tCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGaa                                                                            <  2:236148/149‑1 (MQ=255)                                                                              ggCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACAcc                                                                      <  2:294843/149‑1 (MQ=255)                                                                                       gCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCAtt                                                             >  5:326109/1‑149 (MQ=255)                                                                                        cAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGGCCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTAATGTGCCTGATAGATCTGGTCGAACACACCGCCATCAtt                                                             >  1:120316/1‑148 (MQ=255)                                                                                        cAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTgg                                                                                 <  5:272847/128‑1 (MQ=255)                                                                                         aGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCAtt                                                             <  8:53442/147‑1 (MQ=255)                                                                                         aGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCAtt                                                             >  7:53442/1‑147 (MQ=255)                                                                                         aGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTga                                                           <  3:297548/149‑1 (MQ=255)                                                                                          gCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTgaa                                                          <  3:171746/149‑1 (MQ=255)                                                                                          gCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTgaa                                                          <  3:108373/149‑1 (MQ=255)                                                                                             cATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTc                                                                               >  5:247808/1‑125 (MQ=255)                                                                                             cATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTc                                                                               <  6:247808/125‑1 (MQ=255)                                                                                               ttCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAAtt                                                     >  2:324062/1‑149 (MQ=255)                                                                                                  gcgACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTc                                                   <  1:101349/148‑1 (MQ=255)                                                                                                  gcgACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTc                                                   <  3:157554/148‑1 (MQ=255)                                                                                                   ccccacaCGGCCGCACCCACAGGAGTTCATGGTTCCTTGTGTAAGCGACATTGTGCCGCTAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCgg                                                 <  2:120316/146‑1 (MQ=255)                                                                                                    ggcacaAGGTCGCTTCCACAAGGAACCATGATCTCCTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGt                                                >  7:292072/3‑149 (MQ=18)                                                                                                     gcacaAGGTCGCTTCCACAAGGAACCATGATCTCCTGGGGGAGCGGCCTTGTGTCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGtt                                               <  8:292072/148‑1 (MQ=18)                                                                                                          aGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCa                                          <  1:324062/149‑1 (MQ=255)                                                                                                           ggCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAg                                         <  5:81615/149‑1 (MQ=255)                                                                                                              cGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTTCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCt                                      >  1:92635/1‑149 (MQ=255)                                                                                                                cTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTc                                    >  1:59158/1‑149 (MQ=255)                                                                                                                   ccACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccac                                 >  4:191090/1‑149 (MQ=255)                                                                                                                    cacaGGAGATCATGGTTCCTTGTGGAAGCGACATTGTGCCGCGAATGGGCAG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccacc                                >  1:145022/1‑149 (MQ=255)                                                                                                                     acaGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccaccc                               <  8:286664/149‑1 (MQ=255)                                                                                                                       aGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccacccac                             >  4:242577/1‑149 (MQ=255)                                                                                                                               aTGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTcc                     <  3:250329/149‑1 (MQ=255)                                                                                                                                tGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccaccca                              <  6:61655/139‑1 (MQ=255)                                                                                                                                tGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTccaccca                              >  5:61655/1‑139 (MQ=255)                                                                                                                                   ttCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGt                 <  2:101851/149‑1 (MQ=255)                                                                                                                                   tgccTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGt                 <  2:196349/147‑1 (MQ=255)                                                                                                                                       ttGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCCCCAAAGTCCTTGTCGAt             >  8:121889/1‑149 (MQ=255)                                                                                                                                       ttGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCaaa                         >  7:12136/1‑137 (MQ=255)                                                                                                                                       ttGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCaaa                         <  8:12136/137‑1 (MQ=255)                                                                                                                                       ttGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGAt             >  7:111862/1‑149 (MQ=255)                                                                                                                                         gtgGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTTCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATg            <  2:92635/148‑1 (MQ=255)                                                                                                                                          tgGGAGCGGCCTTGTGTCGCGAATGGGCTGCCAAGCAGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATg            <  6:257411/148‑1 (MQ=38)                                                                                                                                            gAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCa        >  7:308176/1‑149 (MQ=255)                                                                                                                                             aaGCGAGCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGATTTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCAc       <  6:326109/149‑1 (MQ=255)                                                                                                                                             aaGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCAc       <  3:191090/149‑1 (MQ=255)                                                                                                                                             aaGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCAc       >  7:32788/1‑149 (MQ=255)                                                                                                                                              aGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTcc                     <  1:298331/134‑1 (MQ=255)                                                                                                                                              aGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTcc                     >  2:298331/1‑134 (MQ=255)                                                                                                                                                cGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAg    <  8:111862/149‑1 (MQ=255)                                                                                                                                                 gACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAgg   >  3:218423/1‑149 (MQ=255)                                                                                                                                                  aCCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGt  >  6:288008/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CGCGACATAAATATCCCTGGTTGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTG‑GGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCG‑CACGCGGCCCCGCTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGT  >  NC_002947/5896168‑5896460

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