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	95.0%	101.7 / NA	41	A34V (GCG→GTG)	cheA	chemotaxis protein CheA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/20);  minor base T (0/2);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACC  >  NC_000913/1974983‑1975467                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ggcggcATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGc                                                                                                                                                                                                                                             <  1:170078/250‑1 (MQ=255)         ccGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAAcc                                                                                                                                                                                                                                    >  1:424158/1‑251 (MQ=255)                            ataAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTTATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATg                                                                                                                                                                                                                  <  2:424158/250‑1 (MQ=255)                            ataAGCGCCGTGCCGTTATTGCATGATGGCCTTCGTTTCCAAAAACAGATTGATAAGGTCGGTGTGGAGTTGCATCGCTCCTCGGCTGGCGTCATCGAGCAGGTTGTCCATCAGATGCGTGGTTTCCTTCAAAACGCGGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATGCAATTGGTCGGCATCTGGCACTTCCGGCTGCAAGCCCAGCAAAtgtgtcgccatg                                                                                                                                                                                                                  <  2:22173/250‑11 (MQ=255)                                 cGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGcc                                                                                                                                                                                                            >  1:233870/1‑251 (MQ=255)                                  gTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCCGcc                                                                                                                                                                                                            >  2:352800/1‑250 (MQ=255)                                     gAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCCCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAACTTCCTGCCCCTCCTTTGATCGCGTGGGAAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCCCTTCCGGCTGCAACACCAGCAAATGCTGCTCCATGTCAGCCAAc                                                                                                                                                                                                         >  1:60675/1‑250 (MQ=255)                                        cggTTCTTGCATGATGTCCTTCGTTGCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCAGCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTTATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCAGGTCAGCCAGCAGtt                                                                                                                                                                                                     <  2:428598/249‑1 (MQ=255)                                               tGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTGGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTcc                                                                                                                                                                                                <  2:133717/249‑1 (MQ=255)                                                    gatttctttcttttCCAAAAACATATTGATAATGTCGTTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTGCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTGCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTtcat                                                                                                                                                                                         <  2:11048/240‑1 (MQ=255)                                                         ccTTCTTTTCCAAAAACAGATTGATAATGTCGGTGCTTAGTTGCATCGTACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGGCAGCCAACAGTTCGTCCGCTTCATCaaa                                                                                                                                                                                     <  2:242652/250‑1 (MQ=255)                                                               tttCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTc                                                                                                                                                                              >  1:293654/1‑251 (MQ=255)                                                                         cAGATTGATAATGTCGGTGTTGAGTTGCATCTCTCCTCGTCTGTCTTCAGCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAAt                                                                                                                                                                    <  2:337310/251‑1 (MQ=255)                                                                                  aaggtctgtgtTGATTTGCATCTCATCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCATCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTatat                                                                                                                                                           <  2:172318/244‑1 (MQ=255)                                                                                           gtTGAGTTGCATCCCCCCTCGTCTGGCTTCATCGAGCAGGTTTACCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGCTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCCTGCTCa                                                                                                                                                  <  1:163645/251‑1 (MQ=255)                                                                                              gAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGc                                                                                                                                               >  1:90406/1‑251 (MQ=255)                                                                                                gTTGCATTTCACCGCGGCGGGCTTCATCGAGCAGGGTTTCCATCAGATGCGTTGTTTCCTGCAAAACGATGAAGCAAAAAGTTCCTGCCCCTCTTTTGAACGAGTGGGCAGCCCGAAAGATGGCAGTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCGAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAAGCGCTTATAGCCAGTCTCACGCTg                                                                                                                                             <  2:60675/251‑1 (MQ=255)                                                                                                     aTCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGATCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCAcc                                                                                                                                        >  1:39831/1‑251 (MQ=255)                                                                                                        tcACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCAGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATCTCCATGCTCACGCTGTCACCTCg                                                                                                                                     >  1:150499/1‑251 (MQ=255)                                                                                                          cccTCTTCTGGCTTCATCGAGCAGGTTTTCCATCAGAGGCGTGGTTTCCTGTAAAACGCTGAAGCCAAAAGTGCCTTCCCCTCCTTTGATCGAGTGGGCAGCCCGAAATATGGCATTCAATTGTTCGGCATCTGGCTCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTTTCACCTCGGt                                                                                                                                   <  2:309094/250‑1 (MQ=255)                                                                                                                                   ttttCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAacac                                                                                                          >  2:40805/1‑251 (MQ=255)                                                                                                                                        caacAGATGCGTGGTTTCCTGCAAAACGCTGAAGCAAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGAcc                                                                                                     <  2:293654/248‑1 (MQ=255)                                                                                                                                                    ggtgtacggaaaaaCGCTGAAGCCAAAAGTGCCGGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTGCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAAGCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAAcc                                                                                         <  2:90406/241‑1 (MQ=255)                                                                                                                                                        tCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTc                                                                                       >  1:218793/1‑249 (MQ=255)                                                                                                                                                                     tGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGTGCAGCCCGAAAGATGGTATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAACAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCCTTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAggcc                                                                        >  2:7451/1‑249 (MQ=255)                                                                                                                                                                          ccAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCtt                                                                    >  1:363430/1‑250 (MQ=255)                                                                                                                                                                             aaaGTTCCTGCCCCTCCGTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTg                                                                <  1:332973/251‑1 (MQ=255)                                                                                                                                                                              aaGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTgg                                                               >  1:47857/1‑251 (MQ=255)                                                                                                                                                                                gTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGACCACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCt                                                             >  2:125782/1‑251 (MQ=255)                                                                                                                                                                                         cctcctTTGATCGCGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGTGCGCTTCCTGGCTACTTCTGGc                                                    >  2:19701/1‑251 (MQ=255)                                                                                                                                                                                               ttGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTc                                                <  1:191966/249‑1 (MQ=255)                                                                                                                                                                                                           gCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCa                                  <  1:107762/251‑1 (MQ=255)                                                                                                                                                                                                               cagacttgtgctcctccgatcattCGGCCTCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAAt                              <  1:117996/229‑1 (MQ=255)                                                                                                                                                                                                                 cGAAAGATGGCATTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGATCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCCGGGCACTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGAAAAATgg                            >  2:359647/1‑251 (MQ=255)                                                                                                                                                                                                                       aTGGCATTCAATTGGTCGTCAACTGGCACCTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAGCAGTGCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCGTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACGGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCGCATTCTGGCTTTCGGCGTTTTCCTGCAAAATGGCCTGtt                      <  2:39831/251‑1 (MQ=255)                                                                                                                                                                                                                         ggCATGCAATTGTTCGGCATCTGGCACTTCCGGCTGAAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCGCCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTc                     <  2:83910/250‑1 (MQ=255)                                                                                                                                                                                                                            aTTCAATTGTTCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCt                 >  1:305538/1‑251 (MQ=255)                                                                                                                                                                                                                                aaTTGTGCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCttgtt             <  2:179001/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTTGTTCGGCATCTGGCTCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCCTGTCATCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCGCCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCACTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCttgttt            <  2:420232/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttCGGCATCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAAt                              >  2:117996/1‑229 (MQ=255)                                                                                                                                                                                                                                            tCTGGCACTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGCTAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAACTGGCCTGTTCGGCTTGTTTGTTCAGTaac  >  1:314867/1‑248 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    GGCGGCATCCGGCTCTTGCGACTGTTTATAAGCGTCGAGCTGTTCTTGCATGATGTCCTTCGTTTCCAAAAACAGATTGATAATGTCGGTGTTGAGTTGCATCTCACCTCGTCTGGCTTCATCGAGCAGGTTTTCCATCAGATGCGTGGTTTCCTGCAAAACGCTGAAGCCAAAAGTTCCTGCCCCTCCTTTGATCGAGTGGGCAGCCCGAAAGATGGCATTCAATTGTTCGGCATCTGGCGCTTCCGGCTGCAAAACCAGCAAATGCTGCTCCATGTCAGCCAACAGTTCGTCCGCTTCATCAAAAAATGTCTGATAAAAATCGCTTATATCCATGCTCACGCTGTCACCTCGGTTCGGCTGATGGCATTGTGGGAACACTGACCTGTGGTGCAACCTCAGGTTTTTCCAGGGCGCTTACTGGCTCATTCTGGCTTTCGGCGTTTTCATGCAAAATGGCCTGTTCGGCTTGTTTGTTCAGTACC  >  NC_000913/1974983‑1975467

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