Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,358,901	G→T	P96Q (CCG→CAG)	puuP ←	putrescine:H(+) symporter PuuP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,358,901	0	G	T	95.0%	123.1 / ‑6.4	40	P96Q (CCG→CAG)	puuP	putrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (21/17);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGCAGACG  >  NC_000913/1358663‑1359139                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               gATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCACCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGGGAAAGATTTTCGCCCACAACACGTTAATCATCGGCCAAAAAAGACAATCCAACCGCGCAT‑ACCAGCCGACCCTAAAACCGACGGACTGGTTAATCTACt                                                                                                                                                                                                                                    >  1:15419/1‑251 (MQ=255)                gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGGCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGGCCATAAATCCGACGTGCTGGTTAATCGGCTTTTGCGCGGAGGta                                                                                                                                                                                                                      <  2:266213/250‑1 (MQ=255)                   ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACAGCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAg                                                                                                                                                                                                                   <  1:182042/250‑1 (MQ=255)                          tGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAAc                                                                                                                                                                                                            <  1:514262/250‑1 (MQ=255)                            cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATGGGGACGAAGGGTACCACCCCCACCCACCGCGGCACTTTCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGCCCATAAATCCGACGTGCTGGTTAATCGCCTTTTTGGCGTGGGAATAGGCCGGACCgg                                                                                                                                                                                                         >  2:468126/1‑251 (MQ=255)                                                    gTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGtt                                                                                                                                                                                 <  1:80143/251‑1 (MQ=255)                                                            caaGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACACCCACGGCGGCTCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCACCAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAAGCCGACGTGCTGGTTAAGCGCCTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAg                                                                                                                                                                         <  2:450365/249‑1 (MQ=255)                                                                 ttGACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGACGAAGGGTACCCCCCACAACCACGGCGGGACTTCCGGGAAGAGGGCGGAGAGAAAGATTTTCGCCAACAAGACGGTAATCAACTGGAAAAAAAGACAATCCAACAGCG‑ATGACCAGACGACCATAAATCCGACGTGCTGGTTAAACGACTTTTGCGCGTTGGTAAAAGCCGAACCGGCCTCCGGGAACTGGCGAACCCGATTGCCGTAGGTGa                                                                                                                                                                     >  2:312569/1‑250 (MQ=255)                                                                             aGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAg                                                                                                                                                        <  1:531599/251‑1 (MQ=255)                                                                                  cgcgGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCTTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCacac                                                                                                                                                   <  1:240707/251‑1 (MQ=255)                                                                                     ggcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACAcccc                                                                                                                                                >  1:92756/1‑250 (MQ=255)                                                                                              aaTGGCGACGAAGGTTAACACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgcc                                                                                                                                        >  1:170767/1‑250 (MQ=255)                                                                                                     acgaAGGTTACCACCCACACCAACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGCTTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTTGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGCCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                                                                 >  2:289458/1‑250 (MQ=255)                                                                                                      cgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGCTTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                                                                 >  2:163180/1‑249 (MQ=255)                                                                                                      cgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                                                                 >  1:360264/1‑249 (MQ=255)                                                                                                                   ccacacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCAAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACACGCATCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGCTAGCGGTAAAACGCAAACCCGCCAGCGCCAGCAAATAGGACGCCGGAACgg                                                                                                                  >  1:371335/1‑250 (MQ=255)                                                                                                                       ccgccccgcggccatcccGGGAAGAGGGCGGAGAGATAGATTTCCGCCAAAAATACGTGAATCCTCGGAAAAAAGGTACATTCCAGCTGCG‑AGGCCCAGCCGACCAAAAAGCCGACGCGCTGGTAACTCGAGTTGGGCGGGTAGGTATCGGCCGACCCGGCCTCGGGAAACTGGCGAACAATTTGGCCGGAGCTGAGAGCGGAAACCACCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGAcc                                                                                                              <  2:459664/235‑1 (MQ=255)                                                                                                                              gcggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTAAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa                                                                                                       <  1:43381/251‑1 (MQ=255)                                                                                                                                       ccGGGAAGAGGGGGGAGAGATAGATTTTCCCCAACAAGACGTAAATCATCGGAAAAAAGAGACAATCCAGCAGCT‑ATGCCCAGCCGACCATAAATCCGGCGGGCTGGTTAATCGACTTTTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGAAACTGTCGAACCAGTTTGCCGGGGCTGATAGCGGTAAACAGCACACCCGCCAGCTCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgag                                                                                              <  1:280632/251‑1 (MQ=255)                                                                                                                                           gAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                                                                                                                                                                                   <  2:89190/146‑1 (MQ=255)                                                                                                                                           gAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                                                                                                                                                                                   >  1:89190/1‑146 (MQ=255)                                                                                                                                                     gagagagaGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGCCCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAACCTGGCGAACCAGTTTGTCGTAGTTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGt                                                                                <  1:163180/251‑1 (MQ=255)                                                                                                                                                      ggagaAAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGTTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAAACGGGCGAACCATTTTGCCGGAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTa                                                                               <  1:417691/250‑1 (MQ=255)                                                                                                                                                          atatattttctcaaaaaaGACGTTAACCATCGGCAAAAAGATATAACCCAGCAGCG‑ATGACCAGCCGACAATAAATCCGAGGGGTTGGTTAATCGACTTTTGCGCGTAGGTTAAGGCCGAACCGGCCTCCGGAATCTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGGTAATGCTGGAGACAATACCAAAAGTATCaa                                                                           <  1:161728/235‑1 (MQ=255)                                                                                                                                                               ttttCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACCATACCAAAAGTATCAAATACg                                                                      >  2:549444/1‑251 (MQ=255)                                                                                                                                                                    gCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAACCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCAt                                                                 <  1:218748/251‑1 (MQ=255)                                                                                                                                                                     ccAACAAGACGTTAATCATCGGCGAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGCCGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCACAAGTATCAAATACGGCCATc                                                                >  2:74964/1‑251 (MQ=255)                                                                                                                                                                             aCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGCCGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                                                            >  2:381869/1‑215 (MQ=255)                                                                                                                                                                             aCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACAATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                                                            <  1:381869/215‑1 (MQ=255)                                                                                                                                                                                   atcatcGGCAAAAAGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGc                                                  >  1:92451/1‑251 (MQ=255)                                                                                                                                                                                         ggaaaaaaGAGACAATCCAGCAGCG‑AGGACCAGCCGCCAAAAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAcc                                            <  1:53917/248‑1 (MQ=255)                                                                                                                                                                                               aaGAGACAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCcatcat                                      >  1:191498/1‑251 (MQ=255)                                                                                                                                                                                                        cctatctggcaaTGACCAGCCGACCAAAAACCCGACGTGCTGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAACCGGGCGACCCAGTTTGCCGTACCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCt                              <  2:276906/240‑1 (MQ=255)                                                                                                                                                                                                             cagcG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGCAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGACAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGACCAGACCCATCATCACCACCTGcccc                         >  1:539298/1‑248 (MQ=255)                                                                                                                                                                                                              agcG‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCt                              >  1:276906/1‑244 (MQ=255)                                                                                                                                                                                                                 g‑ATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTCATGCCGGAGACAATACCAAAAGAATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAAttt                    >  2:197581/1‑251 (MQ=255)                                                                                                                                                                                                                  gaggcccttccgcccaaaaaTCCGACGTTCGTTTAAATTGACTTTTGCGCGAAGGTAAAGCCCGACCCGGCCTCCGGAACCTGGCGACCCAGTTTGCCTAAGCGGATCGCGGTAAACATCACCCCCGCCAGCGCCAGCAAATAGGACGCCTGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAAttt                    <  1:327499/234‑1 (MQ=255)                                                                                                                                                                                                                       ccAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCAGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTg               >  1:520304/1‑251 (MQ=255)                                                                                                                                                                                                                               ccATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGc       >  2:233696/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaTCCGACGTGCTGGTTAATCGACTTTTGCGAGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCACATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGCAGACg  >  1:27273/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               GATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCG‑ATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGCAGACG  >  NC_000913/1358663‑1359139

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