Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,246,285	A→G	9.5%	intergenic (+31/+60)	ibpA → / ← PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,246,285	0	A	G	9.5%	169.8 / 5.8	61	intergenic (+31/+60)	ibpA/PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (2/4);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAG  >  NC_002947/2246136‑2246428                                                                                                                                                      |                                                                                                                                                 aaGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCa                                                                                                                                                  >  2:147142/1‑150 (MQ=255)   gTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCACCGAAAGAAAGGGCGCCTCCagc                                                                                                                                                >  1:230700/1‑149 (MQ=255)   gTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCaga                                                                                                                                                <  1:98918/150‑1 (MQ=255)    tggaGGCTGCCGGCCTGGCCACTGGCCTGCTCAGCACCGACCAGCTGCGCATCGTGCCTGAAGAGGCCAAGGCTAAGCGCATTCTGATCAATGGCGAAAAACTCGGGCTCAACTGATCAAACTCCAACGACAAAAAAGGCCCCTCCagag                                                                                                                                               >  2:1559337/4‑150 (MQ=255)    tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                               >  2:1559336/1‑150 (MQ=255)    tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                               <  2:1341009/150‑1 (MQ=255)    tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                               <  1:1523346/150‑1 (MQ=255)    tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                               <  2:58055/150‑1 (MQ=255)      aaGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGt                                                                                                                                             <  2:1461018/150‑1 (MQ=255)          ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                         <  2:1209255/150‑1 (MQ=255)          ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCGGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                         <  1:7533/150‑1 (MQ=255)            gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                       >  1:711849/1‑150 (MQ=255)              cgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACAAAAAAAAGGGCACCTCCAGAGGTGCCCtttt                                                                                                                                     <  2:679827/150‑1 (MQ=255)               gtccTGTCCAATAGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAATCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTTCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCttttt                                                                                                                                    <  2:711849/148‑1 (MQ=255)                gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                   <  2:279415/150‑1 (MQ=255)                  cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttttttc                                                                                                                                 <  2:938864/150‑4 (MQ=255)                   tGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                                <  2:1059806/150‑1 (MQ=255)                    ggCCAATGGCCTGCTCAGCATCGACCTGTTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGc                                                                                                                               <  2:802919/150‑1 (MQ=255)                       cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTtc                                                                                                                            >  1:778009/1‑150 (MQ=255)                         aTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTGTCGTGCTtctc                                                                                                                          <  1:556598/150‑1 (MQ=255)                               tGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑Ag                                                                                                                   >  1:829281/1‑150 (MQ=255)                                   cAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTccc                                                                                                               >  1:909091/1‑150 (MQ=255)                                        tCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTa                                                                                                          <  1:671608/150‑1 (MQ=255)                                           aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGt                                                                                                       >  2:1434384/1‑150 (MQ=255)                                             ctgctgCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTgc                                                                                                     >  1:1307883/1‑150 (MQ=255)                                               gctgcGCATCGGGCCTGAAGAAGCCACGGCCAAGCGCCTTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTgcgc                                                                                                   <  2:230700/150‑1 (MQ=255)                                                   cgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGcc                                                                                               <  2:1307883/150‑1 (MQ=255)                                                    gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCCAC‑CGTGCTAGTCCCGCGTAGGTGCGCGGCCg                                                                                              >  1:820341/1‑150 (MQ=255)                                                      aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTc                                                                                            >  1:172711/1‑150 (MQ=255)                                                         gTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                         <  2:288358/150‑1 (MQ=255)                                                          tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGACAGAAAGGGCACCGCCGGAGGGGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                        <  1:1434384/150‑1 (MQ=255)                                                            ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGa                                                                                      >  2:624621/1‑150 (MQ=255)                                                             cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                     <  2:909091/150‑1 (MQ=255)                                                              tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                    >  1:1447595/1‑150 (MQ=255)                                                                    agCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                              >  1:313239/1‑150 (MQ=255)                                                                     gccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATc                                                                             <  2:412514/150‑1 (MQ=255)                                                                             caagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                     >  2:1032278/1‑150 (MQ=255)                                                                              aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                    <  1:1559336/150‑1 (MQ=255)                                                                              aagCCCATTCTGCTCAATGACCAAAAGCCGCCGCCAAACTGATCGAGCTGCATGGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTGTCGTGCTTCTCTCGC‑AGCCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                    <  1:1559337/150‑1 (MQ=255)                                                                                 cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                 <  2:172413/150‑1 (MQ=255)                                                                                 cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAAATGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCACGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                 <  2:881211/150‑1 (MQ=255)                                                                                  gcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGt                                                                                                       >  2:488882/1‑111 (MQ=255)                                                                                  gcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGt                                                                                                       <  1:488882/111‑1 (MQ=255)                                                                                  gcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTcc                                                                <  2:1447595/150‑1 (MQ=255)                                                                                     ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                              >  1:1080839/1‑133 (MQ=255)                                                                                     ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                              <  2:1080839/133‑1 (MQ=255)                                                                                         gATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                            <  2:629180/131‑1 (MQ=255)                                                                                         gATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                            >  1:629180/1‑131 (MQ=255)                                                                                                   gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCt                                                                                                                                        <  2:1239423/62‑1 (MQ=37)                                                                                                   gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCt                                                                                                                                        >  1:1239423/1‑62 (MQ=37)                                                                                                    aaaaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAAc                                              <  2:1126091/150‑1 (MQ=255)                                                                                                              gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                    <  2:1513742/150‑1 (MQ=255)                                                                                                              gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                    <  1:624621/150‑1 (MQ=255)                                                                                                                       aTCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                         <  1:687895/136‑1 (MQ=255)                                                                                                                       aTCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                         >  2:687895/1‑136 (MQ=255)                                                                                                                              tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                              <  2:264534/140‑1 (MQ=255)                                                                                                                              tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                              >  1:264534/1‑140 (MQ=255)                                                                                                                                cAACGAAAAAAAGGGCACCTCCAGAGGTGCGCTTTTTAGTGCTTCTCACGG‑AGTCACGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGACACGCCGGTAACGGGAATGGCCGGGTATGCa                  >  2:673070/1‑150 (MQ=255)                                                                                                                                    gAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                               <  1:825831/133‑1 (MQ=255)                                                                                                                                    gAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                               >  2:825831/1‑133 (MQ=255)                                                                                                                                      gaaaaaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                               >  2:1017763/2‑115 (MQ=255)                                                                                                                                      gaaaaaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                               <  1:1017763/114‑1 (MQ=255)                                                                                                                                       aaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAcc           >  1:1152336/1‑150 (MQ=255)                                                                                                                                         aaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt         <  2:778009/150‑1 (MQ=255)                                                                                                                                           aGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg       <  2:829281/150‑1 (MQ=255)                                                                                                                                             ggCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTa     >  2:410568/1‑150 (MQ=255)                                                                                                                                                aCCTCCAGAGGTGCCCTATTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAg  >  1:55678/1‑150 (MQ=255)                                                                                                                                                  cTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                    >  1:240596/1‑114 (MQ=255)                                                                                                                                                  cTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                    <  2:240596/114‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGC‑AGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAG  >  NC_002947/2246136‑2246428

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