Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,681,622	T→A	13.9%	L216Q (CTG→CAG)	folM →	dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,681,622	0	T	A	13.9%	109.0 / 9.4	43	L216Q (CTG→CAG)	folM	dihydromonapterin reductase, NADPH‑dependent; dihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base A (2/4);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGT  >  NC_000913/1681328‑1681922                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                              aTGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGa                                                                                                                                                                                                                                                                                                         <  1:472796/301‑1 (MQ=255)      cTGGTAAGATTACTGCTTGGTCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcact                                                                                                                                                                                                                                                                                                    <  1:167679/301‑1 (MQ=255)          tttcagacgtgtgctcttccgatctcgccgcCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttctt                                                                                                                                                                                                                                                                                                <  2:26240/276‑4 (MQ=255)          aaaTATTACTGCGTGGGCACTGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttctt                                                                                                                                                                                                                                                                                                <  2:530539/301‑4 (MQ=255)           aaGATTACTGCGTGGGCACGGTCACGCCGCCTGCGATATCAGTCACTTTACCGATTATGTGGGGGAGTGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTCGAtcacttcttt                                                                                                                                                                                                                                                                                               <  1:407001/301‑5 (MQ=255)                           gaCGGACACGCCGCCAGGGATATCATTCACTTTACGGAGTATGTGGGGGAGCGCGGTAGCGAAAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCa                                                                                                                                                                                                                                                                               <  1:396339/300‑1 (MQ=255)                                   cgccgcCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTAATCGACCAGGTAGAtcacttcttaaatcggaaaagcgtcggggaaggg                                                                                                                                                                                                                                                                       >  1:26189/1‑273 (MQ=255)                                         aggCGATATCATTCACTTTACCGATGATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTAGGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTc                                                                                                                                                                                                                                                                 <  1:105410/299‑1 (MQ=255)                                             tttattatTCACTTTACCGATTGTGTGGTGGGGCGCGGTAGCGACAAACATTTTGCGTATGCTGCAAGCAAAGCGGCACTGGTTAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCAc                                                                                                                                                                                                                                                             <  1:339968/299‑1 (MQ=255)                                               tatCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGATTTGTCACAGGACGCAGTTTCCCaaaa                                                                                                                                                                                                                                                           >  2:340431/1‑298 (MQ=255)                                                     tctcTTTACTGATTGTGTGGTGGGGGGCGGTAGCGACAAACATTTTGCTTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGc                                                                                                                                                                                                                                                     <  1:112846/298‑1 (MQ=255)                                                            aCCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCGAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGgtcgtc                                                                                                                                                                                                                                              >  2:299075/1‑301 (MQ=255)                                                                     gtggtgGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTCATGAACTTGATGATGCCGCATCTCTACCACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGCCCTGGTCGATTACTTACTTACCAGTTGCTTTGTCCCCGGACGCAGTTTCCCCCTTTATGGCGATCGTCCTCTCCGtt                                                                                                                                                                                                                                     >  1:556552/1‑301 (MQ=255)                                                                     gggggtgggcggGGTAGCGACAAACATATTGTGTATGCTGAAAGCAAAGCGGCACTGGATTATATGACCCGCTCGTTTGCCCGCATGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGtt                                                                                                                                                                                                                                     <  1:165435/293‑1 (MQ=255)                                                                      tggtggAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTa                                                                                                                                                                                                                                    >  2:165587/1‑301 (MQ=255)                                                                                    aGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGTCTGTCGTCATCTGCGTTAATACAGCTTTAtac                                                                                                                                                                                                                      >  1:106498/1‑299 (MQ=255)                                                                                        aCAAACATATTGCGTATGTTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑Ca                                                                                                                                                                                                                 <  2:541294/301‑1 (MQ=255)                                                                                           aaCATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATATGCGTTAATGCAACTTTAACCAA‑Caaat                                                                                                                                                                                                              >  1:514219/1‑301 (MQ=255)                                                                                                                         cACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAGGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAaca                                                                                                                                                                                <  2:148028/301‑1 (MQ=255)                                                                                                                                taataTGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACa                                                                                                                                                                         <  2:123707/301‑1 (MQ=255)                                                                                                                                           gCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGTATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCAAG‑CATATAATCAACACCAATGCGCTTATCCCCCAACACGAAACAAAAACaaaaaa                                                                                                                                                              >  1:49440/1‑301 (MQ=255)                                                                                                                                                   gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATACGCTTATCCCCAAACACGCAACAGAACCAACAAAAAcggcgg                                                                                                                                                      >  1:477272/1‑301 (MQ=255)                                                                                                                                                   gCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAACGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAG‑ACCAACAAAAGCGGCggg                                                                                                                                                     >  2:636534/1‑301 (MQ=255)                                                                                                                                                            cTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGANAGAAGTGATCGANNNNNNNNATTACNNNNNNNNCAGTTGCTTNNNCACCGGACGCANNNNNNNACTTGATGGAGNNNNNNNTCTACNNNNNNACA‑CATTATNCNNACATATAATACAAACCCAACCGCATATCCCACAACACACCACAAAACCAACaaaaacccacgcacaacgac                                                                                                                                             >  1:131616/1‑285 (MQ=255)                                                                                                                                                                cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGcatcat                                                                                                                                         >  2:236562/1‑301 (MQ=255)                                                                                                                                                                         tgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGCAAAGCGCATCATTCCGGTGaa                                                                                                                                >  2:517912/1‑301 (MQ=255)                                                                                                                                                                                  tctttccctacactacgctcttccgttctAATGAACATGATGATGCCGAAATTCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAg                                                                                                                       <  1:209159/273‑1 (MQ=255)                                                                                                                                                                                   ttCTATTGAGCCATGGCGTATCCTGGTTAATGAACATGATGATGCAGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGGGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTGCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGc                                                                                                                      <  1:236296/301‑1 (MQ=255)                                                                                                                                                                                      tATTGCGCCTTCGCTGTTCCTTTTTAATGAACATGTTGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGGGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGac                                                                                                                   <  1:298671/301‑1 (MQ=255)                                                                                                                                                                                               aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTaa                                                                                                          >  2:468237/1‑301 (MQ=255)                                                                                                                                                                                                 cGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATAACACAGCGAAACCAGGTaaaa                                                                                                        >  2:171567/1‑300 (MQ=255)                                                                                                                                                                                                               aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCAAAgagatcggaaaagcacacgtctgaactcca                                                                                          >  2:209402/1‑275 (MQ=255)                                                                                                                                                                                                               aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGCTGAAAACCCACAGCGACACCAAGAAAACAAAATAGGAaataa                                                                                          >  1:52404/1‑301 (MQ=255)                                                                                                                                                                                                                  gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACGTATGAAAACCGCGCCGGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCg                                                                                       <  1:517240/301‑1 (MQ=255)                                                                                                                                                                                                                   aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTAGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTCCTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGATAATGCAGCCTTATCCAG‑CATATAATCACCACCCATGCGCTTATCCCTCAACACGCAACAGACCCAACAAAAGCGGCGGGCACCACCACCATTCCGTTGCAATACCACAGCAACACCACGTAAACAAAAACGGGacaaaccgc                                                                                      >  1:488363/1‑293 (MQ=255)                                                                                                                                                                                                                                 gAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGCAAACAAAAAAGGGAATAATCAACAAAAAACTaaaaa                                                                        >  1:123327/1‑301 (MQ=255)                                                                                                                                                                                                                                  aatgtggACAACAGGTGCTGAATAAATCTCTGATGAAATCCGCGTCTGGCGGGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAat                                                                       <  1:467673/295‑1 (MQ=255)                                                                                                                                                                                                                                            acaGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTgcg                                                             <  2:514891/301‑1 (MQ=255)                                                                                                                                                                                                                                                        tAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAACGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATAGTTGCGCGTAAGACTTCa                                                 >  1:426402/1‑301 (MQ=255)                                                                                                                                                                                                                                                                            cgcCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCCCTGGcaac                             >  1:414612/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                      aaaagaaGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgc                   <  2:49508/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                         gtaagTAATCGACCAGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgcgca                >  2:392724/3‑300 (MQ=255)                                                                                                                                                                                                                                                                                               gATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGCTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCAATAAgc          >  2:592961/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                   gtccTGGTCTATTACTTACTTACCATTTGCTTTGTCACCGGACGCAGTTTCCCACTGGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGCTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCaa      <  1:592200/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       tGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgcgcttatcgcgcattata  <  2:477873/301‑16 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                              ATGCGCTGGAAAGATTACTGCGTGGGCACGGACACGCCGCCAGCGATATCATTCACTTTACCGATTATGTGGTGGAGCGCGGTAGCGACAAACATATTGCGTATGCTGCAAGCAAAGCGGCACTGGATAATATGACCCGCTCGTTTGCCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGCTTTGTCACCGGACGCAGTTTCCCACTTGATGGCGGTCGTCATCTGCGTTAATGCAGCTTTATCCAG‑CATATAATCAACACCCATGCGCTTATCCCCCAACACGCAACAGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGT  >  NC_000913/1681328‑1681922

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