Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 813,480 A→C intergenic (+533/‑46) bioD → / → uvrB dethiobiotin synthetase/exision nuclease of nucleotide excision repair, DNA damage recognition component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913813,4800AC97.6% 696.8 / 2.2 207intergenic (+533/‑46)bioD/uvrBdethiobiotin synthetase/exision nuclease of nucleotide excision repair, DNA damage recognition component
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base C (107/95);  total (110/97)
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.76e-01

Only 100 of 213 total aligned reads displayed.

CATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATC  >  NC_000913/813384‑813574
                                                                                                |                                                                                              
    aTTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAAt                                                                                         >  2:3647726/1‑100 (MQ=255)
      tGACAGCGGAGTTTACGCTGTATCAGAAATATTATGTTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTa                                                                                       <  2:3915912/100‑1 (MQ=255)
      tGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTa                                                                                       >  2:5468777/1‑100 (MQ=255)
      tGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTa                                                                                       <  2:5401648/100‑1 (MQ=255)
      tGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTa                                                                                       <  2:2120026/100‑1 (MQ=255)
        aCAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACa                                                                                     <  1:1075614/100‑1 (MQ=255)
         cAGCGGAGTTTCCTCTGTACCAGAAATATAATGGTGAGGAATTGTTTTTTTATCAAGTATAATTTTTTGACATAATTAAGTACGACGCGTAAAATTACAt                                                                                    <  2:81001/100‑1 (MQ=255)
         cAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACAt                                                                                    >  2:950814/1‑100 (MQ=255)
         cAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACAt                                                                                    <  2:471880/100‑1 (MQ=255)
          aGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATa                                                                                   >  2:3137061/1‑100 (MQ=255)
          aGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATa                                                                                   <  2:611872/100‑1 (MQ=255)
          aGCGGAGTTTACGCTGTATCAGAAAGATTATGGTGATTAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACaca                                                                                   <  2:5098255/100‑3 (MQ=255)
           gCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATAc                                                                                  >  1:5076126/1‑100 (MQ=255)
           gCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATAc                                                                                  >  1:5237265/1‑100 (MQ=255)
           gCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATAc                                                                                  >  2:4963425/1‑100 (MQ=255)
            cGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATAcc                                                                                 <  1:3023443/100‑1 (MQ=255)
               aGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTgc                                                                              >  1:2140332/1‑100 (MQ=255)
                gTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTgcc                                                                             <  1:827153/100‑1 (MQ=255)
                gTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTgcc                                                                             >  2:4467612/1‑100 (MQ=255)
                gTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTgcc                                                                             <  1:3795764/100‑1 (MQ=255)
                     cGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTgcccgccc                                                                        >  1:1515611/1‑100 (MQ=255)
                     cGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTgcccgccc                                                                        >  2:44098/1‑100 (MQ=255)
                      gCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCa                                                                       <  1:636730/100‑1 (MQ=255)
                       cTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCaa                                                                      <  2:5347261/100‑1 (MQ=255)
                         gTATCAGAAATATTATGGTGATGAACTGTTTTTTTATTTATTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACt                                                                    >  2:1942005/1‑100 (MQ=255)
                         gTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACt                                                                    <  2:1925245/100‑1 (MQ=255)
                         gTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACt                                                                    <  1:1789801/100‑1 (MQ=255)
                          tATCAGAAATATTATGGTGATGAACTGTTTTTTTATTTATTTTTATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTc                                                                   >  2:1753728/1‑100 (MQ=255)
                          tATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTc                                                                   >  1:2768292/1‑100 (MQ=255)
                            tCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTccn                                                                 <  1:4438696/100‑2 (MQ=255)
                             cAGAAATATTATGGTGATGAACTGTTTTTTTTTTTTTTTTAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCtt                                                                >  2:4113656/1‑100 (MQ=255)
                             cAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCtt                                                                >  1:2177445/1‑100 (MQ=255)
                             cAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCtt                                                                >  2:3067618/1‑100 (MQ=255)
                              aGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTc                                                               >  2:2397471/1‑100 (MQ=255)
                               gAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCa                                                              <  1:4507367/100‑1 (MQ=255)
                               gAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCa                                                              <  2:4473282/100‑1 (MQ=255)
                                 aaTATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAgg                                                            <  1:4764345/100‑1 (MQ=255)
                                  atatTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCACCTCCTTCAGGt                                                           <  2:2234940/100‑1 (MQ=255)
                                  atatTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGt                                                           >  1:3458346/1‑100 (MQ=255)
                                  atatTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGt                                                           <  1:5070685/100‑1 (MQ=255)
                                    attatGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAg                                                         <  1:462305/100‑1 (MQ=255)
                                     ttatGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGc                                                        <  2:2739625/100‑1 (MQ=255)
                                      tatGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCg                                                       <  1:1630121/100‑1 (MQ=255)
                                      tatGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCg                                                       <  2:1714396/100‑1 (MQ=255)
                                       atGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGa                                                      >  1:1292312/1‑100 (MQ=255)
                                       atGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGCAGCgg                                                      >  1:4810603/1‑99 (MQ=255)
                                         ggTGATGAACTGTTTTTTTTTTTTTTTTTTTTTGTTGGCATAATTAATTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACt                                                    >  2:753481/1‑100 (MQ=255)
                                         ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTACAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACt                                                    >  2:3805537/1‑100 (MQ=255)
                                         ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACt                                                    >  2:2817759/1‑100 (MQ=255)
                                         ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACt                                                    <  1:3643195/100‑1 (MQ=255)
                                         ggTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACt                                                    >  1:1949138/1‑100 (MQ=255)
                                          ntgatgaACTGTTTTTTTTTTTTTTTTTATTTGTTGGAATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGAcac                                                   >  2:102948/2‑98 (MQ=255)
                                          gTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTc                                                   >  1:5530906/1‑100 (MQ=255)
                                          gTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTc                                                   >  1:4041317/1‑100 (MQ=255)
                                          gTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTc                                                   >  1:1204753/1‑100 (MQ=255)
                                           tgatgaACTGTTTTTTTTTTTTTTTTTTTTTGTTGGCATAATTAAGTACGACTCGTAAATTTACATACTTGCCCGCCCAACTCCTTCAGGTAGCGACTCa                                                  >  2:2925415/1‑100 (MQ=255)
                                           tgatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCa                                                  <  1:5388745/100‑1 (MQ=255)
                                           ngatgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCa                                                  >  1:2401763/2‑100 (MQ=255)
                                             atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATg                                                <  2:344447/100‑1 (MQ=255)
                                             atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATg                                                <  1:2286598/100‑1 (MQ=255)
                                             atgaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATg                                                <  2:2498265/100‑1 (MQ=255)
                                               gaACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAg                                              <  1:2844839/100‑1 (MQ=255)
                                                 aCTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTa                                            <  1:4388842/100‑1 (MQ=255)
                                                  cTGTTTTTTTTTTTATTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTaa                                           >  2:3040429/1‑100 (MQ=255)
                                                  cTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTaa                                           >  2:2010971/1‑100 (MQ=255)
                                                  cTGTTTTTTTATCCAGTAGACGTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTaa                                           <  2:3579383/100‑1 (MQ=255)
                                                    gTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAAc                                         >  2:1063879/1‑100 (MQ=255)
                                                       tttttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGt                                      >  1:4285087/1‑100 (MQ=255)
                                                          ttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCa                                   <  1:3218779/100‑1 (MQ=255)
                                                          ttATCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCa                                   <  1:4115424/100‑1 (MQ=255)
                                                            aTCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCaaa                                 >  1:1412302/1‑100 (MQ=255)
                                                             tCCAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAAc                                >  1:2631743/1‑100 (MQ=255)
                                                             tCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAAc                                >  2:3505511/1‑100 (MQ=255)
                                                              ccAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCACAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAGCt                               >  2:1387806/1‑100 (MQ=255)
                                                               cAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACtn                              <  1:5328192/100‑2 (MQ=255)
                                                               cAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTg                              <  1:2783749/100‑1 (MQ=255)
                                                               cAGTATAATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTg                              >  1:372547/1‑100 (MQ=255)
                                                                  tataATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAAt                           >  1:3847728/1‑100 (MQ=255)
                                                                  tataATTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAAt                           <  1:5473513/100‑1 (MQ=255)
                                                                      aTTTGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCg                       >  1:1104115/1‑100 (MQ=255)
                                                                       tttGTTGGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGc                      <  1:3999757/100‑1 (MQ=255)
                                                                          gttgGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCttt                   >  1:3287856/1‑100 (MQ=255)
                                                                          gttgGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCttt                   <  2:2158535/100‑1 (MQ=255)
                                                                          gttgGCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCttt                   <  1:5113198/100‑1 (MQ=255)
                                                                             ggCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTaa                >  1:2042157/1‑100 (MQ=255)
                                                                             ggCATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTaa                >  2:2062765/1‑100 (MQ=255)
                                                                               cATAATTAAGTACGACGCGTAAGATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAc              >  1:449680/1‑100 (MQ=255)
                                                                               cATAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAc              >  1:3392319/1‑100 (MQ=255)
                                                                                aTAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCATCTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAcc             >  1:4741188/1‑100 (MQ=255)
                                                                                aTAATTAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAAcc             <  2:4529754/100‑1 (MQ=255)
                                                                                     tAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTg        >  2:4613147/1‑100 (MQ=255)
                                                                                     tAAGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTg        <  2:1602341/100‑1 (MQ=255)
                                                                                       aGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGc      <  1:3539431/100‑1 (MQ=255)
                                                                                       aGTACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGc      <  1:1011636/100‑1 (MQ=255)
                                                                                         tACGACGCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGa    >  1:2867187/1‑100 (MQ=255)
                                                                                          acgacgCGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGAt   >  1:4935334/1‑100 (MQ=255)
                                                                                           cgacgcGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATc  <  2:1472064/100‑1 (MQ=255)
                                                                                           cgacgcGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATc  <  2:1335922/100‑1 (MQ=255)
                                                                                           cgacgcGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATc  <  2:4003887/100‑1 (MQ=255)
                                                                                           cgacgcGTAAAATTACATACCTGCCCGCCCAACTCCTGCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATc  <  1:1784740/100‑1 (MQ=255)
                                                                                                |                                                                                              
CATGATTGACAGCGGAGTTTACGCTGTATCAGAAATATTATGGTGATGAACTGTTTTTTTATCCAGTATAATTTGTTGGCATAATTAAGTACGACGAGTAAAATTACATACCTGCCCGCCCAACTCCTTCAGGTAGCGACTCATGAGTAAACCGTTCAAACTGAATTCCGCTTTTAAACCTTCTGGCGATC  >  NC_000913/813384‑813574

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