Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,276	G→A	A16A (GCG→GCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,276	0	G	A	84.6%	84.3 / 4.8	41	A16A (GCG→GCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (14/19);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGT  >  NC_000913/4618034‑4618520                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      tCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCg                                                                                                                                                                                                                                                <  1:267468/250‑1 (MQ=255)  ccAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGATGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGAt                                                                                                                                                                                                                                              >  1:113592/1‑251 (MQ=255)     gCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATga                                                                                                                                                                                                                                            >  1:110091/1‑250 (MQ=255)         gctgGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATgaagaa                                                                                                                                                                                                                                        <  1:263516/250‑1 (MQ=255)           gggctaGCCTGTTGAAAGCTCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAAtt                                                                                                                                                                                                                                     <  1:207655/246‑1 (MQ=255)             gCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCg                                                                                                                                                                                                                                   <  1:225642/251‑1 (MQ=255)                  catgttgAAAGCGCTGGTCCACAGGGACGATAAGAGCGCCAGCAGCTACTAAGTAAGAAGCTTTACGCCTGGTTCGTTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGCTCGGATAAGGGGTAGACGCATCCGGCAAAAGCCTCAGCATACGCTTTTCCTCGGGAGTTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTTATGTACTTAGCCATGTAGAAATTCGTTTCt                                                                                                                                                                                                                              <  2:53209/246‑1 (MQ=255)                                             cGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTa                                                                                                                                                                                                   <  1:70211/251‑1 (MQ=255)                                               gTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAGGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtc                                                                                                                                                                                                 >  1:116576/1‑251 (MQ=255)                                                tAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGCATTCGCGCCAAAACTAtct                                                                                                                                                                                                >  1:50402/1‑251 (MQ=255)                                                                       taagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTCCCTTGTTTCTCGCACAAGAAATTATTC‑GTACAAAACGTGATGGTCATGCGATGAGCGATGAAGAAATTCGTTTCTTTATCAACAGAATTCGCGCCAACACTCACTCCGATGGGCAGATTGCCGCCCTc                                                                                                                                                                         >  2:197277/1‑251 (MQ=255)                                                                        aagATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcg                                                                                                                                                                        <  1:78225/251‑1 (MQ=255)                                                                                    gCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAttt                                                                                                                                                            >  1:17340/1‑251 (MQ=255)                                                                                       gtgtgctctgcgctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGGAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCTTCGCGATGACCAGTttct                                                                                                                                                         <  2:231664/244‑1 (MQ=255)                                                                                                  gctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTAtt                                                                                                                                                                                                                >  2:6593/1‑185 (MQ=255)                                                                                                  gctTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTAtt                                                                                                                                                                                                                <  1:6593/185‑1 (MQ=255)                                                                                                       ccaatcataagagtcGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACCCTATCTCCGAAGGGCAGATTGCCGCACTCGCGATGCCCATTTTCTTCCACGATATGACAAt                                                                                                                                         <  2:161852/237‑1 (MQ=255)                                                                                                                                       ggATAAGGCGTAGACGCATTCGGCAAAAGCCGCCTCATACTCTTGTCCTCGGGAGGTTACCGTGATGCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCCCTGAGCGATGAAGAAATGCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGTCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGc                                                                                                         <  2:108886/251‑1 (MQ=255)                                                                                                                                                           cGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGATCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTtctc                                                                                    <  1:235992/251‑1 (MQ=255)                                                                                                                                                                     cgtcTCATACTCTTGTCCTCGGGAGGTTGCCTTGTTTCTCTAGCAAGAAATTTTTCTGTAAAAT‑CGTGATGGTCATTCACTGAGCGATGAAGAAATTCGATTCTTTATCTACGGCATTCGCGACAATACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCGGAGCGTGTCTCGCTGACCATGGCGATTCGAGA‑TTCAGGATCCGTTCTCTCCTGGaaaa                                                                          <  2:257366/248‑1 (MQ=255)                                                                                                                                                                          cATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCt                                                                      <  1:282715/250‑1 (MQ=255)                                                                                                                                                                          cATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTg                                                                     <  1:104507/251‑1 (MQ=255)                                                                                                                                                                           ctACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAATACCCGGGATTACAATCACCTGAGCGATGAAGAAATTCATGTCTTTAGAAACGGCATTCGCGACAACACTGTCTCCTAAGGGCAGGTTGACGCCATAGCGAGTACCATTTTCTTCCACGATAAGACCATGCCCGCGCGAGTCTCGCTAACATTGGCGATAC‑AGACTTCAGGAACCGGTCTCGACTGGAAAAGccgg                                                                     >  2:88629/2‑248 (MQ=255)                                                                                                                                                                            tACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGc                                                                    >  1:189315/1‑250 (MQ=255)                                                                                                                                                                                     ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATgg                                                          >  1:33187/1‑251 (MQ=255)                                                                                                                                                                                     ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATgg                                                          >  1:181613/1‑251 (MQ=255)                                                                                                                                                                                        cGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTCTCTCCGACGGGAAGCTTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAAGGCCTGCACTTGTCTCGCTGACCTTGGCGATGCGAGC‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCAGCTGCATGGccc                                                       >  2:99295/1‑251 (MQ=255)                                                                                                                                                                                                tACCGTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgtt                                                <  1:77393/250‑1 (MQ=255)                                                                                                                                                                                                 aCCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCAGGCACTGAGCGATGAAGAAATTCGTTTCTTTTTCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgttg                                               <  2:134946/250‑1 (MQ=255)                                                                                                                                                                                                            cGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTcc                                    <  1:374673/250‑1 (MQ=255)                                                                                                                                                                                                               acaAGAAATTATTC‑GTAAAAAACGTGATGGGCATGCGCTGAGCGATGAGGAGATTCGTTTCTTTATCAACGGTATCCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACAATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACAATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAACCCCTGAATCTGAATGGCCCGATTGTTGATAAACACTCCACCg                                >  2:152873/1‑251 (MQ=255)                                                                                                                                                                                                                        tattC‑GTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTc                        <  1:92958/250‑1 (MQ=255)                                                                                                                                                                                                                        tattC‑GTAAAAAACGTGAGGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAGCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTc                        <  2:181613/250‑1 (MQ=255)                                                                                                                                                                                                                               tAAAAAACGTGATGGTCATGCACTGAGAGATGAAGCAATTCGTTTCTTTATCAAAGGCATTCGCGACAACACTATCTCCAAAGGGCAAATTGACGCCCTCGCGATGACCATTTTCTTCCACGATATGAACATGCCTGAGCGTGTCTCGCTGACCATGACGATGCGAGA‑TTCAGGACCCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGGATGTTGATAAACATTCCACTGGTGTCGTAGGCgctg                 >  2:384639/1‑248 (MQ=255)                                                                                                                                                                                                                                  aaaaCGTGATGGTCTTTCGCTGAGCTATTAAGAAATGCTTTTCTTAATAAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAGTTTCTTCCACGATATGACAATGCCTGAGCGTGCCGCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGa              <  2:110091/251‑1 (MQ=255)                                                                                                                                                                                                                                  aaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAA‑GCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGAc             >  1:63836/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaaCGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGAc             <  1:135826/251‑1 (MQ=255)                                                                                                                                                                                                                                   aaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGc                           >  2:178028/1‑237 (MQ=255)                                                                                                                                                                                                                                   aaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGc                           <  1:178028/237‑1 (MQ=255)                                                                                                                                                                                                                                     aCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACtt           >  1:197205/1‑251 (MQ=255)                                                                                                                                                                                                                                              tCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGt  >  1:348669/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTC‑GTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGA‑TTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGT  >  NC_000913/4618034‑4618520

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: ‑