Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,170 Δ1 bp coding (39/1818 nt) ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,1680C.99.5% 849.8 / ‑3.3 190coding (37/1818 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (45/144);  total (45/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 277 total aligned reads displayed.

TGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑‑G‑‑GGCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCT  >  NC_000913/3951085‑3951256
                                                                                         |                                                                                           
cactgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcgccac                                                                                         <  1:3375115‑M2/5‑3 (MQ=255)
cactgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcggggggtcgtgggccggggcgcgcccc                                                                                        >  1:5240097‑M2/89‑92 (MQ=255)
 actgggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcgccacc                                                                                        <  2:4782515‑M2/6‑4 (MQ=255)
    gggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCCCCGGa                                                                                   <  1:6960371‑M2/9‑1 (MQ=255)
     ggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtggtcgtgggctgggggggGCCCCCCGGa                                                                                   >  2:3122327‑M2/84‑93 (MQ=255)
     ggagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGaa                                                                                  <  1:5139639‑M2/10‑1 (MQ=255)
      gagtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggtgcGCCACCGGAAt                                                                                 <  1:7196342‑M2/11‑1 (MQ=255)
        gtaaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgggctgtggggcGCCCCCGggaaga                                                                                >  1:1283542‑M2/81‑91 (MQ=255)
         taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGAc                                                                             <  2:3033054‑M2/14‑1 (MQ=255)
         taaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGAc                                                                             <  2:5052314‑M2/14‑1 (MQ=255)
          agataaagtaggcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGAcc                                                                            <  2:2492763‑M2/15‑1 (MQ=255)
          aaataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcgtgcgctgtggggcGCCACCCGGATGAcc                                                                            >  2:5124531‑M2/79‑93 (MQ=255)
            ataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggctcgtgcgctgtggggcGCCCCCGGGATGAAcgc                                                                          >  2:2087177‑M2/77‑92 (MQ=255)
            ataaagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggggcgcgtgggctgtggggcGCCACCGGGCTGAACGa                                                                          >  2:1862358‑M2/77‑93 (MQ=255)
               aagtatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGc                                                                       <  2:1715031‑M2/20‑1 (MQ=255)
                  tatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGa                                                                    <  1:3350959‑M2/23‑1 (MQ=255)
                  tatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGa                                                                    <  2:4833670‑M2/23‑1 (MQ=255)
                  tatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGa                                                                    <  2:6038259‑M2/23‑1 (MQ=255)
                  tatgcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgccgtggcgcGGCCCCCGGAAGACCGCCGCCGa                                                                    >  2:399217‑M2/71‑93 (MQ=255)
                    cTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGGGCTGT‑‑‑GGCGGGCCCCCCGAATAAACGACCCCCATTTCg                                                               >  1:6813386/1‑93 (MQ=255)
                     gcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGAttt                                                                 <  1:7813406‑M2/26‑1 (MQ=255)
                     gcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGAttt                                                                 <  1:3224008‑M2/26‑1 (MQ=255)
                     gcctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGAttt                                                                 <  2:4599566‑M2/26‑1 (MQ=255)
                      cctaagtgccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTc                                                                <  1:5092584‑M2/27‑1 (MQ=255)
                      cctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTc                                                                <  1:6682216‑M2/27‑1 (MQ=255)
                      cctaagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTc                                                                <  1:1831671‑M2/27‑1 (MQ=255)
                      aaGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGTGCGCCGTGGGG‑‑GGCCCCCCGGAAGACCGGCGCCCCTTTCGGt                                                             >  2:7303874/1‑93 (MQ=255)
                       aGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTGCGCCGTG‑‑G‑‑GGCGCCACCCGGATGACCGACGCCGCTTTCggga                                                            >  2:6673763/1‑91 (MQ=255)
                        taagtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCgg                                                              <  2:5823191‑M2/29‑1 (MQ=255)
                        gTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCTTTCGCTGTG‑‑G‑‑GGCGCCACCGGCAAGGCCGACCCCGCGTTTGGTac                                                           >  2:527022/1‑92 (MQ=255)
                        gTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGT‑‑‑G‑‑GCGCGCCACCCGGAAGACCGACCCCCATTTCGGtta                                                           >  1:5785017/1‑91 (MQ=255)
                          agtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTa                                                            <  2:3405234‑M2/31‑1 (MQ=255)
                          aCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGGCGGGGGCTGT‑‑‑G‑‑GCGCGCCACCCGGATGAACGAAGCCGCTTTTTGTaatc                                                         >  2:1001260/1‑91 (MQ=255)
                           gtaccgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTaa                                                           <  2:1101704‑M2/32‑1 (MQ=255)
                            cGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCCGTG‑‑G‑‑GGCGCCACCCGGATGACCGACGCCCATTTCTGTGAGACg                                                       >  2:2981891/1‑93 (MQ=255)
                             gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGGCGT‑‑‑G‑‑GCGCGCCCCCGGGATGAGCGACGCCCCGTTTGGTAAGCCgc                                                      >  2:7713714/1‑92 (MQ=255)
                             gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCCCCTGCGCCGTG‑‑G‑‑GGCGCCACCGGAAAGACCGACGCCGCCTTTGGTAAGCCGa                                                      >  1:668936/1‑93 (MQ=255)
                             accgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGc                                                         <  2:3206998‑M2/34‑1 (MQ=255)
                             accgttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGc                                                         <  1:24153‑M2/34‑1 (MQ=255)
                             accgttccgccaccaccactcatggtcgtaatatggcgggggctcgtgcgctgtggcgcGCCACCCGGAAGACCGACGCCGACTTCGGTAAGc                                                         >  1:3455072‑M2/60‑93 (MQ=255)
                              ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑‑G‑‑GCGCGCCACCGGGATGAACGCCCCCCAATTCGGTAAACCggt                                                     >  1:5990716/1‑91 (MQ=255)
                              ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCCTGCGCTGTGGG‑‑GGGCCCCCGGAATGACCGCACCCCATTTTTGTAAACCGAt                                                     >  2:3309840/1‑93 (MQ=255)
                              ttCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCCGCGCCGT‑‑‑G‑‑GCGCGCCCCCGGCATGACCCACGCCGACTTCGGTAAGGCgct                                                     >  1:3968045/1‑91 (MQ=255)
                               tCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCCTGCGCCGTG‑‑G‑‑GGCGCCCCCGGGATGACCCGACGCCACTTTGGTAAGcccgtt                                                    >  2:2280274/1‑89 (MQ=255)
                                gttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                                                      <  1:2589270‑M2/37‑1 (MQ=255)
                                cttccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa                                                      <  2:4328033‑M2/37‑1 (MQ=255)
                                   ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                                   <  2:413120‑M2/40‑1 (MQ=255)
                                   ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                                   <  2:2144531‑M2/40‑1 (MQ=255)
                                   ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                                   <  2:3712906‑M2/40‑1 (MQ=255)
                                   ccgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                                                   <  1:2553403‑M2/40‑1 (MQ=255)
                                   ccaccaccaCTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTG‑‑G‑‑GGCGCCACCCGGATGAACGCCCCCCACTTTGGTAAGCCGGTTAttc                                                >  2:722569/1‑91 (MQ=255)
                                    cgccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTAt                                                  <  1:1626269‑M2/41‑1 (MQ=255)
                                     gccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATc                                                 <  1:6562893‑M2/42‑1 (MQ=255)
                                      ccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcg                                                <  1:1225639‑M2/43‑1 (MQ=255)
                                      ccaccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcg                                                <  1:7888568‑M2/43‑1 (MQ=255)
                                       caccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgc                                               <  2:122449‑M2/44‑1 (MQ=255)
                                         ccaccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCgg                                             <  1:1320888‑M2/46‑1 (MQ=255)
                                          caccactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                                            <  1:7174768‑M2/47‑1 (MQ=255)
                                          cacTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTG‑‑G‑‑GGCGCCCCCCGGATGACCGACGCCGATTTCGGGAAGGCGATTATCGCGGGtgg                                         >  2:5515193/1‑92 (MQ=255)
                                             tCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGT‑‑‑G‑‑GCGCGCCACCGGGATGAACGACGCCCATTTCGGGTAGGCGATTATCTCGGTTGTGaa                                      >  2:23522/1‑93 (MQ=255)
                                             cactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                                         <  1:1468338‑M2/50‑1 (MQ=255)
                                             cactcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                                         <  2:4779095‑M2/50‑1 (MQ=255)
                                                tcatggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGaa                                      <  1:5381206‑M2/53‑1 (MQ=255)
                                                  atggtcgtaatatggcggatgctcgtgcgctgtggcgcGACACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt                                    <  1:5041573‑M2/55‑1 (MQ=255)
                                                   tggtcgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc                                   <  1:4314685‑M2/56‑1 (MQ=255)
                                                       cgtaatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc                               <  1:7663020‑M2/60‑1 (MQ=255)
                                                       cgtaatatggcggggggtcgggcgccgtggggcGCCACCCGAATGACCGACGCCGATTTTGGTAAGCCGATTATCGCGGTTGTGAACTCgttt                               >  2:6958253‑M2/34‑92 (MQ=255)
                                                         taatatggcgggtgctcgtgggctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCCGTCAc                             >  2:6343380‑M2/32‑93 (MQ=255)
                                                         taatatggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAc                             <  1:4428691‑M2/62‑1 (MQ=255)
                                                            tatggcgggtgctcctgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATTGCGGTTGTGAACTCGTTTACCCa                          >  2:3089430‑M2/29‑93 (MQ=255)
                                                             atggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCaa                         <  2:7403341‑M2/66‑1 (MQ=255)
                                                              tggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAt                        <  2:2834105‑M2/67‑1 (MQ=255)
                                                              tggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCCCGGTTGTGAACTCGGTTACCCAAt                        >  2:6196827‑M2/27‑93 (MQ=255)
                                                               ggcgggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCCCCCAAtt                       >  1:2622375‑M2/26‑93 (MQ=255)
                                                                gcgggtgctcgtgcgctgtggcgcGCCACCGGAATCACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAAttt                      <  1:3807225‑M2/69‑1 (MQ=255)
                                                                   ggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTa                   >  2:5731333‑M2/22‑93 (MQ=255)
                                                                   ggtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTa                   <  2:581893‑M2/72‑1 (MQ=255)
                                                                    gtgctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAc                  >  2:3807809‑M2/21‑93 (MQ=255)
                                                                     tgctcgtgcgctgtggcgcGCCGCCAGAATGACCGACGCAGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc                 <  1:733463‑M2/74‑1 (MQ=255)
                                                                     tgctcgtgcgctgtggcgcGCCACCGGAATGACCGGCGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTAcc                 >  1:4021883‑M2/20‑93 (MQ=255)
                                                                      gctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                <  2:1379994‑M2/75‑1 (MQ=255)
                                                                      gctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                >  1:2208602‑M2/19‑93 (MQ=255)
                                                                      gctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                <  2:7199724‑M2/75‑1 (MQ=255)
                                                                      gctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                <  1:1164147‑M2/75‑1 (MQ=255)
                                                                      gctcgtgcgctgtggcgcACCCCCCGAATGCCCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                <  1:2046333‑M2/75‑1 (MQ=255)
                                                                       ctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCCATTTGTACCgg               >  1:408525‑M2/18‑93 (MQ=255)
                                                                       ctcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCgg               <  1:1150155‑M2/76‑1 (MQ=255)
                                                                        tcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCggg              >  2:6392042‑M2/17‑93 (MQ=255)
                                                                        tcgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCggg              >  1:4949025‑M2/17‑93 (MQ=255)
                                                                         cgtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGt             >  1:6437952‑M2/16‑93 (MQ=255)
                                                                          gtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGACCTCGTTCACCCAATTTGTACCGGGTc            >  1:2273306‑M2/15‑93 (MQ=255)
                                                                          gtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTc            >  2:6805789‑M2/15‑93 (MQ=255)
                                                                          gtgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTc            <  1:6251377‑M2/79‑1 (MQ=255)
                                                                           tgcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCa           <  1:4531503‑M2/80‑1 (MQ=255)
                                                                            gcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAc          <  2:275799‑M2/81‑1 (MQ=255)
                                                                            gcgctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCAc          <  2:2471960‑M2/81‑1 (MQ=255)
                                                                               ctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTNACCCAATTTGTACNGGGTCACGTc       >  2:135‑M2/10‑93 (MQ=255)
                                                                               ctgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTc       >  2:4116562‑M2/10‑93 (MQ=255)
                                                                                tgtggcgcGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTcc      >  1:5338369‑M2/9‑93 (MQ=255)
                                                                                     ggcGCCGCACCGGAATGACCGACGCCGATTTCGGTAAGCCTATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATc   <  1:5925400‑M2/90‑1 (MQ=255)
                                                                                         |                                                                                           
TGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑‑G‑‑GGCCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCT  >  NC_000913/3951085‑3951256

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