Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,849,941	G→A	intergenic (‑92/+33)	hycB ← / ← hycA	formate hydrogenlyase subunit HycB/regulator of the transcriptional regulator FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,849,941	0	G	A	97.9%	147.1 / NA	48	intergenic (‑92/+33)	hycB/hycA	formate hydrogenlyase subunit HycB/regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/25);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGT  >  NC_000913/2849698‑2850173                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         cAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACg                                                                                                                                                                                                                                   >  1:404841/1‑251 (MQ=255)                              caGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAt                                                                                                                                                                                                      <  1:288361/251‑1 (MQ=255)                                  ctCACTCTCAGGCGCGGCATTGATTGCAGGCGGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACCGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGCCGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATACCCAAACAGACGCAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgctg                                                                                                                                                                                                  >  1:361459/2‑251 (MQ=255)                                  atccctctCCGGCGCGGCATTGATTGCAGGCCGTGTTGGCGATGCGTCTTTGAACAGGCGGCCTCACAAGTGTGGCCGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAGGCCTCAGCTCAtgctg                                                                                                                                                                                                  <  2:387184/247‑1 (MQ=255)                                  aTCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCctgctg                                                                                                                                                                                                  >  2:391076/1‑251 (MQ=255)                                  aTCACTCTCAGGAGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgctg                                                                                                                                                                                                  >  1:361460/1‑251 (MQ=255)                                            ggCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAAc                                                                                                                                                                                                                         >  2:71309/1‑218 (MQ=255)                                            ggCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAAc                                                                                                                                                                                                                         <  1:71309/218‑1 (MQ=255)                                                    attgattgCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTa                                                                                                                                                                                <  1:218614/251‑1 (MQ=255)                                                          ggCAGGCCGTGCTGGCGATGCGTCTCGTAACAGGCGGCCTAAAAAGTTTGGCAGCCGATAAAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTGCGACACTCATCGACACGCCCAGCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGtt                                                                                                                                                                          <  2:414740/250‑1 (MQ=255)                                                           gatggCCGGGCTGGCGATGGGTCTCTGAACTGGCGGCCTCACAAGTGTGGCAGCCGTTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCAGTTTCGACACTCATCGACACGCCCATCCCCAAACATACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTg                                                                                                                                                                         <  1:460551/248‑1 (MQ=255)                                                               gCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCAtgctgc                                                                                                                                                                                                 >  1:233604/1‑223 (MQ=255)                                                                ccggtcTGGCGATGCGTCTCTGAACAGGCGGCCACACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTa                                                                                                                                                                    <  2:404841/246‑1 (MQ=255)                                                                           aTGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGc                                                                                                                                                         >  1:12196/1‑251 (MQ=255)                                                                              cGTCTCTGAACAGGCGGCCTCACAAGGGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCGCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGcaca                                                                                                                                                      <  1:328881/251‑1 (MQ=255)                                                                                           gcggcCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAg                                                                                                                                         <  1:337501/251‑1 (MQ=255)                                                                                               ccTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgt                                                                                                                                     <  1:369729/251‑1 (MQ=255)                                                                                                  cacaGGTGTGGCAGCCGATACATAGCGTGGAGTCAGCAATTACAAAACGATGCACCAGGCATTCCGCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTGCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCGTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGgtg                                                                                                                                    <  1:141285/249‑1 (MQ=255)                                                                                                      aGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCCTTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGtata                                                                                                                                >  1:242135/1‑249 (MQ=255)                                                                                                             cAGCCGCTACAGAGCGTGGAGTCAGCAATTACACCACGATTCACCAGGCATTCATCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCACAAACAGACGTAACACCTGCAACACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATCGCGGTGTATAAGCTGTCGt                                                                                                                       >  2:436849/1‑251 (MQ=255)                                                                                                                ccGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGCTGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTaa                                                                                                                     >  1:71032/1‑250 (MQ=255)                                                                                                                 cgaaccAGAGCGTGGGGTCAGAAATTCAAAAACGATGTACCAGGCATTCCGCAGGTTATTGTCAATTTGGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCGCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAAc                                                                                                                   <  2:250440/246‑1 (MQ=255)                                                                                                                   aTACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGc                                                                                                                 <  1:283510/251‑1 (MQ=255)                                                                                                                           ggTGGAGTCTGCAATTACAAAACAATCCCCCAGGCATTCCTCAGGGTCTTGGAATTTTTGCCGAAACAATGCCGATGAAATGTCATTTTCGACACGCATCGACACGCCCATCCCCAAACAGACGTATCGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTGTACCAGTTGGCTTAAATGTACCAGCACATTGGTTTCAATAGCGGTGTATAAGCTGTCGTAATCGCCAGAGAtt                                                                                                         <  2:314188/250‑1 (MQ=255)                                                                                                                            gTGGAGTCAGCCATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAACCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCATGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCCCATTGTTTTCAATAGCGGTGTATAAGCTGCCGTAAACGCCAGCGAttt                                                                                                        >  2:272752/1‑251 (MQ=255)                                                                                                                            gTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAttt                                                                                                        >  1:83638/1‑251 (MQ=255)                                                                                                                                              aaaCGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTcc                                                                                       >  1:137819/1‑250 (MQ=255)                                                                                                                                                     tCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTg                                                                               >  2:44068/1‑251 (MQ=255)                                                                                                                                                      cgccAGGCATTCCTCGTTTTATTTTCATTTTGGACGAAAACATGCAGATGAAATTTCATTTTCGAGACCCATCGACACGCCCATCCCAAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGTTGTATAAGCGGTCGTAATCGCCAGCGATTTTTTCGAGATAGTGTTCCAGTACCTg                                                                               <  2:235075/248‑1 (MQ=255)                                                                                                                                                      ccccAGGCATCCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTg                                                                               <  2:101839/248‑1 (MQ=255)                                                                                                                                                           ggCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAGCGCCAGCGATTTTTGCGAGATAGTGTTCCAGAACCTGTTcgcg                                                                         <  1:418839/251‑1 (MQ=255)                                                                                                                                                                    cAGGTGATTGTCTTTTGTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCTCGATCGCGGa                                                                 <  2:385577/250‑1 (MQ=255)                                                                                                                                                                    cAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCgg                                                                  <  1:80376/249‑1 (MQ=255)                                                                                                                                                                                            aaCATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGcc                                        >  1:331711/1‑251 (MQ=255)                                                                                                                                                                                                   cGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATg                                 <  1:93185/251‑1 (MQ=255)                                                                                                                                                                                                      tGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAAAGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCAGTCGTAAACGCAAGCGATTTTTTCGAGATAGTGTTCCAGCACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCCATCAGGCATTCAATGcac                              >  2:381035/1‑249 (MQ=255)                                                                                                                                                                                                         aaTGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGt                           <  1:168591/251‑1 (MQ=255)                                                                                                                                                                                                            gTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTgc                        >  1:204509/1‑251 (MQ=255)                                                                                                                                                                                                               aTTTTCGACAATCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCACTCATGCCGTCAATGCTCAGTTGCGcc                     >  1:95381/1‑251 (MQ=255)                                                                                                                                                                                                                        acccaacGACACGCCCATCCCCAAACATACGTAACGCCTGCAAAACGGGCAAAGCCGCAGCTCATGCTTCCGGGCTTAGTCCCTTTACGAGTTGGCTTAATTTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATc            <  2:361459/245‑1 (MQ=255)                                                                                                                                                                                                                        acTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATc            <  2:361460/251‑1 (MQ=255)                                                                                                                                                                                                                           cATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGc         >  1:238712/1‑251 (MQ=255)                                                                                                                                                                                                                            aTCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCt        >  1:15355/1‑251 (MQ=255)                                                                                                                                                                                                                               gACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGtt     >  1:102927/1‑251 (MQ=255)                                                                                                                                                                                                                                 cacGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAGAGCCTCAGCTCATGCTGCCGGGCTTTGTCACTTTACCAGTTCGCTTAAATTCACAGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCCGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGTCCTTGCCATCAATCATGCCGTCACTGATCAGTTGCGCCTGCGCAATCCGCTGTTTg   >  1:457767/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGTGCAAAGCCTCAGCTCATGCTGCCGGGCGTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCTCATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCTCGATCGCGGATGTTGACCTTGGCAGCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGt  <  1:436849/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTgcg                       <  1:98361/228‑1 (MQ=255)                                                                                                                                                                                                                                    gCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTgcg                       >  2:98361/1‑228 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                         CAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGT  >  NC_000913/2849698‑2850173

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑