Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,538	A→G	F56F (TTT→TTC)	cyoE ←	heme O synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,538	0	A	G	92.1%	113.7 / 0.4	38	F56F (TTT→TTC)	cyoE	heme O synthase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (19/16);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

AATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGAT  >  NC_000913/447292‑447779                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aaTCAACGTGCCGTAGACAGAGTGGCGTTGCATGTACAGGCTAGAAACGCCGACATAAACCACAAAGCCCATCGCCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCATCACCCGAGTCTTCGGCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAaca                                                                                                                                                                                                                                                <  2:323727/251‑1 (MQ=255)          gCCGTAGACAGAGTGGCGTTTCATTTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCGGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCg                                                                                                                                                                                                                                       <  1:196665/251‑1 (MQ=255)            cGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGAc                                                                                                                                                                                                                                     <  1:92268/251‑1 (MQ=255)               aGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGcc                                                                                                                                                                                                                                  >  1:305777/1‑251 (MQ=255)                 acagatgggCGTTGCATGTACAGGATAAAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCATCAGGCCAGTTGATGCGCGCCAAACCACAGCAGCATAAAGCCAGTAATGCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACTAGCTGCCGATTCGTCGCCCGTGCCAGCTGTTTGGCGATATCCCTGTCGAGTTAGTTGTTGAAcacacacccagatgcctc                                                                                                                                                                                                                                <  2:236899/244‑13 (MQ=255)                                                                    cccATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAAtgct                                                                                                                                                                            >  1:302376/1‑251 (MQ=255)                                                                      cATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAAtgctgc                                                                                                                                                                          >  1:131886/1‑251 (MQ=255)                                                                           ccaccatccagcAGGCCAGCGGATTCGCGCCAAACCCCAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGTCAGCATGAGAGATCAGGCCTTGCACCAGCACCCGATTCGTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACTCACAACCCGACGCGACAACCAGTGACTCCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGccc                                                                                                                                                                        <  2:207670/241‑1 (MQ=255)                                                                                    agcagGCCAGCGGATTCGCGCCAAACCACAGCAGCAGAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATGCTTCGTCCTTTCCATCTGTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACA‑CCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCcagca                                                                                                                                                           <  2:346445/251‑1 (MQ=255)                                                                                          ccAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAAt                                                                                                                                                      >  1:149459/1‑251 (MQ=255)                                                                                          ccAGCGGATTCGCGCCAAACCACAGCAGCAGAAAGCCAGCATTACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCCTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAAt                                                                                                                                                      <  2:343573/251‑1 (MQ=255)                                                                                             gCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATccc                                                                                                                                                   <  1:78941/251‑1 (MQ=255)                                                                                                      cgcCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGAGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCCATCAcc                                                                                                                                          >  1:40417/1‑251 (MQ=255)                                                                                                      cgcCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAAAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCAcc                                                                                                                                          >  1:263601/1‑251 (MQ=255)                                                                                                       gcCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGCTAATCAATGCTGCCCTTTGAGGCCAGCAGGCATCCCCCAATCACCg                                                                                                                                         >  2:28197/1‑251 (MQ=255)                                                                                                               acaGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAgg                                                                                                                                 >  1:346359/1‑251 (MQ=255)                                                                                                                      caTAAAGCCAGCAATACCCAGCAACGTGGCGTACAAAAGAGAGACCGCCGGAGAGCTCAGGCCTTTCAACAGCACCCGATTCTTCGTCCTTTCAATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCCATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCa                                                                                                                           >  1:178049/1‑250 (MQ=255)                                                                                                                            gCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAgatg                                                                                                                     >  2:111789/1‑250 (MQ=255)                                                                                                                                 cAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCt                                                                                                               >  1:257278/1‑251 (MQ=255)                                                                                                                                 cAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGACAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCt                                                                                                               >  1:25305/1‑251 (MQ=255)                                                                                                                                      cccAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                                                                                          >  1:168044/1‑251 (MQ=255)                                                                                                                                      cccAGCAACGGGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTGCCATCTTTCTGGCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGttt                                                                                                          <  2:131886/251‑1 (MQ=255)                                                                                                                                         aGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCGCAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATTCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGt                                                                                                       <  1:152112/251‑1 (MQ=255)                                                                                                                                           cAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTa                                                                                                     <  1:264993/251‑1 (MQ=255)                                                                                                                                             aCGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGGCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGATGCCAGCAGGAATCCCCCAAGCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACt                                                                                                   <  2:20061/251‑1 (MQ=255)                                                                                                                                                  gCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAg                                                                                              <  1:111789/251‑1 (MQ=255)                                                                                                                                                      actccAGCGAGACAGCAGGGGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTTCAGGTAt                                                                                          <  2:31528/248‑1 (MQ=255)                                                                                                                                                        acCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTg                                                                                        <  2:62384/251‑1 (MQ=255)                                                                                                                                                            gCGAGACAGCAGGAGAGATCAGGCCTGTCCCCAGCACCCGATTCTTCGTCCTTTCCAGCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAAGCGCCGAGATCAGGGTGCCAAAGATGATGCCTGGTTTCTTGACTTGCAGTTATTGCTTa                                                                                    <  2:169868/251‑1 (MQ=255)                                                                                                                                                             cGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGCGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCCGGTATTGCTTa                                                                                    >  2:327204/1‑250 (MQ=255)                                                                                                                                                                                aGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATc                                                                                                                                    <  2:125161/183‑1 (MQ=255)                                                                                                                                                                                aGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATc                                                                                                                                    >  1:125161/1‑183 (MQ=255)                                                                                                                                                                                                         tcGTCCTTTCCATCTTTCTGCCGATACCCCTGGCGATGTAGTTGTTAAACACACACCCCGCCGCCACAACCAG‑GACCCCCCCACCCAGCGTGAAGATAAACGGGGGAAAATCAATGCTGCCCTTGGAGGCCAGAAGAAACCCCCAAATCACCGAGATCAGGTTCCAAAAGTTGATGCCTGTTTTCGTTCCTTGCAGGTATTGCTAAAACATCAAAACCGCCGCTCTTAGTGCATCATCATGTTGTAGt                                          <  1:268137/248‑1 (MQ=255)                                                                                                                                                                                                            tCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCA‑CACCCGCGCGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGACAGCAGGCATCACCCAATAACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATAGCTTAAACATAATAACCGCCGCTCTTAGTGCATCATCATGTTGTCGTTGCGGt                                    >  1:221375/1‑251 (MQ=255)                                                                                                                                                                                                                      cTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCACAACAAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGACTCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGGTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAAt                          >  1:161947/1‑251 (MQ=255)                                                                                                                                                                                                                                     ccTGTCGATGTAGTTGTTGAACACACAACCCGACGCAACAACCAGTGACACCCC‑AACCAGCGTGAAGATAAACAGGGGATAATCAATGATGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATACTCCAGATGGAGCCTAc           >  2:156769/1‑251 (MQ=255)                                                                                                                                                                                                                                              gTAGTTGTTGAATACACAACCCGACGCCACAACCAGTGACACCCC‑AACAAGCGTGTAGATAAACATGGGCTAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCACCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTACAGGTATTGCTTAAAAATCATAAACGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGACTACAACAAGGAt  >  2:309191/1‑251 (MQ=255)                                                                                                                                                                                                                                               tAGTTGTTGAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGAt  >  2:269044/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   AATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCACAACCAGTGACACCCC‑AACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGAT  >  NC_000913/447292‑447779

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