Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,089	C→A	intergenic (+57/‑545)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,089	0	C	A	98.6%	207.0 / NA	70	intergenic (+57/‑545)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (21/48);  total (22/48)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACA‑TCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTT  >  NC_000913/88845‑89329                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ttGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAATTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGAc                                                                                                                                                                                                                                              <  1:944618/250‑1 (MQ=255) ttGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGAc                                                                                                                                                                                                                                              <  1:599650/250‑1 (MQ=255)     cGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCAGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATAACCGCCAGGCGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACA‑TCac                                                                                                                                                                                                                                        >  2:945380/1‑250 (MQ=255)      gATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAat                                                                                                                                                                                                                                       <  1:65234/251‑1 (MQ=255)           cGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGTTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAatatcct                                                                                                                                                                                                                                  <  1:841440/251‑4 (MQ=255)           cGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAatatcct                                                                                                                                                                                                                                  <  1:130050/251‑4 (MQ=255)                   tCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGG‑CACTCAATATCCTTAAAGTa                                                                                                                                                                                                                           >  2:974290/1‑250 (MQ=255)                    cGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTaaa                                                                                                                                                                                                                         <  1:915117/251‑1 (MQ=255)                               aGTGGCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTc                                                                                                                                                                                                              <  1:28109/251‑1 (MQ=255)                                 tgtCCGGTGCTGGCCGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCt                                                                                                                                                                                                            >  1:224965/1‑251 (MQ=255)                                  gtCCGGTGCTGGCAGTGGCTCAACGTCACCGCGGTGTCGAAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGCCGAACCGCGTATGCCAAAACCTG‑GGTTGACGCGCATTACACGTAATCTCTATCGCCGCGCCGTGCTCAGCCGTAGCTAAGCCGCGACCAAAAATAAGCGCCCGGTGAATTTCACTCTGGCGCGTAGAGTCCGGGAATGGACA‑TAAATCTCCTTAACGTAAATCCGACTATTCCTg                                                                                                                                                                                                           >  2:583100/1‑251 (MQ=255)                                       gTGCTGGCAGTGGATCAACGTCAACGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATAGCCGCGGCGTGCTCAGCCGTAGCTCAGCCGCGAACAAAAATCAGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTgccact                                                                                                                                                                                                      >  2:993345/1‑246 (MQ=255)                                        ggCTGGCAGTGGCTCAACGGCACCGCGAGGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACtatt                                                                                                                                                                                                     <  2:9543/250‑1 (MQ=255)                                         gCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACtatta                                                                                                                                                                                                    <  1:868274/251‑1 (MQ=255)                                         gCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACtatt                                                                                                                                                                                                     <  1:774629/250‑1 (MQ=255)                                           tGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTattatt                                                                                                                                                                                                  >  1:654191/1‑251 (MQ=255)                                           tGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTattatt                                                                                                                                                                                                  >  1:135887/1‑251 (MQ=255)                                            ggCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTattatt                                                                                                                                                                                                  <  1:878975/250‑1 (MQ=255)                                                gTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaa                                                                                                                                                                                             <  1:133663/251‑1 (MQ=255)                                                      cAACGTCCCCGCGATGTCTCAGATCGTGTGCTGGAGATTGTTCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCTCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATaagttctt                                                                                                                                                                                       <  1:318562/251‑7 (MQ=255)                                                            cACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAAAAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAAc                                                                                                                                                                                 >  2:758860/1‑251 (MQ=255)                                                            cACCGCGATGTCGCAGAGCGAGTGCTGGAGATTGTCCTGGCAAGCCTTGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTaaa                                                                                                                                                                                  <  1:790076/250‑1 (MQ=255)                                                                       cGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGtta                                                                                                                                                                       <  1:717748/250‑1 (MQ=255)                                                                           gagctggtgCAGGAGAGTGTCATGGCAATCCTGGACGAACCTCGTAAGCCAAAACTTGTGTTTG‑CGCGCATTAAACGTAATCTCTATCGCCGTGACGGGCTCATCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCTTCTGGCGCGTATAGTACGGTAATGGACA‑TCAATATCCTTAAAGTAAATAATGCTAGACTTGACTATTGTTGATAAGTTCTTTTAACcccgcactggaattaa                                                                                                                                                                  <  2:595927/245‑12 (MQ=255)                                                                              cGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTc                                                                                                                                                               <  1:718133/251‑1 (MQ=255)                                                                              cGTGGGCTGGAGTTTGTCATGGCAAGACTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTc                                                                                                                                                               <  2:2343/251‑1 (MQ=255)                                                                                      ggAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATAGCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTg                                                                                                                                                       <  1:628740/251‑1 (MQ=255)                                                                                       gagaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCACAACCTG‑ATTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAACTACGCGCAAGGTGAATTTCCATCTGGCGCGTAGAGTACAGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTg                                                                                                                                                       >  2:390998/1‑250 (MQ=255)                                                                                        agaTTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCAGTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                                                                                                                                     <  2:409812/251‑1 (MQ=255)                                                                                        agaTTGTCCTGGAAAGCCTGGTCGAACCGCGTAAGCAAAAACCTG‑TTTTAACGCGCATTAAACGTAATCTCTAGCGCCGCGGCGTGCTCAGCCGTAGCTAAGCAGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTCCGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTAATAAGTTATTTTAAACCCGCCCGTTAATTAACTCACCAGCTGaa                                                                                                                                                     <  2:584646/251‑1 (MQ=255)                                                                                            tGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAAtt                                                                                                                                                  <  1:739778/250‑1 (MQ=255)                                                                                             gTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTc                                                                                                                                                 <  1:893815/250‑1 (MQ=255)                                                                                               ccTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTcac                                                                                                                                               <  1:385175/250‑1 (MQ=255)                                                                                                  gccAAGCCTGGACGAACCGCGTAAGCCAAATCCTG‑GTGTAACGCGCAATAAACGTAATCTCTACCGCCGCGGCGTGCTAGGCCGTAGCTAAGACGCGAACAAAAATACGCGCCAGGGGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACAATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTTAAATTCACaa                                                                                                                                             <  2:768161/247‑1 (MQ=255)                                                                                                    cAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTACAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACaataat                                                                                                                                         >  2:587110/1‑251 (MQ=255)                                                                                                      aGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTAACGCCGCGGCGTGCGCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGGAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACACACCAGCTGAAATTCACAATAATTa                                                                                                                                       <  1:793413/251‑1 (MQ=255)                                                                                                      aGCCTGGAACAACCGCGTAAGTCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTa                                                                                                                                       <  1:525411/251‑1 (MQ=255)                                                                                                        ccTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAACACCGCCCGTTAATTAACTCACCAACTGAAATTCACAATAATTa                                                                                                                                       >  2:805962/1‑249 (MQ=255)                                                                                                              cGAACCGCGTGAGCCAAAACCTG‑GTTTAACGCGCATTAAACGAAATCTATATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCACTCTGGCGCGTAGAGGACGGGAATGGACA‑TACATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGCTAAGTTCTTTTAAACCCGCCCGTTATTTAACTCACCAGCTGAAATTCACAATAATTAAGTGatat                                                                                                                               >  2:612447/1‑251 (MQ=255)                                                                                                                  ccGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc                                                                                                                           <  1:120147/251‑1 (MQ=255)                                                                                                                  ccGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAc                                                                                                                           <  1:382104/251‑1 (MQ=255)                                                                                                                          gCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGAGATCGACAGCGCGtt                                                                                                                   <  2:629567/251‑1 (MQ=255)                                                                                                                              aaaCCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGACAGGTGAATTTTCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTg                                                                                                               <  1:382646/251‑1 (MQ=255)                                                                                                                                           aCGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGtt                                                                                                   <  1:811498/251‑1 (MQ=255)                                                                                                                                              cgcATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGCGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGAATATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGCTATCGACAGCGCGTTTTTGCATTATTTTGTTAc                                                                                                 >  1:557415/1‑250 (MQ=255)                                                                                                                                                    aaaCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATgc                                                                                             >  2:458376/1‑248 (MQ=255)                                                                                                                                                    aaaCGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGAAAGCGCGTTTTTGCATTATTTTGTTACATgc                                                                                             >  1:42745/1‑248 (MQ=255)                                                                                                                                                          aaTCTCTCTCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGa                                                                                     <  2:145928/250‑1 (MQ=255)                                                                                                                                                               ctATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAATAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTAAATGCAGAGATGAATTGcc                                                                               >  2:753490/1‑251 (MQ=255)                                                                                                                                                                  acgccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAAGACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTGGTTACATGCAGAGATGAATTGCCGAt                                                                            <  2:138653/250‑1 (MQ=255)                                                                                                                                                                   cgccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGAt                                                                            <  1:291423/250‑1 (MQ=255)                                                                                                                                                                   cgccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGAt                                                                            <  1:987486/250‑1 (MQ=255)                                                                                                                                                                    gccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGCCTATTATTGCTAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGAttt                                                                          >  2:12922/1‑251 (MQ=255)                                                                                                                                                                     ccgcGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGAGAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTa                                                                         <  1:606104/251‑1 (MQ=255)                                                                                                                                                                       gcggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaac                                                                       >  1:229572/1‑251 (MQ=255)                                                                                                                                                                        cggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTGTGTTACATGCAGAGATGAATTGCCGATTTaaca                                                                      <  1:285426/251‑1 (MQ=255)                                                                                                                                                                          gcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaacaag                                                                    <  1:458135/251‑2 (MQ=255)                                                                                                                                                                               cTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCCGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTCAATAAGACTATTACTGACTATTCTTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAACTTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGAAGAGATGAATTGCCGATTTaacacgaactt                                                               >  2:738021/1‑244 (MQ=255)                                                                                                                                                                                 cAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACtttt                                                             <  1:607162/251‑1 (MQ=255)                                                                                                                                                                                   gCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGGGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCt                                                           <  2:452821/251‑1 (MQ=255)                                                                                                                                                                                     cGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCttt                                                         <  1:221974/251‑1 (MQ=255)                                                                                                                                                                                         acTAAGCTGCGCACAAAAATACACGCCAGTTGCATTTCCCTCTGGCGCGTAGTGTACGGGAATGAACA‑TCACTATCCTTAAAGTAAATAAGACTTTTCCTGACTATTATTGATAAGTTCTTTTAAACCAGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACGGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGACGATTAAACAAGAACTTTTCTTTGCtt                                                     >  2:736325/2‑251 (MQ=255)                                                                                                                                                                                            aaGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCtttt                                                   <  1:756358/250‑1 (MQ=255)                                                                                                                                                                                            aaGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCttt                                                    <  2:135887/249‑1 (MQ=255)                                                                                                                                                                                            aaGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCGTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATAGCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTg                                                  <  1:890448/251‑1 (MQ=255)                                                                                                                                                                                                          aaTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTgg                                    <  1:3046/251‑1 (MQ=255)                                                                                                                                                                                                                            ttCCCTCTGGCGCGTAGAGTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCTACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCATACCCGCTGGCATCAAGCGCCGCACAg                   <  1:753490/250‑1 (MQ=255)                                                                                                                                                                                                                              cccTCTGGCGCGTAGAGTACGGGAATGGACA‑TAAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAGCCGCTGGCATCAAGCGCCACACAGACg                >  2:680917/1‑251 (MQ=255)                                                                                                                                                                                                                                 tctGGCGCGTAGAGTACGGGAATGGAAA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACCCAGACgca              >  2:905929/1‑248 (MQ=255)                                                                                                                                                                                                                                            gagTACGGGAATGGACA‑TCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGtt  <  1:986919/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TTGGCGATAACGAACTGCTCGACTTCTTACAGTGTCCGGTGCTGGCAGTGGCTCAACGTCACCGCGATGTCGCAGAGCGTGTGCTGGAGATTGTCCTGGCAAGCCTGGACGAACCGCGTAAGCCAAAACCTG‑GTTTAACGCGCATTAAACGTAATCTCTATCGCCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACA‑TCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTT  >  NC_000913/88845‑89329

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