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% 142.1 / NA 43D305A (GAC→GCC) aesacetylesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/22);  total (21/22)

AAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTC  >  NC_000913/498821‑499301
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
aaaaagtctgagtccATATTCAGGTAAGGTTCGCGGGCTTTGCTGCTGATTGTCTTGAGTCGCTGGAAGTGATTGCCGAGCAAAATCGTGAGGTCTGCCTCGGTTCCGGCGAGAAAAAATATGAATAGATTCCGGCGCTTAATGCCACGCCGGAACATATCGAACTGATGGCTGCTCTTGTTGCCTCGTAGCTGTAAGGTTGAGCGTTAAAGAACGGATCGCCGTCGCGAAGTGCCTCGGCGGGGGTTTtc                                                                                                                                                                                                                                        <  2:303584/237‑1 (MQ=255)
    aGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTtcatca                                                                                                                                                                                                                                    <  2:158652/251‑1 (MQ=255)
                 atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCTCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGGGGCGGTTTTCATCATCCGTGtagtcgt                                                                                                                                                                                                                       >  1:352750‑M1/1‑244 (MQ=255)
                      ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggca                                                                                                                                                                                                                  <  1:144591‑M1/251‑13 (MQ=255)
                                      tttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGGCACGCCGGAACATATCGAAAAGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGGGGGTTTTATCATCCGTGtagtcgtagggctttgacatttgtgcct                                                                                                                                                                                                  >  1:217282‑M1/1‑223 (MQ=255)
                                         gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCCGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgat                                                                                                                                                                                               <  1:33256‑M1/251‑32 (MQ=255)
                                                                                                        gCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACAGATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacacccggtgccgcatccggcaattggtgcacaat                                                                                                                                <  2:142103‑M1/251‑95 (MQ=255)
                                                                                                           ggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCACCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGCAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgctgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcc                                                                                                                             >  1:77609‑M1/1‑154 (MQ=255)
                                                                                                             cggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttttgcctgatgcgacgcttgacgcgtcttatcattcctacaaccggttccgcatccggcaattggtgcacaatgcctg                                                                                                                           <  1:57443‑M1/251‑100 (MQ=255)
                                                                                                                 aaaaaaTATGAATATATTCCGGCGCTTAATACCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAAAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgc                                                                                                                       <  2:97‑M1/251‑104 (MQ=255)
                                                                                                                        atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                          <  2:446302‑M1/209‑69 (MQ=255)
                                                                                                                        atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                          >  1:446302‑M1/1‑141 (MQ=255)
                                                                                                                            atatatTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacg                                                                                                            >  1:101530‑M1/1‑137 (MQ=255)
                                                                                                                                atTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatc                                                                                                        <  1:309423‑M1/251‑119 (MQ=255)
                                                                                                                                atTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAGCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatgtgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatcctgcaattggtgcacaatgcctgatgcgatgcttgacgcatc                                                                                                        <  1:451415‑M1/251‑119 (MQ=255)
                                                                                                                                 acctgccgcttaatgCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatct                                                                                                       <  1:380030‑M1/245‑120 (MQ=255)
                                                                                                                                   ccGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGAGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                >  2:99318‑M1/1‑130 (MQ=255)
                                                                                                                                   ccGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGAGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                                <  1:99318‑M1/240‑111 (MQ=255)
                                                                                                                                      gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                                   >  1:337378‑M1/1‑127 (MQ=255)
                                                                                                                                        gcTTAATGACACGCCGGAACCTATCGAAATGATGGCTAATATTGTTGCCGCGTATCGCTACAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcggaggcattagacatttgtgcctgatgcgacgcttgaagcgtcttctccggcctacaaccggtgccgcatccggccattggtgcataatgccctatgcgatgcttgacgccacttatcagg                                                                                                >  2:285303‑M1/1‑125 (MQ=255)
                                                                                                                                           tAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGGGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGTCGGTTTTCATCATCCGTGgagtcggaggcattagacaggtgtgcctgatgcgacgcttgacgggtgtgatcaggcctacaaccggtaccgcatcctgcaatgtgtgcactatgcctgattcgatgcttgacgcagcttatcaggc                                                                                               <  2:403690‑M1/249‑128 (MQ=255)
                                                                                                                                            aaTGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCAGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggagcacactgcctgatgcgatgcttgacgcatcttatcaggcct                                                                                             >  1:223625‑M1/1‑121 (MQ=255)
                                                                                                                                                        gAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccg                                                                                <  2:337378‑M1/251‑143 (MQ=255)
                                                                                                                                                          aCATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtcggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggg                                                                              >  2:202060‑M1/1‑107 (MQ=255)
                                                                                                                                                              atCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcctcttatcaggcctacaatgggtaccggatcgg                                                                          >  1:70516‑M1/1‑103 (MQ=255)
                                                                                                                                                                cGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggt                                                                         >  1:452739‑M1/1‑101 (MQ=255)
                                                                                                                                                                         ggCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggat                                                                >  1:145050‑M1/1‑92 (MQ=255)
                                                                                                                                                                          gCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGGGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCCTCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgaagcttgacgcgtcttctcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggata                                                               >  1:406960‑M1/1‑91 (MQ=255)
                                                                                                                                                                          gCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccgggta                                                               >  2:16749‑M1/1‑91 (MQ=255)
                                                                                                                                                                                   gttgCCGCGTATAGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCAGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcaccttatcaggcctacaatgggtcccggctcggaaggccggctaaggcgattac                                                     >  2:166840‑M1/1‑82 (MQ=255)
                                                                                                                                                                                                     aGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcctcttatcaggcctacaatgggtaccggatcgg                                                                          >  2:358192‑M1/1‑64 (MQ=255)
                                                                                                                                                                                                     aGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcgg                                                                          <  1:358192‑M1/212‑149 (MQ=255)
                                                                                                                                                                                                               tAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCGGGGtagtcggaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcacgaatagagcacccg                         >  2:200481‑M1/1‑54 (MQ=255)
                                                                                                                                                                                                                aaaGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccag                         <  2:145050‑M1/250‑198 (MQ=255)
                                                                                                                                                                                                                   gAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgctttacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgtctgatgcgatgcgtgacgcatcttatcaggcctacaatggttaccggatcggtaggccggataatgcatttacgccgcatccggcaagaatagagcaccagttaa                     <  1:35109‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                                    aaCTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaac                    >  1:4264‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                                        gAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccg                  <  1:174320‑M1/249‑205 (MQ=255)
                                                                                                                                                                                                                           cgcCGGCGCGAGGAGCCTCGGCGGCGGTTTTCATCATCCGTGtggtcgtaggcattagacatttgtgcctggtgcgacgctgtacgcgtcgtatcaggcctacaaccgggtccgcatccggcaattggtgcgcaatgcctgatgcgatgcttgacgcagcttatcatgcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcttgcggcaagaagatagcaccatttaaccgaactt             <  2:62875‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                                              cGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttataaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaactttac          >  2:346388‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                              cGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcggaggcattagacatttgtgcctgatgcgacgcgtgacgcgtcttatcaggcctacatccggtggcgcatccggcaattggtgcacaatgcgtgatgcgatgcgtgacgcatcttatcaggcctacaatgggtgccggagcggtaggccggataaggcgtttagggcgtttccggcaagaatagagcaccagttaaccgaacttact          <  1:407323‑M1/251‑213 (MQ=255)
                                                                                                                                                                                                                               gTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactc         >  2:79124‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                                  gcgAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgc      <  1:121508‑M1/251‑217 (MQ=255)
                                                                                                                                                                                                                                      agagCCGCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgccc  <  2:377558‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
AAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTC  >  NC_000913/498821‑499301

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