Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	499,060	T→G	D305A (GAC→GCC)	aes ←	acetylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,060	0	T	G	98.2%	171.1 / ‑5.7	55	D305A (GAC→GCC)	aes	acetylesterase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (29/25);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCG  >  NC_000913/498816‑499302                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCTGTAAAGATCTGAGCGCCGGCGCGAATAGCCgcggcggcg                                                                                                                                                                                                                                               <  2:29696/250‑1 (MQ=255)     gaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAAAACTGAGCGGCGTCGCGGAAAGCCGcggcggcggggtt                                                                                                                                                                                                                                          >  1:321241/1‑247 (MQ=255)                 ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCGTCGTCGGCGGTTTTCATCATCCGTGt                                                                                                                                                                                                                              <  1:120028‑M1/251‑2 (MQ=255)                     atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATAACGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCCCGGCGGCGGGTTTCATCATCCGTGtagtc                                                                                                                                                                                                                          >  1:427266‑M1/1‑246 (MQ=255)                           aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGGGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAAAGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGGCGCGTATCGCTAAAGCTGAGGGGTAAAGAACTTGGCGCCGGCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGCGtagtcgtcggg                                                                                                                                                                                                                    >  1:98157‑M1/1‑240 (MQ=255)                                        gggCTTTGCTGCGGATTGTATGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAAAAAATTCCGGCGCTTAATGACACGCCGGAACATATCGAAATGGTGGCCCATCTTGTTGCCGCGGATCGCCTAAGCTTAGCGGGAAAAAACTGGGGGCCCTCGCGGAAAGCCTCGGCGGCGGGTTTTATCATCCGTggcggcgtggggccttgacattttt                                                                                                                                                                                                       >  2:428575/1‑227 (MQ=255)                                                       ttGTCTGGAGCCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAAAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATACGTGtaggcgtcggcatttgccatttttgcctgatgcgacgct                                                                                                                                                                                        >  1:109602‑M1/1‑212 (MQ=255)                                                                                   gAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAGAAAATATCAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGAGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCCTCATACGTGtagtcgtaggcattagacatttgtgcctgatgcggcgcttgacgcgtcttttcagggctacaaacgg                                                                                                                                                            >  1:180192‑M1/1‑184 (MQ=255)                                                                                            cGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATAGTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatc                                                                                                                                                   <  2:33235‑M1/251‑77 (MQ=255)                                                                                                      tCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattg                                                                                                                                          <  1:188744‑M1/250‑86 (MQ=255)                                                                                                      tCCTCGGTGCCGGCGGGAAAAAATATGAATATATGCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCTGTAAAGAACTGAGCGCCGTCGCGAAGAGCCGCGGCGGCGTTTTTCATCATCCGTGtagtcgtaggcatgaggcatttgtgcctgatgcgacgcttgatgcgtcttatcaggcctacacccggtgccggatccggcaat                                                                                                                                            <  2:200824‑M1/248‑84 (MQ=255)                                                                                                                cggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATTTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGCAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagctgtaggcattagacatttgtgcctgatgcgccgcttgacgcgtcttattaggcctccaaccggtgccgcatccggcaattggtgcacaacgc                                                                                                                               >  1:237532‑M1/1‑155 (MQ=255)                                                                                                                cggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgc                                                                                                                               <  1:257560‑M1/251‑97 (MQ=255)                                                                                                                cggcggGAAAAAATATGAATATATGCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgc                                                                                                                               <  1:293192‑M1/251‑97 (MQ=255)                                                                                                                 ggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGGGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggggcacaatgcc                                                                                                                              >  2:329926‑M1/1‑154 (MQ=255)                                                                                                                      gAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatg                                                                                                                         <  1:353104‑M1/251‑103 (MQ=255)                                                                                                                       aaaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgc                                                                                                                        <  1:9286‑M1/251‑104 (MQ=255)                                                                                                                         aaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACCTATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAACGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcg                                                                                                                       >  1:404213‑M1/1‑146 (MQ=255)                                                                                                                            atatGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacctttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                   >  1:425573‑M1/1‑143 (MQ=255)                                                                                                                              atGAATATCTTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTACTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaactggtgcacaatgcctgatgcgatgctt                                                                                                                 >  2:133584‑M1/1‑141 (MQ=255)                                                                                                                              atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                 >  1:57906‑M1/1‑141 (MQ=255)                                                                                                                               ggaATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCAGCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttg                                                                                                                <  1:284783‑M1/250‑112 (MQ=255)                                                                                                                                gAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGTGGAAAAGTACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCACCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggggccgcatccggcaattggtgcacaaagcctgatgcgatacttg                                                                                                                >  1:386023‑M1/1‑139 (MQ=255)                                                                                                                                 aaTATATTCCGGCGCTTAATGCCACGCCGGCACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGAGCCGTAGCGAAGCGCCTCGGCGGCGGTTTTCATCATCCGTGtattcgaaggtattagacatttgtgcctgatgcgccgcttgaagcgccttatcaggcctacaaccggtgcagcttccggcaattggagcaccatgcctgatgcgatgcttga                                                                                                               >  2:342353‑M1/1‑138 (MQ=255)                                                                                                                                 aaTATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAACTTGTTTCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACGGATCGCCGTCGCGAAGAGCCTCGGCGGCGGTGTTCATCAGCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcatgcctacaaccggtgccgcatccggcaattggtgcacaatgccggatgcgattcttgac                                                                                                              <  2:180192‑M1/251‑114 (MQ=255)                                                                                                                                    atatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgc                                                                                                            >  1:176551‑M1/1‑135 (MQ=255)                                                                                                                                      atTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatc                                                                                                         <  2:425573‑M1/251‑119 (MQ=255)                                                                                                                                        tCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcttaggcattagacatgtgttccttatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatct                                                                                                        <  1:201818‑M1/250‑120 (MQ=255)                                                                                                                                        tCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgtcgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatct                                                                                                        <  1:299187‑M1/250‑120 (MQ=255)                                                                                                                                            gctttgaatgccacGCCTGATCTTATCGAAATTATGGCTAATCTTGGTGCCGCGTATCGCTAAAGCTGAGTGGTAAATAACTTAGCGCCGTCGCGAAGAGCCTCGGCGGCTGTTTTCATCATCCGTGtagtcgtaggcattagacatttgttcctgatccgacgcttgactcgttttatcaggtctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgccgcatcttatca                                                                                                   <  1:433330‑M1/245‑125 (MQ=255)                                                                                                                                              gcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcagg                                                                                                 >  1:96867‑M1/1‑125 (MQ=255)                                                                                                                                                   aTGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctac                                                                                            <  1:336685‑M1/251‑132 (MQ=255)                                                                                                                                                     gCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaca                                                                                           >  1:113898‑M1/1‑118 (MQ=255)                                                                                                                                                        acgtcggaaCATATCGAAAGGATGGCTAATCTTGTGGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGAGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtcgtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgg                                                                                       <  2:98157‑M1/247‑137 (MQ=255)                                                                                                                                                        aCGCCGGACCATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggaattagacatttgcgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggtcattgtggcacaatgcctgttgcgatgcttgacgcatcttatcaggcctacaatgg                                                                                       >  1:321618‑M1/1‑115 (MQ=255)                                                                                                                                                            cGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtac                                                                                   >  1:112326‑M1/1‑111 (MQ=255)                                                                                                                                                             ggAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTAACGCTAACGCTGAGCGGTAAAGAACTGAGCGCGGTCGCGAAGCGCCTCGGCGGCGGTTTTCATCCTCCGTGtagtagtaggcattagacatttgtgcctgatgcggagcttgaggcgtcttatcaggcctacaaccggggccgcctccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtacc                                                                                  >  1:282447‑M1/1‑110 (MQ=255)                                                                                                                                                              gagcataTCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCTGTAAAGAGCGGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTGTTCATCATCCGTGtagtcgtaggcattagacattggtgcctgatgcgacgcttggcgcgtcttatcaggcttataaccggttccgcatccggcaattggtgcgcaaggccttaggcgatgcttgacgcatcttatcaggcctacaatgggtaccg                                                                                 <  1:254045‑M1/248‑143 (MQ=255)                                                                                                                                                               aaCATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAGGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcttgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccgg                                                                                >  1:334097‑M1/1‑108 (MQ=255)                                                                                                                                                                      cGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcacttggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggt                                                                          >  2:165294‑M1/1‑101 (MQ=255)                                                                                                                                                                          atgatgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggcc                                                                     <  1:3874‑M1/251‑155 (MQ=255)                                                                                                                                                                                    aTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgccggatgcgatgcttgacgcatcttatcaggcctaaaatgggtaccggatcggtaggccggataaagcg                                                           >  2:129034‑M1/1‑87 (MQ=255)                                                                                                                                                                                        tgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtt                                                         <  1:126815‑M1/249‑167 (MQ=255)                                                                                                                                                                                            gCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggct                                                                 >  2:352758‑M1/1‑79 (MQ=255)                                                                                                                                                                                            gCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggat                                                                 <  1:352758‑M1/237‑159 (MQ=255)                                                                                                                                                                                                 gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatc                                              >  1:81085‑M1/1‑74 (MQ=255)                                                                                                                                                                                                 gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgcc                                                   <  1:58709‑M1/246‑173 (MQ=255)                                                                                                                                                                                                 gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgcc                                                   >  2:58709‑M1/1‑74 (MQ=255)                                                                                                                                                                                                      gCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgagtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcctccggc                                          >  2:288258‑M1/1‑69 (MQ=255)                                                                                                                                                                                                       cTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggca                                         >  1:179991‑M1/1‑68 (MQ=255)                                                                                                                                                                                                               gAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagag                                <  1:7047‑M1/251‑192 (MQ=255)                                                                                                                                                                                                                   ggTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcacc                            >  1:88366‑M1/1‑56 (MQ=255)                                                                                                                                                                                                                                 cgcCGTCGCGAAGAGCCTCGGCGGCGGTTTGCATCAGCCGTTtagtcgtaggctgtagacatttgtgcctgatgcgacgcgtgacgcgtcttagcaggccttcaaccggggccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgagtaccggatcggtaggccggataatgcgtttacgccgcatccggcaagaatagagcaccagttaaccgaactt              <  2:404213‑M1/251‑210 (MQ=255)                                                                                                                                                                                                                                  gcCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaactta             <  1:260023‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                                    cGTCGCGAAGAGCCTCGGCGACGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttact           >  1:309285‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                       cgcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactct         <  1:428575‑M1/250‑215 (MQ=255)                                                                                                                                                                                                                                             gagCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgccca  <  1:129034‑M1/251‑222 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCG  >  NC_000913/498816‑499302

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