Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,016,180 A→G I448T (ATA→ACA)  b0955 ← putative ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,016,1800AG82.6% 577.3 / 119.8 271I448T (ATA→ACA) b0955putative ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base A (20/27);  new base G (114/109);  total (134/137)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

Only 100 of 271 total aligned reads displayed.

AAACCACCGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCT  >  NC_000913/1016090‑1016270
                                                                                          |                                                                                          
  accaccGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAg                                                                                        >  1:5831814/1‑93 (MQ=255)
  accaccGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAg                                                                                        >  2:2839675/1‑93 (MQ=255)
    caccGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAggg                                                                                      >  1:4179647/1‑93 (MQ=255)
    caccGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAggg                                                                                      >  2:3953924/1‑93 (MQ=255)
    caccGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAggg                                                                                      >  1:1854291/1‑93 (MQ=255)
     accgaccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGgc                                                                                     <  2:6549542/93‑1 (MQ=255)
     accgaccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGgc                                                                                     <  1:4676786/93‑1 (MQ=255)
     accgaccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGgc                                                                                     <  1:4944004/93‑1 (MQ=255)
       cgaccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGgcgc                                                                                   >  2:1883136/1‑93 (MQ=255)
        gaccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCa                                                                                  >  2:2015920/1‑93 (MQ=255)
         accgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCag                                                                                 <  1:1926698/93‑1 (MQ=255)
          ccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCaga                                                                                <  1:5912680/93‑1 (MQ=255)
          ccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCaga                                                                                >  2:1076558/1‑93 (MQ=255)
          ccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCaga                                                                                <  2:5945206/93‑1 (MQ=255)
          ccgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCaga                                                                                <  1:7049007/93‑1 (MQ=255)
           cgGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCagag                                                                               >  2:5354671/1‑93 (MQ=255)
            ggCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGt                                                                              <  2:6322622/93‑1 (MQ=255)
            ggCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGt                                                                              >  2:1067068/1‑93 (MQ=255)
             gCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGtt                                                                             >  2:861462/1‑93 (MQ=255)
             gCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGtt                                                                             >  1:168482/1‑93 (MQ=255)
             gCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGtt                                                                             >  2:323813/1‑93 (MQ=255)
              cTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGATGGCGCTTTTCAGGGCGCAGAGTTc                                                                            <  2:7359815/93‑1 (MQ=255)
              cTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTc                                                                            <  2:2248214/93‑1 (MQ=255)
              cTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTc                                                                            <  1:6910906/93‑1 (MQ=255)
               tGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCa                                                                           >  2:4714873/1‑93 (MQ=255)
                gggCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAg                                                                          >  1:4052684/1‑93 (MQ=255)
                 ggCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGc                                                                         <  1:1468718/93‑1 (MQ=255)
                    gcgACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCAt                                                                      <  2:2966609/93‑1 (MQ=255)
                     cgACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATc                                                                     <  1:840730/93‑1 (MQ=255)
                      gACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCg                                                                    >  2:1172879/1‑93 (MQ=255)
                      gACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCg                                                                    >  1:6455387/1‑93 (MQ=255)
                       aCCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGa                                                                   <  1:2522944/93‑1 (MQ=255)
                       aCCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGa                                                                   <  2:6010770/93‑1 (MQ=255)
                       aCCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGa                                                                   >  2:6106405/1‑93 (MQ=255)
                        ccGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAg                                                                  <  1:5461062/93‑1 (MQ=255)
                         cGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAgg                                                                 <  1:5595572/93‑1 (MQ=255)
                         cGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAgg                                                                 <  1:1398759/93‑1 (MQ=255)
                         cGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAgg                                                                 >  2:4266472/1‑93 (MQ=255)
                         cGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAgg                                                                 >  2:4134930/1‑93 (MQ=255)
                         cGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAgg                                                                 >  2:6616320/1‑93 (MQ=255)
                          gAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGc                                                                <  1:772572/93‑1 (MQ=255)
                           aaCTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCa                                                               <  2:4999589/93‑1 (MQ=255)
                            aCTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCAc                                                              >  2:2380843/1‑93 (MQ=255)
                            aCTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCAc                                                              >  2:4482560/1‑93 (MQ=255)
                              tGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTa                                                            >  1:7120954/1‑93 (MQ=255)
                              tGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTa                                                            >  1:3014510/1‑93 (MQ=255)
                              tGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTa                                                            <  2:6219367/93‑1 (MQ=255)
                              tGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTa                                                            >  1:6425453/1‑93 (MQ=255)
                               gATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTAt                                                           <  1:5285912/93‑1 (MQ=255)
                               gATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTAt                                                           >  2:3694364/1‑93 (MQ=255)
                                 tCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTAtct                                                         <  2:7499286/93‑1 (MQ=255)
                                 tCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTAtct                                                         >  1:2000015/1‑93 (MQ=255)
                                   gACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTcc                                                       >  1:445303/1‑93 (MQ=255)
                                     cTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCAt                                                     >  1:6530508/1‑93 (MQ=255)
                                          ccTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACt                                                <  2:4078137/93‑1 (MQ=255)
                                          ccTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACt                                                <  1:5159504/93‑1 (MQ=255)
                                          ccTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACt                                                <  2:5161878/93‑1 (MQ=255)
                                           cTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACtt                                               >  2:1705773/1‑93 (MQ=255)
                                              tgATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCAc                                            >  2:4218237/1‑93 (MQ=255)
                                               gATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACt                                           >  2:7240715/1‑93 (MQ=255)
                                                aTAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTg                                          >  1:2399991/1‑93 (MQ=255)
                                                   gCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTAt                                       <  1:6962615/93‑1 (MQ=255)
                                                   gCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTAt                                       >  1:1988446/1‑93 (MQ=255)
                                                   gCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTAt                                       <  1:83433/93‑1 (MQ=255)
                                                     gATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGa                                     >  1:400297/1‑93 (MQ=255)
                                                      aTACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGAc                                    <  2:2576692/93‑1 (MQ=255)
                                                       tACTCTGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACt                                   <  1:500718/93‑1 (MQ=255)
                                                       tACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACt                                   <  1:2226797/93‑1 (MQ=255)
                                                         ctctGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGc                                 <  1:6191096/93‑1 (MQ=255)
                                                          tctGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCt                                >  1:4120612/1‑93 (MQ=255)
                                                            tGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCa                              >  2:1574548/1‑93 (MQ=255)
                                                            tGATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCa                              >  2:6883585/1‑93 (MQ=255)
                                                            tGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCtnn                              >  1:3419826/1‑91 (MQ=255)
                                                             gATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCaa                             >  2:3403579/1‑93 (MQ=255)
                                                             gATTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCaa                             >  2:2267566/1‑93 (MQ=255)
                                                              aTTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCaaa                            <  2:5555151/93‑1 (MQ=255)
                                                              aTTCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCaaa                            >  1:511429/1‑93 (MQ=255)
                                                               ttCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAAt                           <  1:1049244/93‑1 (MQ=255)
                                                                tCACCGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATg                          <  1:3692863/93‑1 (MQ=255)
                                                                   ccGGCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCa                       <  1:273722/93‑1 (MQ=255)
                                                                     ggCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGc                     <  1:966222/93‑1 (MQ=255)
                                                                     ggCACATCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGc                     <  2:4954922/93‑1 (MQ=255)
                                                                        acaTCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGAt                  >  2:3449395/1‑93 (MQ=255)
                                                                        acaTCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGAt                  <  2:5822526/93‑1 (MQ=255)
                                                                          aTCGGCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGc                <  2:2722/93‑1 (MQ=255)
                                                                             ggCGAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTga             <  1:1070957/93‑1 (MQ=255)
                                                                                gAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGaa          >  1:3175947/1‑93 (MQ=255)
                                                                                gAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGaa          <  1:7339723/93‑1 (MQ=255)
                                                                                gAGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGaa          <  2:6760724/93‑1 (MQ=255)
                                                                                gAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGaa          <  2:5402593/93‑1 (MQ=255)
                                                                                 aGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAg         <  1:1750319/93‑1 (MQ=255)
                                                                                 aGGGCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAg         <  2:7071626/93‑1 (MQ=255)
                                                                                  gggCGCTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAgg        >  1:6237506/1‑93 (MQ=255)
                                                                                     cgcTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGa     <  1:4396469/93‑1 (MQ=255)
                                                                                      gcTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGAt    >  2:1390660/1‑93 (MQ=255)
                                                                                      gcTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGAt    <  2:1440452/93‑1 (MQ=255)
                                                                                       cTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATc   >  1:3910267/1‑93 (MQ=255)
                                                                                       cTTGTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATc   >  1:5659618/1‑93 (MQ=255)
                                                                                       cTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATc   >  2:3795643/1‑93 (MQ=255)
                                                                                        ttATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCt  >  1:3516743/1‑93 (MQ=255)
                                                                                          |                                                                                          
AAACCACCGACCGGCTGGGCGCGACCGAACTGATCGACTGAACCTGTGATAGCGATACTCTGATTCACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCT  >  NC_000913/1016090‑1016270

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