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_000913499,0600TG100.0% 178.3 / NA 52D305A (GAC→GCC) aesacetylesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/27);  total (25/27)

GGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498817‑499305
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ggAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAAGGTGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGGCGGCGCGAAGAGCCTcggcggcggg                                                                                                                                                                                                                                                >  1:16985/1‑250 (MQ=255)
                  cATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtag                                                                                                                                                                                                                               <  2:170194‑M1/250‑4 (MQ=255)
                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                                         <  2:94554‑M1/250‑10 (MQ=255)
                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGAGTGCCGAGCAAAACCGTGAGGTCTGCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTGGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAACGCCGTCGCCAAGATCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                                         <  1:350058‑M1/250‑10 (MQ=255)
                              aTGTGCCCGGGCTTTGCTGCGGATTGCCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAAAATATTCCGGCGCTTAATGCCCCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAAATGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCCTCATCCGTGgagtcgcaggcatt                                                                                                                                                                                                                    >  1:476826‑M1/1‑235 (MQ=255)
                                                                cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgactcttgacgcgtctt                                                                                                                                                                                <  1:470936‑M1/251‑51 (MQ=255)
                                                                        gagaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcc                                                                                                                                                                        <  2:297699‑M1/251‑59 (MQ=255)
                                                                         agaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatgtgtgcctgatg                                                                                                                                                                                                  <  1:468615‑M1/224‑33 (MQ=255)
                                                                                   gCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtg                                                                                                                                                             <  1:261490‑M1/251‑70 (MQ=255)
                                                                                              aGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCCCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCAGCCGTGtagtcgtaggcattagacatttgggcctggtgcggcgcgtggcgcgtctgatcaggcgttcaaccggtgccgcatccggc                                                                                                                                                  <  2:581501‑M1/251‑81 (MQ=255)
                                                                                                     ccTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattgg                                                                                                                                            <  2:174595‑M1/250‑87 (MQ=255)
                                                                                                                gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcct                                                                                                                                >  1:446672‑M1/1‑153 (MQ=255)
                                                                                                                     aaaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgc                                                                                                                           >  1:354725‑M1/1‑148 (MQ=255)
                                                                                                                        aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                                        <  1:517277‑M1/251‑107 (MQ=255)
                                                                                                                        aaaTATGAATATATTCAGGCGCTTAATGCCACGCCGGCACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATACGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatccggcctacaacaggagccgcaccacgcaatgggtgcacaatgcctgctgcgct                                                                                                                        >  2:127466‑M1/1‑145 (MQ=255)
                                                                                                                          atatGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAAGGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGTTTTTCATCATCCGTGtagtcgtatgcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaacctgggccgcatccggtaattggtgcacaatgcctgatgcgatgc                                                                                                                      <  1:229937‑M1/251‑109 (MQ=255)
                                                                                                                            atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATAACGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCAGCAGCCGTGtagtcggaggcattagacatttgtgcctgatgcgacgcttgacgcgttttatcaggcctacaaccggggccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                    <  1:540030‑M1/251‑111 (MQ=255)
                                                                                                                                       ccGGCGCTTAATGCCACGCCGGAACATAGCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatctt                                                                                                          <  1:128155‑M1/250‑121 (MQ=255)
                                                                                                                                           cgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTACAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcctcttatc                                                                                                       >  1:410502‑M1/1‑126 (MQ=255)
                                                                                                                                                    ccACGCCGTAACATATCGAAATGATGGCTAATCTGGTGGCCGCGTATGACTAAAGCGGAGGGGTAAAGAACTGAGCGCCGGCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagccatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggttccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaat                                                                                            <  2:495514‑M1/251‑135 (MQ=255)
                                                                                                                                                         ccGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggggccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaaaatgggt                                                                                        >  2:455202‑M1/1‑112 (MQ=255)
                                                                                                                                                           ggAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtacc                                                                                     >  1:112562‑M1/1‑110 (MQ=255)
                                                                                                                                                               aagatcgaAATGAGGGCTAATCTTGTGGCCGCGTATCGCTAAAGCGGAGCGGTAAAGAACTGATCGCCGTCGCGAAGAGCATCGGCGGCGGTTTTCATCACCCGGTtagtcgtaggcactagacatttgtgcctgatgcgacgcttgacgcgtcgtagcaggcctgcaaccggtgccgcagccggcaattggtgcacaatgcctgatgctatgcgtgacgcatcttatcaggcagacaatgggtaccggat                                                                                 <  2:348251‑M1/248‑146 (MQ=255)
                                                                                                                                                                atatCGAAATGATGGCTAATCTTGTTGCCGCGTCTCGCTAAAGCTGAGCGGTAAAGCACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtcgtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaaaatgggtaccggatc                                                                                >  2:276479‑M1/1‑105 (MQ=255)
                                                                                                                                                                atatCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatc                                                                                <  1:102207‑M1/251‑147 (MQ=255)
                                                                                                                                                                   tCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggt                                                                             >  2:179588‑M1/1‑102 (MQ=255)
                                                                                                                                                                   tCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggaacggg                                                                             >  2:337886‑M1/1‑102 (MQ=255)
                                                                                                                                                                    cGAAATGATGGCTAATCGTGGTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggt                                                                             <  1:276479‑M1/250‑150 (MQ=255)
                                                                                                                                                                        atgatgGCTAATCTTGTTGCCGCGTATCGCTAACGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcaccatgcctgatgcgatgcttgacgcctcttatcaggcctacaattggtcccggctcggtaggcc                                                                        >  2:286419‑M1/1‑97 (MQ=255)
                                                                                                                                                                        atgatgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatccggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggcc                                                                        >  2:473690‑M1/1‑97 (MQ=255)
                                                                                                                                                                        atgatgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcagcttatcaggcctacaatgggtaccggatcggtaggcc                                                                        <  1:361945‑M1/251‑155 (MQ=255)
                                                                                                                                                                             ggCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggata                                                                   <  1:326902‑M1/251‑160 (MQ=255)
                                                                                                                                                                             ggCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggat                                                                    >  2:93722‑M1/1‑92 (MQ=255)
                                                                                                                                                                                      tgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatcaggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtt                                                            <  1:446284‑M1/249‑167 (MQ=255)
                                                                                                                                                                                       gttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcatgcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttac                                                         <  1:56900‑M1/251‑170 (MQ=255)
                                                                                                                                                                                       gttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttac                                                         <  1:179588‑M1/251‑170 (MQ=255)
                                                                                                                                                                                         tgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCAGCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcgtgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggggcacaaggcctgatgcgatgcttgacgcatcttagcagtcctacaatgggtaccg                                                                                    <  2:77134‑M1/222‑143 (MQ=255)
                                                                                                                                                                                         tgCCGCGTATCGCTAAAGCTGAGCGGTAAAAAACTGAGAGCCGCCGCGAAGTGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggccttagacatttgtgcctgatgcgacgcttgaagcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgctgcttgacgcatcttatcaggcctacaatgggtcccg                                                                                    >  1:77134‑M1/1‑80 (MQ=255)
                                                                                                                                                                                            cgcgTATCGCTAAAGCTGAACGGTCAAGAACTGAGCGCCGTCGAGAAGAGCCTCGGCGGCGGTTTTCATCATACGTGtagtcgcacgcattagacatttgtgcccgatgcgacgcttgacgcgtcttatcaggcctccaaccggtgccgaatccggcaattggtgcacaatgcctgctgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggacggatcaggcgtttatgccgc                                                    >  2:604039‑M1/1‑77 (MQ=255)
                                                                                                                                                                                               gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatc                                                 >  2:91143‑M1/1‑74 (MQ=255)
                                                                                                                                                                                                       aaaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggca                                            >  1:86513‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                           gtgAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcagggctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcgtgacgcatcttatcaggcctacaatgggtaccggatcgttagtccggataaggcgtttacgccgcatccggcaagaatag                                     <  1:604039‑M1/250‑190 (MQ=255)
                                                                                                                                                                                                                cGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcac                                <  1:443582‑M1/251‑195 (MQ=255)
                                                                                                                                                                                                                 ggTAAAGAACTGAGCGCCGTCGCGACGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcggaggcattagacatttgtgcctgatgccacgcttgacgcgtcttatcaggcctacaaccggtgccgcctccggcaattggtgcacaatgcctgatgcgatgcttgccgcatcttatcaggcctccaatgggtaccggatcggtaggccggataaggcgtttactccgcatccggcaagactagagcacc                               >  2:103903‑M1/1‑56 (MQ=255)
                                                                                                                                                                                                                   tAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcacc                               >  2:333689‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                          cTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccg                      >  2:464008‑M1/1‑47 (MQ=255)
                                                                                                                                                                                                                                    gcgcgaAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactct            <  1:91143‑M1/250‑215 (MQ=255)
                                                                                                                                                                                                                                     cgcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctg           >  1:376628‑M1/1‑36 (MQ=255)
                                                                                                                                                                                                                                      gcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtaggcgtaggcattagacatttgtgcctgatgcgacggttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgc          <  1:418689‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                       cgAAGAGCCTCGGCGGCGGCTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcg         >  1:156254‑M1/1‑34 (MQ=255)
                                                                                                                                                                                                                                         aaGAGCCTCGGCGGCGGTTTTAATCATACGTGtagtcgtaggcattagacacttgcgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggagcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggacaaggcgtttacgccgcatccggcaagaataaagcaccagttaaccgaacttaatctgcgcc       >  2:405247‑M1/1‑32 (MQ=255)
                                                                                                                                                                                                                                          agagCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgccc      >  1:220763‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                                              ccTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagccatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtatgccggcaaaggcgtttacgccgcatccggcaagcatagagcacccgttaaccgaacttactcttgcccccact  >  2:167153‑M1/1‑27 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498817‑499305

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