Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,975,224	G→A	A34V (GCG→GTG)	cheA ←	chemotaxis protein CheA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,975,224	0	G	A	100.0%	163.9 / NA	51	A34V (GCG→GTG)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/30);  total (21/30)

GAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAG  >  NC_000913/1974977‑1975463                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gtagcgggcggcAGCCGGCCCGTGCGACTGTTTAGAAGCGTAGAGCTGTTCTTGCATGATGTCCTTGGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCACCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACtt                                                                                                                                                                                                                                              <  2:218745/245‑1 (MQ=255)            cATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCaaa                                                                                                                                                                                                                                   <  1:430516/251‑1 (MQ=255)                cGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAAcc                                                                                                                                                                                                                                >  1:368330/1‑250 (MQ=255)                       tGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCatat                                                                                                                                                                                                                         <  2:560981/250‑4 (MQ=255)                        gCGACTGTTTATAAGCTTCGATCTGTTCTTGCATGATGTTCTTAGTGTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATTTCACCTCGTCTTGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGTCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGTACTTCCGGCTGCAAAACCAGCAAAtgc                                                                                                                                                                                                                       <  1:402877/251‑1 (MQ=255)                          gACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAAtgctg                                                                                                                                                                                                                     >  2:245913/1‑251 (MQ=255)                                           gAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAAc                                                                                                                                                                                                     >  1:237320/1‑250 (MQ=255)                                           gAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACa                                                                                                                                                                                                    <  2:69404/251‑1 (MQ=255)                                           gAGCTGTGCTTGTATGATGTCCTTCTTTTCCAAAAAAAGATGTATAAGGTCGGTGTTGAGTTGCAGCTCACCTCGTCTTGCTGCATCGAGCAGGTTTTGCAGCAGATGCGTGGTTTACTGCAAACCGCTGAATCCAAAAGGTCCTGCTCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGTATTCAATTGTTCGGCATCTGGCACTTCAGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCcaca                                                                                                                                                                                                    <  2:539642/251‑1 (MQ=255)                                               tGTTCTTGCATGATGTCCTTCGTTTACAAAAACAGCTTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCAATGTCAGCCAACCGTTc                                                                                                                                                                                                >  1:303522/1‑251 (MQ=255)                                                     tGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTGTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTcc                                                                                                                                                                                            <  2:436124/249‑1 (MQ=255)                                                      gCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCt                                                                                                                                                                                         <  1:27681/251‑1 (MQ=255)                                                        atgatgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTtc                                                                                                                                                                                       <  1:111405/251‑1 (MQ=255)                                                          gatgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGTTTTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGTGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTtcat                                                                                                                                                                                     <  2:108038/251‑1 (MQ=255)                                                           atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGTTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTtcatc                                                                                                                                                                                    <  2:302206/251‑1 (MQ=255)                                                           atgTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTtcatc                                                                                                                                                                                    <  1:10048/251‑1 (MQ=255)                                                              ccctgcttttCCAAAAACAGACTGAAAATGTCGGTGTTGAGTAGCATCTCACCACGTCTTGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTCCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCTCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAAGGTTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCaaa                                                                                                                                                                                 <  1:564875/244‑1 (MQ=255)                                                                cTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCaaaa                                                                                                                                                                                >  1:57729/1‑250 (MQ=255)                                                                      ttCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGACCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCACCAGTTCGTCCGCTTCATCAAAAAATGTCt                                                                                                                                                                         >  1:539147/1‑251 (MQ=255)                                                                       tCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCt                                                                                                                                                                         <  1:285262/250‑1 (MQ=255)                                                                          aaaaaCAGATTGATAATGACGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAATTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCGCCATGTCAGCCAACAGTTAGTCCGCTTCATCAAACAATGTCTGAt                                                                                                                                                                      >  1:244089/1‑250 (MQ=255)                                                                           aaaaCAGATTGTTAATGTCGGTTGTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATaa                                                                                                                                                                    <  2:116777/251‑1 (MQ=255)                                                                               cAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAACCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATaaaact                                                                                                                                                                >  1:43398/1‑249 (MQ=255)                                                                                    tGATAATGTCGGTGTGTAGTGTCATCTCACCCCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAATCGCTGAAGCCAAAAGTTCCTGCCCCTCCTTGGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCt                                                                                                                                                            <  2:268569/250‑1 (MQ=255)                                                                                              ggTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCACTTATATCCATGc                                                                                                                                                 >  1:552580/1‑251 (MQ=255)                                                                                                     aGTTGCATCTCACCTCGTCTTGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCATAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTTATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCActct                                                                                                                                          <  2:405509/251‑4 (MQ=255)                                                                                                       ttGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATACAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGt                                                                                                                                        >  2:493954/1‑251 (MQ=255)                                                                                                              tcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGCGTGGGCAGCCCGAACGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCACAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAACTGTCTGATAAACATCGCTTATATCAATGCTCACGCTCTCAccacg                                                                                                                                 >  2:245462/1‑248 (MQ=255)                                                                                                               cccctcgacgggctccaccGAGCAGGTTTTCCATCAGATGCGTGGTGTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCGTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGGTCGTTCGCTGCATCAAAAAATGTCTGATAAAAATCGCTTATAGCCAGGCTCACGCTGTCACCTc                                                                                                                                  <  2:74010/231‑1 (MQ=255)                                                                                                               caccccctCTGGCTTCACCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAACCGCTGAAGCCAAAAGTACCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAGTCGCTTATATCCATGCTCACGCTGTCACCTCg                                                                                                                                 <  2:447781/243‑1 (MQ=255)                                                                                                               cACCTCGTCTGGCTTCATCGCGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAATCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTc                                                                                                                                  <  1:325119/249‑1 (MQ=255)                                                                                                                  cTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCAATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTc                                                                                                                             >  1:477865/1‑251 (MQ=255)                                                                                                                         ggCTTCATCGGGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCGGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTg                                                                                                                        <  2:447285/249‑1 (MQ=255)                                                                                                                                   aGCAGGTTTTCCATCAGATGCGTGGTCTCCTGCAACCCGCTGCACCCAAACGATCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCAAACAGTTCGTCCGCTTCATCAAAAAGTGTCTGATAAAAATAGCTTATATCCATGCTCACGCTGTCACATCGGTTCGGCTGCTGCCCTTGTgg                                                                                                            >  1:227456/1‑251 (MQ=255)                                                                                                                                            tCCATCAGATGCGTGGTTTCCTGCAAACCGCTGAAGCCAAAAGTGCCTGCCCCTCCTTGGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATAGCCATGCTCACGCTGTCACCTCGGTTCGTCTGATGGCATTGTGGGAACACTg                                                                                                    <  2:323454/250‑1 (MQ=255)                                                                                                                                                    aTGCGTGGTTTACTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTgtggt                                                                                           <  2:465907/251‑1 (MQ=255)                                                                                                                                                       cGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTGGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCa                                                                                        <  2:237320/251‑1 (MQ=255)                                                                                                                                                          ggattCCTGCAAAACGCTGAAGCCAAATGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCAATTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATAAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTAGGAACCCTGACCTGTGGTGCAAcc                                                                                     >  1:416708/4‑251 (MQ=255)                                                                                                                                                               ccTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAg                                                                                 <  2:539147/250‑1 (MQ=255)                                                                                                                                                                   cAAAACGCTGAAGCCAAAGGTCCCTGCCCCTCCTTGTATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTGCGGCAGCTGGCACTTGCGGCTTCAAAACCAGCAAATGCTGCTCCATTTCAGCTAACAGGCCGTCCGCTGCATCAAAGAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCGCCTCAGTGCGGTTGATGGCATGGTGGGAAGTCTGGCCTGTGGTTCAACCTCATGtttt                                                                            <  2:217813/251‑1 (MQ=255)                                                                                                                                                                                        tCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTc                                                        <  2:70878/250‑1 (MQ=255)                                                                                                                                                                                             ccccTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGCCCTGTGGTGCAACCTCAGGTTTTTCCCGGACGCTTACTGGCTCATTCTg                                                  >  1:211974/1‑251 (MQ=255)                                                                                                                                                                                              cccTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTgg                                                 <  1:404816/251‑1 (MQ=255)                                                                                                                                                                                                    tttGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCg                                           >  1:316560/1‑251 (MQ=255)                                                                                                                                                                                                        aTCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGt                                       >  2:42767/1‑251 (MQ=255)                                                                                                                                                                                                          cGAGTGGGCAGCCCGAAAGATTGCATTCAATTGTTAGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCATCCGCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGtt                                      <  2:414241/250‑1 (MQ=255)                                                                                                                                                                                                            aGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTc                                   >  1:533425/1‑251 (MQ=255)                                                                                                                                                                                                                       cGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATgg                        >  1:181580/1‑251 (MQ=255)                                                                                                                                                                                                                             aTGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGtt                  >  2:304428/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCAATTGTTCGGCATCTGGAACTTCCGGCTGCACCACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGGGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTATGGCTTTCGGCGTCTTCATGCAAAATGGCCTTTTCGGCttg           >  2:384573/1‑251 (MQ=255)                                                                                                                                                                                                                                             cGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAg  <  1:457229/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GAAGCTGGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAG  >  NC_000913/1974977‑1975463

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