Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,279,034	Δ2 bp	coding (556‑557/1350 nt)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,279,033	0	T	.	97.5%	173.6 / ‑3.6	40	coding (555/1350 nt)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (21/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
*	NC_000913	4,279,034	0	A	.	97.4%	167.8 / ‑2.7	40	coding (556/1350 nt)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (20/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGT  >  NC_000913/4278786‑4279260                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                   aaCATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAttgtc                                                                                                                                                                                                                                  <  2:303182/251‑4 (MQ=255)   cATATTGGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGTTCTGGTGATTAAAAACCCGCTTGATGGTCTTCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAttgtcac                                                                                                                                                                                                                                <  2:58600/251‑6 (MQ=255)     tatTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGAGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTCGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCCTTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAttgtcactg                                                                                                                                                                                                                              >  1:390002/1‑244 (MQ=255)            gTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTc                                                                                                                                                                                                                     >  1:375945/1‑251 (MQ=255)             tACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTc                                                                                                                                                                                                                     >  2:82111/1‑250 (MQ=255)              aCCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCat                                                                                                                                                                                                                   <  1:166092‑M1/251‑3 (MQ=255)                  gTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacc                                                                                                                                                                                                                <  1:212986‑M1/250‑6 (MQ=255)                  gTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatac                                                                                                                                                                                                                 >  1:72787‑M1/1‑245 (MQ=255)                  gTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTGGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatac                                                                                                                                                                                                                 <  2:72787‑M1/249‑5 (MQ=255)                               ccGTGGTCCAGATGGGTGTGTTGTTTGCGTAAATTTCCGCTCA‑GGGTATCCGTAGCTGTATCTTGCGCAATTTGCCTCCCTGTGTGGCGCATGGCACGGGGATTGGGATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCTGTTGAGGTTCTGCCCGGTGCGCTGAGTGATTTCTGGACCTCCACGGTGAT‑‑TGTCACTGGTAGGTcatacctaacgacaacagg                                                                                                                                                                                                   <  2:175244/250‑19 (MQ=255)                                        tGATGCGTGTGCTGTTTACGCTAATTTCTGC‑CAAGGGTATCCGTAGTTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCAAGGCACGGGGATTGGTCTAGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGCTGGTCTGCCCGTTGCGCTGGGAGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaccc                                                                                                                                                                                         >  2:201382‑M1/1‑223 (MQ=255)                                          aTGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccat                                                                                                                                                                                       >  2:316497‑M1/1‑221 (MQ=255)                                                         aCGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtg                                                                                                                                                                        >  1:71638‑M1/1‑206 (MQ=255)                                                           ggTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccataccccaaaggtagtga                                                                                                                                                                       >  2:122880‑M1/1‑204 (MQ=255)                                                                  tCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCACCTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTATGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCataactaacgacaacaggaaaaaaaaccatacccaaacggcagtgaatcgtttg                                                                                                                                                               >  2:397945‑M1/1‑197 (MQ=255)                                                                   tctgccgagggtatccgtaGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGGTTGGTATCGGTCTGTTCATGCTGCTCATTGCCGCGTACGGTGTCGGGCTGGTGATTAAAAACCCGCTTGAGGGTCTGCCCGTTGCGCTGGGTGATTTCGCAGCCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggataaaaaaccagacccaaaagggagtgaatcgttt                                                                                                                                                                <  2:35245‑M1/242‑54 (MQ=255)                                                                          cGGGTATCCGTAGCTGGATTGTGCGCAACTTGCCTCACGGGGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCTCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGGCTGGCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTG‑TT‑TGTCACTGGTAGGTCatacctaacttcagcgggataaaaaaccatacccaaaaggtagtgaatcgtttgcttttaa                                                                                                                                                        <  1:180046‑M1/251‑62 (MQ=255)                                                                           gggTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGGCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaac                                                                                                                                                       <  1:166202‑M1/251‑63 (MQ=255)                                                                               aTCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggg                                                                                                                                                   >  1:410237‑M1/1‑185 (MQ=255)                                                                                     aGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgt                                                                                                                                             >  1:255715‑M1/1‑179 (MQ=255)                                                                                       cTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTACCGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttc                                                                                                                                           >  2:420358‑M1/1‑177 (MQ=255)                                                                                         ggATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcag                                                                                                                                         <  1:58050‑M1/251‑77 (MQ=255)                                                                                                    accttGCCTCTCGGTGTGGCGCACGGCAGCTGGACTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAAAGGTGTCGGTCTGGTGATTCAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacagggaaaaaaaccataaccaaacggtagtgaatctattgcttttaaatggggtatgttcaggtggcgttgcg                                                                                                                              >  2:245151‑M1/3‑164 (MQ=255)                                                                                                     aCTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGAATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggcacaaaaaccatacccacaaggaagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcgg                                                                                                                             >  2:20109‑M1/1‑163 (MQ=255)                                                                                                      cTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggg                                                                                                                            >  2:121093‑M1/1‑162 (MQ=255)                                                                                                        tGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcgg                                                                                                                             <  1:121093‑M1/248‑89 (MQ=255)                                                                                                         gCCTCACTGTGTGGCGCACGTCACGGGGATGGGTATCGGTCTGTTCCTGCGGCTCATTGCCGCGAACGGTGTCGGTCTGGTGAGTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCTGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtgtgcgtttaacttgggattgttcatgtggcgttgcggtgtt                                                                                                                         <  1:165305‑M1/251‑93 (MQ=255)                                                                                                                 gtgtGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgt                                                                                                                 <  1:25299‑M1/251‑101 (MQ=255)                                                                                                                  gggggcgcgcgtaacgggtgttggtatgggtctggtcctgctTCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCACTTGATGGTCTGCCCGTTTCGCTGGGTGAGTTCGCGACCTTCCCGGTGAT‑‑TGGCACTGGTAGGTCaaacctaacgacaacagtaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgtt                                                                                                                <  2:893‑M1/231‑102 (MQ=255)                                                                                                                   gtgGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgtt                                                                                                                >  2:144797‑M1/1‑149 (MQ=255)                                                                                                                            ggcacggGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgg                                                                                                      >  1:134648‑M1/1‑140 (MQ=255)                                                                                                                                  gggATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggaggga                                                                                                <  2:321969‑M1/251‑118 (MQ=255)                                                                                                                                           aTCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggggttttgaatgttaagactgtgggagggagaatctgtg                                                                                       <  1:397945‑M1/251‑127 (MQ=255)                                                                                                                                            tCGGTCTGTTCCTGCTGCTCAGTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtg                                                                                       <  2:154945‑M1/250‑127 (MQ=255)                                                                                                                                                tCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaacca                                                                                                                                                                                        <  2:28633‑M1/149‑30 (MQ=255)                                                                                                                                                tCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaacca                                                                                                                                                                                        >  1:28633‑M1/1‑120 (MQ=255)                                                                                                                                                     tCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgattgttaagactgtgggagggagaatctgtggccggatccg                                                                             >  2:28751‑M1/1‑115 (MQ=255)                                                                                                                                                           tgctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaactacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgttatcggttaactggggattgttc                                                                                                                                           <  2:192417‑M1/183‑75 (MQ=255)                                                                                                                                                           tgctCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcgtttgcttttaactggggattgttc                                                                                                                                           >  1:192417‑M1/1‑109 (MQ=255)                                                                                                                                                              tCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCAGCTTGATGGTCTGCCCGTTGCGCTGGATGATTTCGCGACCTCCCCGGTGAT‑‑TGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggaagtgaatcgtttgcttttaactggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggagggagaatctgtggcaggaaccgcctctgggt                                                                    >  1:88757‑M1/1‑106 (MQ=255)                                                                                                                                                                  tGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGA‑‑GTGTCACTGGGAGGTCatacctaacgacaacaggaaaaaaaaccatacccaaaaggtagtgaatcggttgcttttcactggggataggtcaggtagcgttgcggtgttttgaatgttaagactgggggagggaggatctgtggcaggaaacgccgctggtatggg                                                                >  1:419315‑M1/1‑102 (MQ=255)                                                                                                                                                                                              aaaaaCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTCCCCGGTGA‑‑ATGTCACTGGTAGGTCatacctaacgacaacaggaaaaaaaaccataccaaaaaggtagtgaatcgtttgctttaaactggggattgttcaggtggcgttgcggtgttttgaatgttaagtctgtgggagggagaatctgtggcaggaaccgcctctggtatagggggaggctaagatagcattatttcctct                                    <  1:316497‑M1/251‑178 (MQ=255)                                                                                                                                                                                                                                gTGATTCCGCCACCTCCCCGGTGAT‑‑TGTCACGGGTAGGTCatcccaaacgacaacagaaaaaaaacccatccccaaaaggtagtgaatcgtttgcttttaacgggggattgttcaggtggcgttgcggtgttttgaatgttaagactgtgggggggaaaatctgtggcaggaaccgcctctggtatagggggaggcgaagatagcattatttcctctccctcttcttcgctgtaccagataaccggcacca  <  1:300976‑M1/251‑212 (MQ=255)                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                   AACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGC‑CACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGGTCTGTTCCTGCTGCTCATTGCCGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGT  >  NC_000913/4278786‑4279260

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