Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,193 C→T 42.1% P1317S (CCG→TCG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,1930CT42.1% 94.7 / 168.2 158P1317S (CCG→TCG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (49/42);  new base T (31/35);  total (80/78)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 160 total aligned reads displayed.

GCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCA  >  NC_000913/4185094‑4185292
                                                                                                   |                                                                                                   
   tACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGATGGCAACCATCAGATGGAGCCgg                                                                                                  <  1:2276457/100‑1 (MQ=255)
    acacCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCggg                                                                                                 <  1:3399413/100‑1 (MQ=255)
    acacCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCggg                                                                                                 <  1:1918154/100‑1 (MQ=255)
     cacCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCggac                                                                                                >  1:108809/1‑98 (MQ=255)
     cacCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGc                                                                                                >  2:2927206/1‑100 (MQ=255)
      acCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCa                                                                                               <  2:1008792/100‑1 (MQ=255)
        ccTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATg                                                                                             >  2:2341942/1‑100 (MQ=255)
         cTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGc                                                                                            >  2:1875037/1‑100 (MQ=255)
         cTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGc                                                                                            <  1:2989985/100‑1 (MQ=255)
         cTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGc                                                                                            >  2:3407404/1‑100 (MQ=255)
          tGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGcc                                                                                           >  1:222448/1‑100 (MQ=255)
          tGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGcc                                                                                           <  2:534661/100‑1 (MQ=255)
            cAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAg                                                                                         <  2:487288/100‑1 (MQ=255)
            cAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCAGGCATGCCAg                                                                                         <  2:1383571/100‑1 (MQ=255)
             aGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGa                                                                                        >  2:3619965/1‑100 (MQ=255)
                aaaTGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATc                                                                                     <  2:1068928/100‑1 (MQ=255)
                 aaTGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATcc                                                                                    <  1:3815539/100‑1 (MQ=255)
                     cTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCCAGATGTATAAAAACCTCGTGGACGGCAACCAGCAGATGGAGTCGGGCATGCTAGAATCCTTCa                                                                                >  2:174108/1‑100 (MQ=255)
                     cTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCa                                                                                >  1:3344652/1‑100 (MQ=255)
                     cTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCa                                                                                >  2:808533/1‑100 (MQ=255)
                      tCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCaa                                                                               >  1:1219086/1‑100 (MQ=255)
                      tCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCaa                                                                               >  2:763763/1‑100 (MQ=255)
                       cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAAc                                                                              <  1:3605127/100‑1 (MQ=255)
                       cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                                                                              <  1:798357/100‑1 (MQ=255)
                       cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                                                                              >  1:2410376/1‑100 (MQ=255)
                       cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAAc                                                                              <  1:3348090/100‑1 (MQ=255)
                        aCCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACg                                                                             <  2:2108262/100‑1 (MQ=255)
                         ccGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGt                                                                            >  1:2769030/1‑100 (MQ=255)
                           gTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTAt                                                                          >  1:2298219/1‑100 (MQ=255)
                           gTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAt                                                                          >  1:3546535/1‑100 (MQ=255)
                             tAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCATGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTAttg                                                                        <  2:2852303/100‑1 (MQ=255)
                              aaGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCATTCAACg                                                                             <  2:2007026/94‑1 (MQ=255)
                              aaGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                                                                       >  1:2494490/1‑100 (MQ=255)
                                gTCTGATGACGTGAACGGTCTTACCAAGGTGTATAAAAACATCGTGGACGGAAACCATCAGATGGATCCGGGCATGCTAGAATCCTGGCACGTAtt                                                                         >  2:3355118/1‑96 (MQ=255)
                                gTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgttg                                                                     >  1:842346/1‑100 (MQ=255)
                                  cTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGaa                                                                   >  2:129782/1‑100 (MQ=255)
                                   tgatgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGaaa                                                                  <  2:473170/100‑1 (MQ=255)
                                   tgatgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGaca                                                                  >  2:2512148/1‑98 (MQ=255)
                                    gatgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAag                                                                 <  1:3374139/100‑1 (MQ=255)
                                     atgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAaga                                                                <  2:2196202/100‑1 (MQ=255)
                                      tgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAagag                                                               <  1:2482718/100‑1 (MQ=255)
                                      tgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAagag                                                               <  1:1081345/100‑1 (MQ=255)
                                       gaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAgaga                                                              <  2:2791192/100‑1 (MQ=255)
                                        aCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAt                                                             >  2:1722865/1‑100 (MQ=255)
                                         cGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAtt                                                            >  1:2499234/1‑100 (MQ=255)
                                          gTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttc                                                           <  2:3605548/100‑1 (MQ=255)
                                          gTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttc                                                           <  1:1368574/100‑1 (MQ=255)
                                           tGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttcg                                                          <  1:3407430/100‑1 (MQ=255)
                                             aaCGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCCTGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttcgtt                                                        >  2:2686113/1‑100 (MQ=255)
                                                 gTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCt                                                    >  2:834216/1‑100 (MQ=255)
                                                   cGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTgg                                                  <  1:683598/100‑1 (MQ=255)
                                                    gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg                                                 <  2:1034147/100‑1 (MQ=255)
                                                    gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATGCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg                                                 <  1:2695093/100‑1 (MQ=255)
                                                    gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg                                                 <  2:2396775/100‑1 (MQ=255)
                                                           gATGTATAAAAACATTGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAAc                                          >  1:1290476/1‑100 (MQ=255)
                                                            aTGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAGAGAGATTCGTTCGCTGGGTATCAACa                                         <  1:2330419/100‑1 (MQ=255)
                                                             tGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACAt                                        >  2:1454910/1‑100 (MQ=255)
                                                              gTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATc                                       >  1:3151873/1‑100 (MQ=255)
                                                              gTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATc                                       >  1:846385/1‑100 (MQ=255)
                                                              gTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATc                                       <  2:3331418/100‑1 (MQ=255)
                                                                 taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGaa                                    >  1:2416017/1‑100 (MQ=255)
                                                                   aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCCTGCCAGAATCCCTCAACGTATTGTTGCAAGAGATTCGTTCGCTGGGTATCAACATCGAACt                                  >  2:1708879/1‑100 (MQ=255)
                                                                   aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACt                                  <  1:2225766/100‑1 (MQ=255)
                                                                   aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACt                                  >  1:2177987/1‑100 (MQ=255)
                                                                    aaaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTg                                 <  2:394778/100‑1 (MQ=255)
                                                                     aaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                <  1:142867/100‑1 (MQ=255)
                                                                     aaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                <  2:719308/100‑1 (MQ=255)
                                                                      aCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGa                               <  1:2149911/100‑1 (MQ=255)
                                                                        aTCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAg                             <  1:472117/100‑1 (MQ=255)
                                                                          cGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGAc                           >  1:1848509/1‑100 (MQ=255)
                                                                          cGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGAc                           <  1:960596/100‑1 (MQ=255)
                                                                           gTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACg                          <  1:2988150/100‑1 (MQ=255)
                                                                           gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACg                          <  1:2675951/100‑1 (MQ=255)
                                                                           gTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACg                          <  2:881299/100‑1 (MQ=255)
                                                                            tGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGa                         >  1:1521727/1‑100 (MQ=255)
                                                                             ggACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAg                        >  2:1306656/1‑100 (MQ=255)
                                                                               aCGGCAACCATCAGATCGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGTGTATCAACATCGAACTGGAAGACGAGTa                      <  2:1558635/100‑1 (MQ=255)
                                                                                cGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTaa                     >  2:1435181/1‑100 (MQ=255)
                                                                                cGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTaa                     >  1:2756806/1‑100 (MQ=255)
                                                                                 ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAt                    <  1:348267/100‑1 (MQ=255)
                                                                                  gCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAtt                   >  1:1215374/1‑100 (MQ=255)
                                                                                   cAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtc                  >  2:241153/1‑100 (MQ=255)
                                                                                       cATCAGATTGAGGAGGGCATGCCATAATGCTGCAACGTATTGTTGAAAGGGATTCGTTCGCTGGGTATCAACATCGAGCTGGAAGACGAGTAGTTCTCGc              <  2:1118914/100‑1 (MQ=255)
                                                                                       cATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGc              >  2:2775863/1‑100 (MQ=255)
                                                                                        aTCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCt             >  2:2729856/1‑100 (MQ=255)
                                                                                        aTCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCt             >  1:1902346/1‑100 (MQ=255)
                                                                                         tCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTc            >  1:929035/1‑100 (MQ=255)
                                                                                         tCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTc            <  1:1520645/100‑1 (MQ=255)
                                                                                          cAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCa           <  1:3155125/100‑1 (MQ=255)
                                                                                            gATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCaaa         >  1:322226/1‑100 (MQ=255)
                                                                                             attgAGCCGAGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAAc        <  2:2689912/97‑1 (MQ=255)
                                                                                             aTGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCaaan        <  2:348910/100‑2 (MQ=255)
                                                                                              tGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTTGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTAGCTCAAACa       >  2:2521898/1‑100 (MQ=255)
                                                                                                gAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg     >  2:1527790/1‑100 (MQ=255)
                                                                                                gAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg     <  2:3546508/100‑1 (MQ=255)
                                                                                                 aGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGt    >  2:2064568/1‑100 (MQ=255)
                                                                                                 aGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGt    >  1:3466487/1‑100 (MQ=255)
                                                                                                  gCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc   >  1:516868/1‑100 (MQ=255)
                                                                                                  gCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc   >  1:2267316/1‑100 (MQ=255)
                                                                                                   ccGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCa  <  1:2927227/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCA  >  NC_000913/4185094‑4185292

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