Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,962,881	A→G	L31P (CTA→CCA)	flhE ←	flagellar protein FlhE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,962,881	0	A	G	100.0%	246.9 / NA	70	L31P (CTA→CCA)	flhE	flagellar protein FlhE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/33);  total (37/33)

CACCGCCTGGCACTTCCCAGATAAATCGCAACGGTTCTGCTGCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTAGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCC  >  NC_000913/1962639‑1963125                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cACCGCCTGGCACTTCCCAGATAAATCGCAACGGTTCTGCTGCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAAtt                                                                                                                                                                                                                                              >  2:45464/1‑251 (MQ=255)         ggCACTTCCCAGATAAATCGCAACGGTTCTGCTGCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTaca                                                                                                                                                                                                                                        <  1:214402/249‑1 (MQ=255)               tCCCAGATAAATCGCAACGGTTCTGCTGCCGCTATGCCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTGCTGAAGGTGTAGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACGACtt                                                                                                                                                                                                                                  <  2:273891/249‑1 (MQ=255)                              aaCGGTTCTGCTGCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGCTAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGACAGAGGCGCAGACGACATCGACTCCCCGCGATGCTTTGGCGCAATTCCTACACTACTTGCCTGCCACCtcccctc                                                                                                                                                                                                                 >  2:96734/1‑251 (MQ=255)                                 ggtgctgctgCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACAtcccctcccc                                                                                                                                                                                                              <  2:264771/247‑1 (MQ=255)                                          gcCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACCCAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGTAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCtt                                                                                                                                                                                                      <  2:284410/250‑1 (MQ=255)                                          gcCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTg                                                                                                                                                                                                     >  2:212048/1‑251 (MQ=255)                                          gcCGCTATACCAGAAAAGGCCACGGTGATTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGGAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACAGCCCCTCCCCGGCGGCt                                                                                                                                                                                                       <  2:265111/249‑1 (MQ=255)                                            cGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTCATTCGACACCACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATACTGAAGCAGGTTGGCGCCGAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTCGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCCGCTTGc                                                                                                                                                                                                    >  1:174403/1‑250 (MQ=255)                                            cGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGc                                                                                                                                                                                                    >  2:304665/1‑250 (MQ=255)                                            cGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTAGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCa                                                                                                                                                                                                   <  1:187965/251‑1 (MQ=255)                                             gCTATAGCGGAAAAGGCCACGGTGGTTCCCTTCTGCACCTGTAATTCGACACAACGAGATTGCGAACACCAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTTTCGAGTAGAAAGAGGCGAAGACGACATCGACTCACCGCGATGAGTTGGCGCAATTCCTACACTCCTTGCCTGCTACATCCCCTCCCCGGCGTCTTGCAc                                                                                                                                                                                                  >  2:74205/1‑251 (MQ=255)                                              cTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCAcc                                                                                                                                                                                                 >  1:110448/1‑251 (MQ=255)                                                    cGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCagcagc                                                                                                                                                                                           <  1:62167/251‑1 (MQ=255)                                                       aaaaGGCGTCGGTGGTTCCGCGCTGCCCCTCTAAGTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACGCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGa                                                                                                                                                                                        <  1:215029/251‑1 (MQ=255)                                                            gCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACa                                                                                                                                                                                    >  1:111892/1‑250 (MQ=255)                                                              cACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCcac                                                                                                                                                                                                                         >  1:139965/1‑211 (MQ=255)                                                              cACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCcac                                                                                                                                                                                                                         <  2:139965/211‑1 (MQ=255)                                                                  gtggtTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTaa                                                                                                                                                                                                                                                >  2:400819/1‑184 (MQ=255)                                                                  gtggtTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTaa                                                                                                                                                                                                                                                <  1:400819/184‑1 (MQ=255)                                                                    ggtTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACCCAAGCGAACCCGCATTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTCCCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATAtt                                                                                                                                                                            <  2:82617/250‑1 (MQ=255)                                                                    ggtTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTg                                                                                                                                                                           <  1:139267/251‑1 (MQ=255)                                                                         cGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCtaat                                                                                                                                                                      >  1:6130/1‑251 (MQ=255)                                                                                                   cGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAAtt                                                                                                                                            <  1:22193/251‑1 (MQ=255)                                                                                                    gagaTTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTAGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTg                                                                                                                                           >  1:168744/1‑251 (MQ=255)                                                                                                      gaTTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTc                                                                                                                                         >  2:161813/1‑251 (MQ=255)                                                                                                      gaTTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTGGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTc                                                                                                                                         <  1:224931/251‑1 (MQ=255)                                                                                                            gAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTc                                                                                                                                   >  1:378114/1‑251 (MQ=255)                                                                                                              acacAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCAt                                                                                                                                 >  2:184491/1‑251 (MQ=255)                                                                                                                   aGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCgg                                                                                                                             >  1:161830/1‑250 (MQ=255)                                                                                                                      gAACCCGGAGTCCTGCAGGCGTCGGGACGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAAACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCCGCAGCGGAAACAATAATCTTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGAatat                                                                                                                         >  2:181833/1‑251 (MQ=255)                                                                                                                       aaCCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATg                                                                                                                        >  1:129033/1‑251 (MQ=255)                                                                                                                             cAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGtt                                                                                                                  <  1:45464/251‑1 (MQ=255)                                                                                                                                 ccTGAAGGTGGCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCGGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATTTCGGTTATCa                                                                                                              <  2:206314/251‑1 (MQ=255)                                                                                                                                      aGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTGGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGAc                                                                                                           >  1:22755/1‑249 (MQ=255)                                                                                                                                         tgtCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTc                                                                                                      >  2:252525/1‑251 (MQ=255)                                                                                                                                            cGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCACTCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATCTGTCGGTTCTCAAACAGTTcc                                                                                                     >  2:371341/1‑249 (MQ=255)                                                                                                                                                ccGATAAGCTGATAACGCCTCGCTACCAGCGTCATCAAGCCTGAAGCAGGTTGGCGCGTAGAAAGCGGCGCAGCCGCCATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTCCTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTc                                                                                               <  2:249820/251‑1 (MQ=255)                                                                                                                                                  gATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATACTATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTGTCAGACAGTTcccgcttcg                                                                                              >  1:242700/1‑243 (MQ=255)                                                                                                                                                   aTAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTACTAAGGTTCTCATTATTTGCCGCAAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGAAAGTTCCAGATTCGAc                                                                                            >  1:188334/1‑251 (MQ=255)                                                                                                                                                   aTAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGAc                                                                                            <  1:237293/251‑1 (MQ=255)                                                                                                                                                    tAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGAcc                                                                                           >  2:185634/1‑250 (MQ=255)                                                                                                                                                             aaCGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGTCAGTTACAGATTTGACAGGACCACt                                                                                   >  2:154775/1‑250 (MQ=255)                                                                                                                                                             aaCGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGGTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACt                                                                                   <  1:73682/250‑1 (MQ=255)                                                                                                                                                               cGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATAtt                                                                                                                                                                            <  1:228131/159‑1 (MQ=255)                                                                                                                                                               cGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATAtt                                                                                                                                                                            >  2:228131/1‑159 (MQ=255)                                                                                                                                                                     gCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGAAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGc                                                                                                                                                                          >  2:74179/1‑155 (MQ=255)                                                                                                                                                                       tACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCa                                                                         >  1:81531/1‑250 (MQ=255)                                                                                                                                                                          cAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACGCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGAGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCt                                                                     <  2:110448/251‑1 (MQ=255)                                                                                                                                                                              gTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTgcgg                                                                 <  1:444644/251‑1 (MQ=255)                                                                                                                                                                                   caATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGgcgc                                                              >  1:384638/1‑249 (MQ=255)                                                                                                                                                                                    aaTCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCa                                                                                                                                                            <  1:372682/154‑1 (MQ=255)                                                                                                                                                                                    aaTCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCa                                                                                                                                                            >  2:372682/1‑154 (MQ=255)                                                                                                                                                                                    aaTCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAgg                                                           <  1:236548/251‑1 (MQ=255)                                                                                                                                                                                      tCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTAATTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGa                                                          >  1:297863/1‑250 (MQ=255)                                                                                                                                                                                          gattcAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCAATGACTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTTCCCCAGCAGCGGAAACAATAATATTGCTAATATGGTTCTCGTTATTTGCCGCCAATTGTTGCCGTCAGGCGGATATGTCGGTTATCAGACATTTCCAGATTCTACAGGACCACTAACTGCGGCATGCTGCGGCGCAGGAAGc                                                       <  1:405730/245‑1 (MQ=255)                                                                                                                                                                                           aaGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTGGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCgag                                                    <  2:187922/251‑1 (MQ=255)                                                                                                                                                                                                ggTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACaat                                               >  1:454632/1‑251 (MQ=255)                                                                                                                                                                                                       gcgTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTc                                        <  1:239489/251‑1 (MQ=255)                                                                                                                                                                                                         gTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCCGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGc                                      >  2:208282/1‑251 (MQ=255)                                                                                                                                                                                                             aaaGAGGCGCAGACGACATCGACTCACCGCGATGATCTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCCGGAAGCGAGCCAATAATGGTCGCAgcg                                  >  2:364813/1‑251 (MQ=255)                                                                                                                                                                                                               agagGCGCAGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGcgcg                                <  1:454921/251‑1 (MQ=255)                                                                                                                                                                                                                       aGACGACATCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGCCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCAc                        <  2:426639/251‑1 (MQ=255)                                                                                                                                                                                                                              aTCGACTCACCGCGATGATTTGGCGTAATTCCTCCACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAAt                 <  2:341186/251‑1 (MQ=255)                                                                                                                                                                                                                              aTCGACTCACCGCGATGATTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAAt                 <  1:208282/251‑1 (MQ=255)                                                                                                                                                                                                                                   cTCACCGCGATGAGTTGGCGAAATTCCTACACTACGTGCCGGCCACATCCACTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGATCCAGAATTGACAGGACAACTAACTGAGGCAAGCTGCGGCACGGGACGCGAGACAAAAATGGTCGCAGCGCATGGATAACCAAAAATCCTgg            >  2:208018/1‑251 (MQ=255)                                                                                                                                                                                                                                           gatgatTTGGCTTAAGTCCTACACTACTTGCCGGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCAc    <  2:242700/251‑1 (MQ=255)                                                                                                                                                                                                                                           gatgatTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTggcggcc       >  2:62833/1‑247 (MQ=255)                                                                                                                                                                                                                                           gatgatTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGgcggcg       <  1:62833/248‑1 (MQ=255)                                                                                                                                                                                                                                             tgatTTGGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCAccc  >  1:40732/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CACCGCCTGGCACTTCCCAGATAAATCGCAACGGTTCTGCTGCCGCTATACCGGAAAAGGCCACGGTGGTTCCGCTCTGCCCCTCTAATTCGACACAACGAGATTGCGAACACAAGCGAACCCGCAGTCCTGAAGGTGTCGGGCCGATAAGCTGATAACGCCACGCTACCAGCGTCATCAATCCTGAAGCAGGTTGTCGCGTAGAAAGAGGCGCAGACGACATCGACTCACCGCGATGATTTAGCGTAATTCCTACACTACTTGCCTGCCACATCCCCTCCCCGGCGGCTTGCACCAGCAGCGGAAACAATAATATTGCTAATAAGGTTCTCATTATTTGCCGCCAATTGTCGCCGTCATGCGGATATGTCGGTTATCAGACAGTTCCAGATTCGACAGGACCACTAACTGCGGCAAGCTGCGGCGCAGGAAGCGAGACAATAATGGTCGCAGCGCGTGGTTCACCAACAATACTGGCGGCGCACCC  >  NC_000913/1962639‑1963125

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