Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,189,448 C→T 79.0% Q1367* (CAG→TAG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,4480CT79.0% 1041.2 / 284.9 488Q1367* (CAG→TAG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (49/53);  new base T (173/211);  total (224/264)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 500 total aligned reads displayed.

GTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTA  >  NC_000913/4189300‑4189595
                                                                                                                                                    |                                                                                                                                                   
gTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTAccacca                                                                                                                                                    <  1:5532092/150‑1 (MQ=255)
gTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTAccacca                                                                                                                                                    <  2:2116044/150‑1 (MQ=255)
   cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGa                                                                                                                                                 >  1:263613/1‑150 (MQ=255)
    ttCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGAt                                                                                                                                                >  2:2731040/1‑150 (MQ=255)
    ttCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGGTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGAt                                                                                                                                                >  2:1753914/1‑150 (MQ=255)
       aTCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGt                                                                                                                                             <  1:4042381/150‑1 (MQ=255)
          tcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATg                                                                                                                                          >  2:1276813/1‑150 (MQ=255)
          tcCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATg                                                                                                                                          >  1:6348347/1‑150 (MQ=255)
           ccGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACGAGGATCGTATGc                                                                                                                                         >  2:4191130/1‑150 (MQ=255)
            cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCg                                                                                                                                        >  2:4554847/1‑150 (MQ=255)
              cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGACCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTc                                                                                                                                      >  2:6497517/1‑150 (MQ=255)
              cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACTGGTTACGCGTACCACTAGGATCGTATGCGTc                                                                                                                                      <  1:4201276/150‑1 (MQ=255)
               ggCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCg                                                                                                                                     <  2:3403941/150‑1 (MQ=255)
                  aTCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCg                                                                                                                                  >  2:3581097/1‑150 (MQ=255)
                    cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtg                                                                                                                                <  2:1670018/150‑1 (MQ=255)
                        ccAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGtgctgc                                                                                                                            <  2:4076544/150‑1 (MQ=255)
                          aGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCgg                                                                                                                          <  2:773859/150‑1 (MQ=255)
                            gagaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTCGGATCGTATGCGTCGCCGTGCTGCGGGt                                                                                                                        >  1:1397115/1‑150 (MQ=255)
                            gagaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGt                                                                                                                        <  2:5312716/150‑1 (MQ=255)
                            gagaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGt                                                                                                                        <  2:5560524/150‑1 (MQ=255)
                            gagaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGt                                                                                                                        <  1:877414/150‑1 (MQ=255)
                             agaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTg                                                                                                                       <  2:3280544/150‑1 (MQ=255)
                             agaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTg                                                                                                                       <  2:4697342/150‑1 (MQ=255)
                              gaCCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGa                                                                                                                      <  1:6281348/150‑1 (MQ=255)
                              gaCCACTCGCGTGCTGACCGAAGCAGCCGGTGCGGGCAAACGCGACGAACTGCGCGGCCTGACCGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTCCGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGa                                                                                                                      <  2:5236520/150‑1 (MQ=255)
                                 cACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGACCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGc                                                                                                                   >  1:2229079/1‑150 (MQ=255)
                                      gcgTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCgg                                                                                                              >  2:4478721/1‑150 (MQ=255)
                                      gcgTGCTGACCGAAGCAGCCGCTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCgg                                                                                                              >  2:492283/1‑150 (MQ=255)
                                      gcgTGCTGACCGAAGCAGCCGCTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCgg                                                                                                              >  2:642807/1‑150 (MQ=255)
                                        gTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                            >  1:3699310/1‑150 (MQ=255)
                                        gTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                            >  2:2225532/1‑150 (MQ=255)
                                        gTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                            <  2:5968300/150‑1 (MQ=255)
                                         tGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTg                                                                                                           <  1:1754023/150‑1 (MQ=255)
                                           cTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCa                                                                                                         <  1:806473/150‑1 (MQ=255)
                                           cTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCa                                                                                                         <  1:6084060/150‑1 (MQ=255)
                                             gACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCAcc                                                                                                       >  2:3629828/1‑150 (MQ=255)
                                              aCCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCg                                                                                                      >  1:5070852/1‑150 (MQ=255)
                                                cGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCa                                                                                                    <  1:6461722/150‑1 (MQ=255)
                                                  aaGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAgg                                                                                                  >  2:5495716/1‑150 (MQ=255)
                                                    gcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTg                                                                                                >  1:1638262/1‑150 (MQ=255)
                                                     cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGa                                                                                               >  2:4838268/1‑150 (MQ=255)
                                                     cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGa                                                                                               <  2:4313793/150‑1 (MQ=255)
                                                         cGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGc                                                                                           <  2:3959970/150‑1 (MQ=255)
                                                           ttGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCag                                                                                         >  1:2545776/1‑150 (MQ=255)
                                                               gggCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaaga                                                                                     <  2:5448982/150‑1 (MQ=255)
                                                                 gCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACg                                                                                   >  1:2315724/1‑150 (MQ=255)
                                                                 gCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACg                                                                                   >  1:4495870/1‑150 (MQ=255)
                                                                   aaaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCAGCTGCACCGCAGGTGACTGCAGAAGACGCa                                                                                 >  2:2462000/1‑150 (MQ=255)
                                                                    aaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCAt                                                                                <  1:2872960/150‑1 (MQ=255)
                                                                       gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTg                                                                             >  2:5757748/1‑150 (MQ=255)
                                                                       gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTg                                                                             >  1:3009096/1‑150 (MQ=255)
                                                                       gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTg                                                                             >  1:956917/1‑150 (MQ=255)
                                                                       gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTg                                                                             <  1:4712943/150‑1 (MQ=255)
                                                                        cgacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGc                                                                            <  1:1357388/150‑1 (MQ=255)
                                                                         gacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGcc                                                                           >  2:3159909/1‑150 (MQ=255)
                                                                         gacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGcc                                                                           >  1:2776957/1‑150 (MQ=255)
                                                                           cgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAg                                                                         <  2:5392855/150‑1 (MQ=255)
                                                                               cTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTg                                                                     >  1:3820464/1‑150 (MQ=255)
                                                                                tGCGCGGCCTGAAAGAGAACGTTATCGTGTGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTgg                                                                    <  1:197686/150‑1 (MQ=255)
                                                                                tGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTgg                                                                    >  1:5466140/1‑150 (MQ=255)
                                                                                  cgcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCa                                                                  <  1:1021685/150‑1 (MQ=255)
                                                                                   gcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAg                                                                 >  1:6386329/1‑150 (MQ=255)
                                                                                     ggCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGaa                                                               <  1:4295072/150‑1 (MQ=255)
                                                                                        cTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTAACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGAGACTCCAGAAGACGCATCTGCCCGCCTGGCAGAActg                                                            >  1:3394196/1‑150 (MQ=255)
                                                                                         tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgc                                                           <  2:5661096/150‑1 (MQ=255)
                                                                                          gAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAActgct                                                          >  2:1377776/1‑150 (MQ=255)
                                                                                             agagAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTAACGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGaa                                                       >  1:4567825/1‑150 (MQ=255)
                                                                                                 aaCGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCa                                                   >  2:752916/1‑150 (MQ=255)
                                                                                                   cGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAgg                                                 >  1:533305/1‑150 (MQ=255)
                                                                                                       aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTg                                             >  1:3122154/1‑150 (MQ=255)
                                                                                                        tCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGTGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTgg                                            >  2:63751/1‑150 (MQ=255)
                                                                                                         cGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTggg                                           <  2:6007590/150‑1 (MQ=255)
                                                                                                          gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                          <  2:4425559/150‑1 (MQ=255)
                                                                                                              gtcgtcTGATCCCGGCAGGTACCGGTTACGCGGACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGtt                                      <  1:4779618/150‑1 (MQ=255)
                                                                                                               tcgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTc                                     >  1:5842691/1‑150 (MQ=255)
                                                                                                               tcgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTc                                     >  1:6582429/1‑150 (MQ=255)
                                                                                                                cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                    <  1:3043109/150‑1 (MQ=255)
                                                                                                                cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt                                    <  2:1917388/150‑1 (MQ=255)
                                                                                                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                  >  2:1091938/1‑150 (MQ=255)
                                                                                                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                  <  1:1166010/150‑1 (MQ=255)
                                                                                                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                  <  2:3820464/150‑1 (MQ=255)
                                                                                                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                  >  1:4990752/1‑150 (MQ=255)
                                                                                                                  tcTGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGAGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                  >  1:1234394/1‑150 (MQ=255)
                                                                                                                    tGATCCCGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATa                                <  2:1497401/150‑1 (MQ=255)
                                                                                                                      aTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAAc                              >  2:1447168/1‑150 (MQ=255)
                                                                                                                        cccGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGa                            <  1:3190942/150‑1 (MQ=255)
                                                                                                                        cccGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGa                            <  2:3060387/150‑1 (MQ=255)
                                                                                                                         ccGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg                           >  2:5039147/1‑150 (MQ=255)
                                                                                                                         ccGGCAGGTACCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg                           <  1:5210733/150‑1 (MQ=255)
                                                                                                                               ggTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCg                     <  2:2365219/150‑1 (MQ=255)
                                                                                                                                  aCCGGTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTa                  <  2:2885878/150‑1 (MQ=255)
                                                                                                                                      gTTACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATcc              <  1:3307391/150‑1 (MQ=255)
                                                                                                                                       ttACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCg             >  1:3477887/1‑150 (MQ=255)
                                                                                                                                       ttACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCg             <  1:650887/150‑1 (MQ=255)
                                                                                                                                        tACGCGTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGc            <  2:5851942/150‑1 (MQ=255)
                                                                                                                                         aCGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGACGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa           <  2:611865/150‑1 (MQ=255)
                                                                                                                                           gcgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCaaa         >  2:370130/1‑150 (MQ=255)
                                                                                                                                           gcgTACCACTAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCaaa         >  1:2386314/1‑150 (MQ=255)
                                                                                                                                           gcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCaaa         >  1:2105920/1‑150 (MQ=255)
                                                                                                                                                 cactaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt   >  2:2628295/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTA  >  NC_000913/4189300‑4189595

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