Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,069,065	G→A	20.2%	L259L (CTG→TTG)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,065	0	G	A	20.2%	123.6 / 4.4	57	L259L (CTG→TTG)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base G (37/7);  new base A (5/7);  total (42/15)
Fisher's exact test for biased strand distribution p-value = 5.87e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAG  >  NC_000913/1068838‑1069299                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            cAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGcc                                                                                                                                                                                                                                  >  1:727651/1‑238 (MQ=255)      tGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGGCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTCAGTTTTCCGCCCCCCAAA‑TCACCCCCCCggg                                                                                                                                                                                                                             >  2:137815/1‑238 (MQ=255)      tGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGGCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACggg                                                                                                                                                                                                                             >  1:636815/1‑238 (MQ=255)                   agcagcaACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCCACGAATGCCCGCCCCCTATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAAATGCCCTAGGTTTTCCCCCACCAAAA‑TCCCCGCCACGGGCGCacaaaacacaa                                                                                                                                                                                                               >  2:60373/1‑229 (MQ=255)                          acaaacaCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTg                                                                                                                                                                                                        >  1:415125/1‑239 (MQ=255)                                acCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCAGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATAACCTAAGCTTCACGCCACCCAAA‑TCACCCCCACGGGCCCAATCACCACAACTGCCCGACCact                                                                                                                                                                                                  >  2:166159/1‑237 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCacaaa                                                                                                                                                                                                >  2:717393/1‑235 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCACCATTGCCTGACCATTaa                                                                                                                                                                                                >  1:758741/1‑239 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAAAA‑TCACCGCCCCGGGCGCAATCAGCACCATTGCCTGACCATaaa                                                                                                                                                                                                >  2:20878/1‑239 (MQ=255)                                     cGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCACCCCACA‑TCACAGCCCCGGTCGCAATCAGCATCATTGCCGGACCATTaaaaa                                                                                                                                                                                             >  1:1015465/1‑238 (MQ=255)                                          aGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTCAGTTTTCCGCCCCCAGAA‑TCACCGCCCCGGGGGCAAGCAGACTAATTGCCGGACCATTaaaaaaaggg                                                                                                                                                                                        >  2:382254/1‑233 (MQ=255)                                          aGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGCCCCTATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCCCCACCAGAA‑TCACCGCCACGGGCGCAACCAGCATCATTGCCTGACCCTTAAAAACAggg                                                                                                                                                                                        >  2:1033884/1‑238 (MQ=255)                                                                     cGATATTTTCAGCATAGGTGGTCACGCCGCTACTGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCACCATCATTGCCTGCCCCTTAAAGACAGGAGGGGAAAAAGTGGGCAAACCGAACCAg                                                                                                                                                             >  1:466130/1‑239 (MQ=255)                                                                     cGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCGAACCGAACCAg                                                                                                                                                             >  1:420925/1‑239 (MQ=255)                                                                            ttCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCAt                                                                                                                                                      >  1:47546/1‑239 (MQ=255)                                                                                            aCGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAGACAGAAGTGGAAAAATGGGGCAGACCGAACCAGGCAGCATGGCGGACGAGAGTaaa                                                                                                                                      >  1:158322/1‑239 (MQ=255)                                                                                                               agCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCACCGGCTTTTCCCat                                                                                                                    >  1:712751/1‑237 (MQ=255)                                                                                                                               tGCCAGACCA‑CCCCTAAGGATGCCCCCCCCAAAAACCGGGCCAATTTGCCGCCCGGTCACCCCGGCTACTGCCTTGAGATGACCTAAGTTTCCCGCCCCCAAAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtca                                                                                                   <  1:382254/239‑1 (MQ=255)                                                                                                                                         tCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACTTTGGTCATCACGCCATa                                                                                         >  1:207454/1‑239 (MQ=255)                                                                                                                                                       ccgcccCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGCCATCACGCAATACAGAAAACAgcgca                                                                           >  2:318339/1‑235 (MQ=255)                                                                                                                                                       ccgcccCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCACCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGGCCTAACCCCAAAAAGCCGAAAGGGGa                                                                           >  2:425191/1‑239 (MQ=255)                                                                                                                                                                     ggTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGGGTAAAATCGACGGCTTTTCCCAACCCCAAAAAGTTGGGCATAACCCCATAAAGCAAACAGGCCACAACAAATCCAAcc                                                             >  2:868733/1‑238 (MQ=255)                                                                                                                                                                     ggTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGCCAATTAATCCTACg                                                             >  1:143098/1‑239 (MQ=255)                                                                                                                                                                     ggTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAAAAAAGCGAAACTTAAACCTAca                                                             >  2:222497/1‑238 (MQ=255)                                                                                                                                                                       tCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACgag                                                           >  1:201944/1‑239 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGCAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGGAAAATCGACGCCTTTTCCCAGCCCTAAAACGTTGGGCATCACCCCAAAAAGCAAAAAGGCGAAAAATAAACCCACggaaa                                                         >  2:504762/1‑235 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCCGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATAAATCCAAAGAGaa                                                         >  2:762367/1‑239 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACgaga                                                          >  1:330233/1‑238 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCAAGCCATACAGCAGACAGGCGACACTTAATCCTACGCGaa                                                         >  1:764674/1‑239 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCAAGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACCTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGaa                                                         >  1:576452/1‑239 (MQ=255)                                                                                                                                                                         cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGCAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGAAAAATCGCCGGCCTTTCCCAGCCCTAAAACGCTGGTCATCACGCCATAAAGCAGACAGGCGACGCTTAATCCTACGAGaa                                                         >  1:1051990/1‑239 (MQ=255)                                                                                                                                                                               tgcccccGGTCCCCCCGGGGACTTCCTTGAGTTGCCCTATGTTTGCCGCCCCCAAAA‑CCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:143098/234‑1 (MQ=255)                                                                                                                                                                               gtcccccGGTCACCCCGGCGCCTCCCTGGAGATGACTTGATTTTCCCCCCACAATGA‑TCCCCGACACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTCTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:415125/234‑1 (MQ=255)                                                                                                                                                                               gccGTCGGGTCACCCGGGCGATTGCCTTGAGAGGACCTAAGTTTTCCGCCACCAAAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  1:955186/237‑1 (MQ=255)                                                                                                                                                                                   tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGACAATTAATCCTACGAGAATCAATAAc                                                 >  1:273999/1‑238 (MQ=255)                                                                                                                                                                                   tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAAcc                                                >  1:47023/1‑239 (MQ=255)                                                                                                                                                                                       gTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAAGCAGCCTGGCTGACGCGAGGAAAATCGACGGCTTTTCCCGGCCCTAAAACGTTGGTCATCACGCCAAACAGCGAAAAGGCGACAATTCATCCTACGAGAAACAATAAcccccg                                            >  2:590637/1‑235 (MQ=255)                                                                                                                                                                                                    cTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGGGAAAATTAATCCTACGAGAATCAATAACCCCTGAATCATGCCCCgg                               >  1:188282/1‑239 (MQ=255)                                                                                                                                                                                                     tGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCACTAACCGCGGGATCATCCCACggg                              >  1:724085/1‑239 (MQ=255)                                                                                                                                                                                                      gCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGt                             >  1:656639/1‑239 (MQ=255)                                                                                                                                                                                                       ccTGGAGATCACTTAATTTTCCCGCCACCAGAC‑TCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTa                            <  1:717393/239‑1 (MQ=255)                                                                                                                                                                                                       catggagaTGACATAATTTTTCCGCCACTAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGa                                                                                                                                 <  2:268438/134‑1 (MQ=255)                                                                                                                                                                                                               cttccaacttttcCCCCCCCCAAAA‑CCCCCGCCCCGGGCGAAATCAGCACCATTGCCTGACCATTAAAAGCAGGAGGGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACAcc                      <  2:194016/228‑1 (MQ=255)                                                                                                                                                                                                                 aCCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGACCCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTcc                                                                                                                       >  1:268438/1‑138 (MQ=255)                                                                                                                                                                                                                 aCCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGACCATTAATCCTACAAGAATCAAAAACCGCTCGATCATGCCACAGGAAAAAAACGcc                   >  1:501964/1‑238 (MQ=255)                                                                                                                                                                                                                  cttAATTTTCCCGCCCCAAAAA‑CCACCGCCACGGGCCCAATCAGCATCATTTCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGccac                 <  2:743167/237‑1 (MQ=255)                                                                                                                                                                                                                  ccTAAGTTTCCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGccac                 <  1:60373/239‑1 (MQ=255)                                                                                                                                                                                                                      aGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGCCAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAg              >  1:743167/1‑238 (MQ=255)                                                                                                                                                                                                                         ttcccGCCACAAGAA‑CCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCa          <  1:222497/236‑1 (MQ=255)                                                                                                                                                                                                                         tatCCGCCCACAGAA‑TAACCGGCACGGGCGCAATCAGCAGCACTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGGGGCAGACCGAACCAGGCAGCATGGCTGACTAGAGTAAAATCGACGGCTTTTCCCATCCCGAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACCCCGCCACCAGCCCa          <  1:590637/237‑1 (MQ=255)                                                                                                                                                                                                                          ttCCCCCACCAAAAGTCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAAGGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCa          <  1:166159/239‑1 (MQ=255)                                                                                                                                                                                                                           tCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGgcagca                                                                                                                                                       <  2:1037609/96‑1 (MQ=255)                                                                                                                                                                                                                           tCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGgcagca                                                                                                                                                       >  1:1037609/1‑96 (MQ=255)                                                                                                                                                                                                                             cGCCACCAAAA‑CCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATac      <  1:425191/239‑1 (MQ=255)                                                                                                                                                                                                                                caccaGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTa                                                                                                               >  1:788697/1‑131 (MQ=255)                                                                                                                                                                                                                                caccaGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAA‑TAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAg  >  1:223544/1‑239 (MQ=255)                                                                                                                                                                                                                                cacaaGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTa                                                                                                               <  2:788697/131‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            CAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGTCCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCACCAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAG  >  NC_000913/1068838‑1069299

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