Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,697,374	Δ1 bp	6.0%	intergenic (+28/‑49)	hslV → / → hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,697,374	0	T	.	6.0%	200.8 / 7.2	67	intergenic (+28/‑49)	hslV/hslU	ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base . (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

GCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCA  >  NC_002947/5697226‑5697516                                                                                                                                                      |                                                                                                                                               gctgAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTg                                                                                                                                               <  2:207284/150‑1 (MQ=255)       cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgggg                                                                                                                                         >  1:322132/1‑150 (MQ=255)       cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgggg                                                                                                                                         <  2:1526977/150‑1 (MQ=255)         aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           <  1:1338552/146‑1 (MQ=255)         aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           >  2:1338552/1‑146 (MQ=255)         aGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCt                                                                                                                                       <  2:187666/150‑1 (MQ=255)           accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCttt                                                                                                                                     <  1:1464166/150‑1 (MQ=255)                   gtcGGCCGGCGAAATCGCCGAAACCGCCCTCAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGGAGGGGCTTTTCCACGCt                                                                                                                             <  2:209398/149‑1 (MQ=255)                           gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGc                                                                                                                              <  2:4769/141‑1 (MQ=255)                           gcgAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGc                                                                                                                              >  1:4769/1‑141 (MQ=255)                               aaTCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGcc                                                                                                                 >  1:1201803/1‑150 (MQ=255)                                   gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           >  2:1632336/1‑120 (MQ=255)                                   gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           <  1:1632336/120‑1 (MQ=255)                                   gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg                                                                                                             <  2:809028/150‑1 (MQ=255)                                       agaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAcc                                                                                                         <  2:632250/150‑1 (MQ=255)                                          ccGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCATCTGCCCAAGGACCGTa                                                                                                      <  1:1498059/150‑1 (MQ=255)                                          ccGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa                                                                                                      <  2:854613/150‑1 (MQ=255)                                                 gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGG‑GAAGGGGcgcc                                                                                                                                     >  2:1490235/1‑109 (MQ=255)                                                 gAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTCACCGG‑GAAGGGGcgcc                                                                                                                                     <  1:1490235/112‑4 (MQ=255)                                                     aTCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                           >  2:182176/1‑150 (MQ=255)                                                      tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           <  2:1051905/101‑1 (MQ=255)                                                      tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAgg                                                                                                                                           >  1:1051905/1‑101 (MQ=255)                                                        gCCGGCGAGATCTGCGTGTTCAACAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAt                                                                                        <  2:892154/150‑1 (MQ=255)                                                         ccGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATg                                                                                       >  1:797532/1‑150 (MQ=255)                                                          cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGa                                                                                      >  1:1283795/1‑150 (MQ=255)                                                          cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGa                                                                                      <  1:182176/150‑1 (MQ=255)                                                            gCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                          <  2:1522133/144‑1 (MQ=255)                                                            gCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTcc                                                                                          >  1:1522133/1‑144 (MQ=255)                                                              gACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAcccc                                                                                  >  1:618073/1‑150 (MQ=255)                                                               aCATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGAccccc                                                                                 >  2:174887/1‑150 (MQ=255)                                                                    tGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            >  1:316916/1‑150 (MQ=255)                                                                      cGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCgaga                                                                          <  2:797532/150‑1 (MQ=255)                                                                         gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCAt                                                                                              <  2:1293628/127‑1 (MQ=255)                                                                         gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCAt                                                                                              >  1:1293628/1‑127 (MQ=255)                                                                           tCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                     >  1:1680956/1‑150 (MQ=255)                                                                            cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTcc                                                                    >  2:1539503/1‑150 (MQ=255)                                                                              ccaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAACGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                                  <  2:322132/150‑1 (MQ=255)                                                                                aaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATGGTCCACGa                                                                >  1:331249/1‑150 (MQ=255)                                                                                          tGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGc                                                      <  2:618073/150‑1 (MQ=255)                                                                                           gACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGcc                                                     >  1:1846125/1‑150 (MQ=255)                                                                                              cctCGAGTAGCAGGACCTGGCCGGCTGATTCAGCTGTTCTGGCTCCCCTTC‑CCGGTGAATGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCTTATTGATCTAATCCAACGCCCCCCGCGAAATCCTACACGAACTCAACCGCcaca                                                  <  2:917899/148‑1 (MQ=255)                                                                                               aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcat                                                 >  2:1302024/1‑150 (MQ=255)                                                                                                  gaggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                              <  2:1846125/150‑1 (MQ=255)                                                                                                    ggagCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCg                                            <  2:1680956/150‑1 (MQ=255)                                                                                                       gCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCaa                                                                                                              <  2:1261257/81‑1 (MQ=255)                                                                                                       gCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCaa                                                                                                              >  1:1261257/1‑81 (MQ=255)                                                                                                           gACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGga                                     <  2:957605/150‑1 (MQ=255)                                                                                                             ccTGGACGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgacg                                   >  2:1536542/1‑150 (MQ=255)                                                                                                                  ccGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGcc                              <  2:1201803/150‑1 (MQ=255)                                                                                                                  ccGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGC‑TTTCCACGCTGGCTCACCTGCCCAAGGACCGTATTGa                                                                                                  <  1:8914/81‑1 (MQ=255)                                                                                                                  ccGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGC‑TTTCCACGCTGGCTCACCTGCCCAAGGACCGTATTGa                                                                                                  >  2:8914/1‑81 (MQ=255)                                                                                                                   cGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATg                                                                                             >  1:1732744/1‑86 (MQ=255)                                                                                                                   cGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATg                                                                                             <  2:1732744/86‑1 (MQ=255)                                                                                                                     aCTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAg                           <  1:1539503/150‑1 (MQ=255)                                                                                                                             aGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            <  2:249532/93‑1 (MQ=255)                                                                                                                             aGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCga                                                                            >  1:249532/1‑93 (MQ=255)                                                                                                                              gCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg                  <  2:1649096/150‑1 (MQ=255)                                                                                                                               cTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAAGGACTGTATTAATCATCTCCTTGTCCCCCCGCGAGGTCGTCCACGAACTCAACCGCCACATCATCTGCCAGGGCGACCCCAAGCGCTCCGTgg                 >  2:647766/1‑150 (MQ=255)                                                                                                                               cTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 <  2:1619634/150‑1 (MQ=255)                                                                                                                                    cTGGCGCTCCTTC‑CCAGTGCAGGGGCTTTTCCACGCTCACTTACCTGCCCAAGGACCGCATTGATCATGCCCATGCCCCCCCGCGAGATCGTCCACCAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt            <  1:634315/150‑1 (MQ=255)                                                                                                                                       gcgcCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgc                          >  1:1821791/1‑133 (MQ=255)                                                                                                                                       gcgcCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgc                          <  2:1821791/133‑1 (MQ=255)                                                                                                                                          ccccTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 <  2:723208/139‑1 (MQ=255)                                                                                                                                          ccccTTCACCGG‑GAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTgg                 >  1:723208/1‑139 (MQ=255)                                                                                                                                           cccTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGTCGCCAAGCGCGCCGTGGCCATTGCCCTgc     <  1:453656/150‑1 (MQ=255)                                                                                                                                             cTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgcgc   >  1:1363237/1‑150 (MQ=255)                                                                                                                                              ttC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCa  >  1:1062639/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                               GCTGAACAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTGTTCTGGCGCCCCTTC‑CCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGCGCA  >  NC_002947/5697226‑5697516

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