Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,606 Δ1 bp 26.5% pseudogene (65/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,6060A.26.5% 454.2 / 170.5 170pseudogene (65/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base A (59/65);  new base . (24/21);  total (83/87)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 171 total aligned reads displayed.

CTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGC  >  NC_000913/3951508‑3951707
                                                                                                   |                                                                                                     
ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                                                                      >  1:2983933/1‑100 (MQ=255)
 tGTTACCAGCATTACAGCAGCCGTTAAACAAGGACCGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACg                                                                                                     <  2:60659/99‑1 (MQ=255)
    tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTa                                                                                                  >  1:1573853/1‑100 (MQ=255)
    tACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTa                                                                                                  <  1:2484433/100‑1 (MQ=255)
      ccAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTa                                                                                              >  1:2401209/1‑100 (MQ=255)
      ccAGCATTACAGCAGCCGTTAAATCAATGACTAGCAGAAACACTACGCGCAGCTGCGTGATGAACATTCCTGGAGTTACGACCATCCCGGTGACGCTATc                                                                                                >  2:304604/1‑100 (MQ=255)
       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACCGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCt                                                                                               <  1:476396/100‑1 (MQ=255)
        aGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTa                                                                                              >  1:2202965/1‑100 (MQ=255)
         gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                                                                             >  1:802680/1‑100 (MQ=255)
         gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                                                                             >  2:665004/1‑100 (MQ=255)
         gCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTAcgc                                                                                           >  2:1024020/1‑100 (MQ=255)
         gCATTACAGCAGCCGTTAAACAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                                                                            <  1:1431329/100‑1 (MQ=255)
          cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAcg                                                                                            >  1:1473787/1‑100 (MQ=255)
          cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTAcgc                                                                                           >  1:315863/1‑99 (MQ=255)
             tACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACgcgc                                                                                         >  1:1212509/1‑100 (MQ=255)
              aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCgt                                                                                      >  2:1434848/1‑100 (MQ=255)
               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCgtt                                                                                     >  2:1737412/1‑100 (MQ=255)
               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCgtt                                                                                     >  2:2047217/1‑100 (MQ=255)
                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCgttgt                                                                                   >  1:738704/1‑100 (MQ=255)
                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                                                                  <  1:866228/100‑1 (MQ=255)
                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCgttgttgt                                                                                <  2:2202949/100‑1 (MQ=255)
                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCgttgttgt                                                                                >  1:2265404/1‑100 (MQ=255)
                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCgttgttgt                                                                                <  2:316273/100‑1 (MQ=255)
                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                                                                 <  2:1762917/100‑1 (MQ=255)
                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttg                                                                                 <  1:1278036/100‑1 (MQ=255)
                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCgttgttgtt                                                                               <  1:3180849/100‑1 (MQ=255)
                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCgttgttgtt                                                                               <  2:258687/100‑1 (MQ=255)
                     ccGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCgttgttgtt                                                                               <  2:1664553/100‑1 (MQ=255)
                       gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCG‑‑GACGCTATCTACGCGCCGTTGTTGTTaa                                                                             <  1:1979176/100‑1 (MQ=255)
                         tAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                                                                              <  1:2814133/100‑1 (MQ=255)
                          aaaTCAATGACTGACAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAaac                                                                          <  1:2454740/100‑1 (MQ=255)
                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                                            >  2:882937/1‑100 (MQ=255)
                           aaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                                                                            >  1:229289/1‑100 (MQ=255)
                           aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAaaca                                                                         >  1:1866697/1‑100 (MQ=255)
                             tCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                                                                          <  1:1936992/100‑1 (MQ=255)
                             tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaaca                                                                         <  1:1434529/100‑1 (MQ=255)
                             tCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAaacaac                                                                       <  2:1083384/100‑1 (MQ=255)
                              cAATGACGGTCAGCAACACTGCAACCGAGTTGCGTAATGAACATTCCTAGTGTTACTACCATCCCTGTGAAGCTATCTTCTCGCCGTTGTTGTTaaaa                                                                           <  2:933625/98‑1 (MQ=255)
                               aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTTTaacaac                                                                       <  1:1030180/100‑1 (MQ=255)
                               aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTg                                                                     >  1:2950255/1‑100 (MQ=255)
                                 tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                                                                     >  2:74558/1‑100 (MQ=255)
                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                                                                  <  2:2804664/100‑1 (MQ=255)
                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                                                                  <  2:619686/100‑1 (MQ=255)
                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCGACGCGCCGTTGTTGTTAAAACAACTGTCg                                                                  >  2:1945547/1‑100 (MQ=255)
                                      ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                                              >  2:2269793/1‑100 (MQ=255)
                                      ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc                                                              >  1:2724044/1‑100 (MQ=255)
                                       gcagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCg                                                             <  1:3188070/100‑1 (MQ=255)
                                         agcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCg                                                             <  1:1173684/100‑1 (MQ=255)
                                          gcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa                                                          <  1:433439/100‑1 (MQ=255)
                                             acacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc                                                       <  2:2401193/100‑1 (MQ=255)
                                              cacTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt                                                      >  1:989132/1‑100 (MQ=255)
                                                 tGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg                                                   <  2:2539952/100‑1 (MQ=255)
                                                 tGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg                                                   <  1:85118/100‑1 (MQ=255)
                                                  gcgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCgg                                                  <  2:2983906/100‑1 (MQ=255)
                                                   cgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg                                                   <  1:394175/100‑1 (MQ=255)
                                                      cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGa                                                 <  2:2628666/100‑1 (MQ=255)
                                                      cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGAtt                                               <  2:1095553/100‑1 (MQ=255)
                                                      cgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGAtt                                               >  2:1651365/1‑100 (MQ=255)
                                                          gCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGcgt                                           >  1:2603723/1‑100 (MQ=255)
                                                              cGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCATAAACCTGCGGATTGcgtcg                                         <  2:1990243/100‑1 (MQ=255)
                                                              cGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                                       <  1:130130/100‑1 (MQ=255)
                                                              cGTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                                       >  2:1201326/1‑100 (MQ=255)
                                                               gTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                                      >  1:3108748/1‑100 (MQ=255)
                                                               gTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTg                                       <  1:1067776/99‑1 (MQ=255)
                                                               gTGATGAACATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                                      <  1:1096933/100‑1 (MQ=255)
                                                                tgatgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAc                                     <  2:802678/100‑1 (MQ=255)
                                                                 gatgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGa                                      >  2:509139/1‑100 (MQ=255)
                                                                  atgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACca                                   <  1:1603432/100‑1 (MQ=255)
                                                                   tgaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGAcc                                    <  1:1009056/100‑1 (MQ=255)
                                                                    gaccATTCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                                  <  1:2816163/96‑1 (MQ=255)
                                                                    gaACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACca                                   >  1:2899133/1‑100 (MQ=255)
                                                                     aaCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                                  <  1:2582401/100‑1 (MQ=255)
                                                                     aaCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcac                                  <  1:1980844/100‑1 (MQ=255)
                                                                      aCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACcaca                                 <  1:1644495/100‑1 (MQ=255)
                                                                      aCATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGa                               <  1:1661661/100‑1 (MQ=255)
                                                                       cATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAt                              >  1:1824336/1‑100 (MQ=255)
                                                                          tCCTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgtg                           <  1:2269809/100‑1 (MQ=255)
                                                                           ccTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg                          <  1:1771779/100‑1 (MQ=255)
                                                                           ccTGGCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTgg                          <  2:3251363/100‑1 (MQ=255)
                                                                            cTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgtg                           >  2:1606248/1‑100 (MQ=255)
                                                                            cTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGAtgtg                           >  2:1605296/1‑100 (MQ=255)
                                                                            cTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg                         >  1:3062979/1‑100 (MQ=255)
                                                                            cTGGCGTTACGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTggg                         <  1:2337318/100‑1 (MQ=255)
                                                                              ggCGTTACGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgc                       <  2:369723/100‑1 (MQ=255)
                                                                                   tACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGgcagc                    >  1:1957503/1‑100 (MQ=255)
                                                                                    aCGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAcc                 <  1:1248895/100‑1 (MQ=255)
                                                                                    aCGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAcc                 <  1:510494/100‑1 (MQ=255)
                                                                                    aCGACCATCCCGGTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAcc                 <  2:1473779/100‑1 (MQ=255)
                                                                                     cGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCAc                  >  2:2892152/1‑100 (MQ=255)
                                                                                     cGACCATCCCGGT‑‑CGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa                >  1:500484/1‑100 (MQ=255)
                                                                                       tccATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa                >  2:2635445/2‑100 (MQ=255)
                                                                                       aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa                >  2:374708/1‑100 (MQ=255)
                                                                                       aCCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCa                <  1:3005093/100‑1 (MQ=255)
                                                                                       aCCATCCCGGT‑‑CACTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGa              >  2:1437927/1‑100 (MQ=255)
                                                                                           tCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAtg            >  2:2914216/1‑100 (MQ=255)
                                                                                           tCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAtg            <  1:1813573/100‑1 (MQ=255)
                                                                                               ggTGACGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCt      >  2:782988/1‑100 (MQ=255)
                                                                                                gTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGAtgtg          <  1:1119704/96‑1 (MQ=255)
                                                                                                 tGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGTATTGCGTCGTGACCACAGATGTGGCGCAGCACCCGAtgttgtct      >  2:3175484/1‑95 (MQ=255)
                                                                                                   aCGCTATCTA‑‑CGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTgcgc  <  1:1762241/100‑1 (MQ=255)
                                                                                                   |                                                                                                     
CTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCGGATTGCGTCGTGACCACAGATGTGGGGCAGCACCAGATGTGGGCTGCGC  >  NC_000913/3951508‑3951707

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