Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 566,770 C→A pseudogene (87/93 nt) renD → protein RenD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,7700CA100.0% 139.9 / NA 40pseudogene (87/93 nt)renDprotein RenD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/22);  total (18/22)

ACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGT  >  NC_000913/566528‑566904
                                                                                                                                                                                                                                                   |                                                                                                                                      
gcttgtccatatgactgcgAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAcc                                                                                                                                  >  2:575291‑M1/10‑247 (MQ=255)
gcttgtccatatgactgcgAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAcc                                                                                                                                  >  1:507132‑M1/10‑247 (MQ=255)
  ttgtccatatgactgCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagc                                                                                                                               >  1:318176‑M1/8‑245 (MQ=255)
        ataTGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatc                                                                                                                        <  1:215581‑M1/250‑13 (MQ=255)
        ataTGACTGCGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAAGCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatc                                                                                                                        <  2:722472‑M1/250‑13 (MQ=255)
         tATGACTGCGAGAATTAACCGTGGTGAGGCGATACCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAACGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCACCAGTAAccagcataaatc                                                                                                                        >  1:428487‑M1/1‑238 (MQ=255)
                 cgagaatTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcgg                                                                                                               <  1:333367‑M1/248‑22 (MQ=255)
                   agaATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCCTGGGCGGTAGCCCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggct                                                                                                             >  2:462025‑M1/2‑228 (MQ=255)
                   agaATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCGCAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGTAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGTCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggat                                                                                                             <  2:556800‑M1/250‑24 (MQ=255)
                   agaATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggct                                                                                                             <  2:314239‑M1/250‑24 (MQ=255)
                   agaATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagca                                                                                                                              <  2:518955‑M1/233‑7 (MQ=255)
                   agaATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagca                                                                                                                              >  1:518955‑M1/2‑228 (MQ=255)
                     aaTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcagaaatcaggcggtggctaa                                                                                                           <  2:274002‑M1/251‑26 (MQ=255)
                      aTTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaa                                                                                                           >  1:332450‑M1/1‑225 (MQ=255)
                        tAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaa                                                                                                         <  1:99370‑M1/250‑28 (MQ=255)
                           ccGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggc                                                                                                     >  1:480489‑M1/1‑220 (MQ=255)
                            cGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggc                                                                                                     >  2:332819‑M1/1‑219 (MQ=255)
                               ggtgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacg                                                                                                  <  2:524069‑M1/250‑35 (MQ=255)
                                 tgAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggc                                                                                               >  1:606793‑M1/1‑214 (MQ=255)
                                      cGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtc                                                                                          >  2:29129‑M1/1‑209 (MQ=255)
                                        atccctgaACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtctt                                                                                        <  2:344145‑M1/248‑45 (MQ=255)
                                           ccTGAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCGGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctg                                                                                     <  1:496107‑M1/251‑48 (MQ=255)
                                              gAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctgg                                                                                    <  1:2386‑M1/249‑49 (MQ=255)
                                              gAACCAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcaagggcaggtcttctggtt                                                                                  >  1:216927‑M1/1‑201 (MQ=255)
                                                  cAGTAAAACAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcg                                                                              >  1:331429‑M1/1‑197 (MQ=255)
                                                      aaaaCAACTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaagg                                                                          <  2:369015‑M1/251‑59 (MQ=255)
                                                           aacTTCCTGCCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaagg                                                                     <  2:216927‑M1/251‑64 (MQ=255)
                                                            acTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCATAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcagggcttctggttatcgaaggtcaggt                                                                    >  1:467128‑M1/1‑187 (MQ=255)
                                                            acTTCCTGTCATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaagttaaggt                                                                    <  1:180954‑M1/251‑65 (MQ=255)
                                                                         ggCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtg                                                       <  1:251491‑M1/251‑78 (MQ=255)
                                                                          gCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCGAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGACGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaagttctggcgaacggtgt                                                      <  2:638803‑M1/251‑79 (MQ=255)
                                                                                       aaaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgc                                         >  2:227532‑M1/1‑160 (MQ=255)
                                                                                        aaTCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgct                                        >  1:729090‑M1/1‑159 (MQ=255)
                                                                                              gCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccag                                   <  2:608027‑M1/250‑98 (MQ=255)
                                                                                                                         atcaaagctAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacg       <  1:283532‑M1/248‑126 (MQ=255)
                                                                                                                          tcaaagcTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacg       >  1:473267‑M1/3‑125 (MQ=255)
                                                                                                                           caaagCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacgaa     <  1:419521‑M1/250‑128 (MQ=255)
                                                                                                                           caaagCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacga      >  1:416149‑M1/2‑124 (MQ=255)
                                                                                                                            aaaGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacgaac    <  2:199011‑M1/251‑129 (MQ=255)
                                                                                                                              aGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAAccagcataaatcaggccgcggctaaaatggcacgggcaggtcttctggttatcgaaggtaaggtctggcgaacggtgtattaccggtttgctaccagggaagaacgggaaggaaagatgagcacgaacct  <  1:227532‑M1/251‑131 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                      
ACTGACCTGTATGACAGCGCGAATTAACCGTGGTGAGACGATACCTGAACCAGTAAAACAACTTCCTGTTATGGGCGGTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATTAAAGCTAAGTTCGGACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCATAAGAGCTTCTGTGCACAGGACGTTGCCGCGGTAACAGGCGCAACCGTAATCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGT  >  NC_000913/566528‑566904

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