Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,896,391	A→G	7.0%	intergenic (‑160/+25)	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,391	0	A	G	7.0%	203.6 / 4.6	71	intergenic (‑160/+25)	cysU/sbp‑II	sulfate ABC transporter permease CysT/sulfate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/37);  new base G (3/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

TGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCGGGAAC  >  NC_002947/5896245‑5896530                                                                                                                                                   |                                                                                                                                             tGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCAcgc                                                                                                                                           >  2:952813/1‑150 (MQ=255)    aCGCACTAAGCCCGGAAACAATCGCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCACACGCGGc                                                                                                                                        <  1:16229/150‑1 (MQ=255)       cACTAAGCCCGGAAACAATCCCGGGCCCGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGAGATACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGcccc                                                                                                                                     <  1:1240910/150‑1 (MQ=255)        aCTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGcccc                                                                                                                                     <  2:1069088/149‑1 (MQ=255)        aCTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGcccc                                                                                                                                     >  1:1069088/1‑149 (MQ=255)        aCTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCg                                                                                                                                    >  2:104914/1‑150 (MQ=255)          tAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑Ct                                                                                                                                 >  1:542541/1‑150 (MQ=255)                 ggAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCAcgc                                                                                                                                           <  1:585545/134‑1 (MQ=255)                 ggAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCAcgc                                                                                                                                           >  2:585545/1‑134 (MQ=255)                  gAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTg                                                                                                                         <  1:540076/150‑1 (MQ=255)                  gAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTg                                                                                                                         <  1:1412223/150‑1 (MQ=255)                       aaTCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCGTGCGGcccattcgc                                                                                                                                <  2:1363606/139‑7 (MQ=39)                       aaTCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCGTGCGGcccattcgc                                                                                                                                >  1:1363606/1‑133 (MQ=39)                       aaTCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTgg                                                                                                                        <  2:1270447/146‑1 (MQ=255)                       aaTCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTgg                                                                                                                        >  1:1270447/1‑146 (MQ=255)                         tCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTa                                                                                                                  <  2:542541/150‑1 (MQ=255)                                 cGGTGCAAGCTTTTGAATCCTGGGGCTTCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGACTGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACGGGGCCt                                                                                                          <  2:132345/150‑1 (MQ=255)                                    tGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGAt                                                                                                       >  2:1292031/1‑150 (MQ=255)                                       aaGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGa                                                                                                        >  2:928992/1‑146 (MQ=255)                                       aaGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGa                                                                                                        <  1:928992/146‑1 (MQ=255)                                           ttttGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTg                                                                                                >  1:827403/1‑150 (MQ=255)                                           ttttGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTg                                                                                                <  1:97464/150‑1 (MQ=255)                                           ttttGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTg                                                                                                >  2:530171/1‑150 (MQ=255)                                           ttttGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTg                                                                                                >  2:533389/1‑150 (MQ=255)                                              tGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTc                                                                                             <  2:861364/150‑1 (MQ=255)                                                  tCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑Ct                                                                                                                                 >  1:213660/1‑110 (MQ=255)                                                  tCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑Ct                                                                                                                                 <  2:213660/110‑1 (MQ=255)                                                  tCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAac                                                                                         >  2:632325/1‑150 (MQ=255)                                                   ccTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAaca                                                                                        <  2:827403/150‑1 (MQ=255)                                                       gggCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACAcc                                                                                    >  2:85455/1‑150 (MQ=255)                                                        gatctctTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGCGATCTTGGTTCCTTGCGGAAGCAACCTTGTGCCGCGAATGCGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCg                                                                                   <  2:518934/145‑1 (MQ=255)                                                          ctgctTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGcc                                                                                 >  1:371396/1‑150 (MQ=255)                                                              ttGGCAGCCTATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑GTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGcagatc                                                                             >  1:518934/1‑145 (MQ=255)                                                              ttGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCcatca                                                                             <  2:371396/150‑1 (MQ=255)                                                                ggCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCAtt                                                                           <  2:1387898/150‑1 (MQ=255)                                                                  cAGCCCATTCGCGACACAAGGCCGCTCCCCCAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTga                                                                         <  2:690433/150‑1 (MQ=255)                                                                   aGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTgaa                                                                        <  1:1035262/150‑1 (MQ=255)                                                                       cATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCca                                                                                >  1:470098/1‑138 (MQ=255)                                                                       cATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCca                                                                                <  2:470098/138‑1 (MQ=255)                                                                       cATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAAt                                                                    >  1:858832/1‑150 (MQ=255)                                                                              gACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATc                                                                                                  <  2:568028/113‑1 (MQ=255)                                                                              gACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATc                                                                                                  >  1:568028/1‑113 (MQ=255)                                                                               acacAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTtg                                                            <  1:781109/150‑1 (MQ=255)                                                                                   aaGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGGCATCAGTGAAGAATTTCGGTTGTGCa                                                        >  2:330740/1‑150 (MQ=255)                                                                                   aaGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCa                                                        >  2:1038819/1‑150 (MQ=255)                                                                                   aaGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCATTGCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCa                                                        >  1:727757/1‑150 (MQ=255)                                                                                    aGGCTGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCGGCTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAg                                                       <  1:330740/150‑1 (MQ=255)                                                                                           tCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTcca                                                <  1:530171/150‑1 (MQ=255)                                                                                               acaGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTgaagaa                                                                     >  1:344897/1‑125 (MQ=255)                                                                                               acaGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTgaagaa                                                                     <  2:344897/125‑1 (MQ=255)                                                                                               acaGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGCCAGTGCAGTCTTccacctt                                            <  2:1022183/150‑3 (MQ=255)                                                                                                      aTCATGGTGCCTTGTGGAAGCGACCTTGTGCCGCGAATGCGCCGCACGCTGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCAGGTCGAACACACCGCCATCATTGAAGACTTTCGGTTGTGCATTCTTCCACCCACCAAAGt                                     >  2:1350627/1‑150 (MQ=255)                                                                                                       tCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTg                                                                                                                         <  2:821769/65‑1 (MQ=255)                                                                                                       tCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTg                                                                                                                         >  1:821769/1‑65 (MQ=255)                                                                                                        cATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTcc                                   >  1:948289/1‑150 (MQ=255)                                                                                                         aTGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCt                                                                                                                    >  2:1361847/1‑68 (MQ=255)                                                                                                         aTGGTTCCTTGTGGAAGCGACCTTGTGACGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCt                                                                                                                    <  1:1361847/68‑1 (MQ=255)                                                                                                            tgtgcgtgtgGAAGCGGACGTGTGCCGCGAATGCGTTGCA‑TCGGCCGCGACTTCCCGTTGGCTCGTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGt                               <  2:727757/144‑1 (MQ=255)                                                                                                                cTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGAt                           <  2:858832/150‑1 (MQ=255)                                                                                                                cTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGAt                           <  2:1382676/150‑1 (MQ=255)                                                                                                                        aGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCa                   <  2:948289/150‑1 (MQ=255)                                                                                                                            aCCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCANTCTTCCACCCACNAAAGTCCTTGTCGATGGTCACCAGGTc               <  1:1038819/150‑1 (MQ=255)                                                                                                                               ttGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAg            >  1:921824/1‑150 (MQ=255)                                                                                                                                    ccGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCg       <  1:952813/150‑1 (MQ=255)                                                                                                                                     cgcgAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCgg      <  1:104914/150‑1 (MQ=255)                                                                                                                                       aagcggggCCGCGTGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCGGGa    >  1:617614/6‑150 (MQ=255)                                                                                                                                         gcggggCCGCGTGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTc                                                  >  1:1421578/4‑102 (MQ=255)                                                                                                                                         gcggggCCGCGTGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTc                                                  <  2:1421578/99‑1 (MQ=255)                                                                                                                                         aaTGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCg       <  2:1104418/145‑1 (MQ=255)                                                                                                                                         aaTGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCg       >  1:1104418/1‑145 (MQ=255)                                                                                                                                         aaTGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCGGGAAc  <  2:1199517/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             TGAACGCACTAAGCCCGGAAACAATCCCGGGCCGGTGCAAGCTTTTGAATCCTGGGGCTGCTTGGCAGCCCATTCGCGACACAAGGCCGCTCCCACAGGAGATCATGGTTCCTTGTGGAAGCGACCTTGTGCCGCGAATGGGCCGCACGCGGCCCCG‑CTTTCCTTTGGCGCTTACTGTGCCTGATAGATCTGGTCGAACACACCGCCATCATTGAAGAATTTCGGTTGTGCAGTCTTCCACCCACCAAAGTCCTTGTCGATGGTCACCAGGTCCAGTTTCGGGAAC  >  NC_002947/5896245‑5896530

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑