Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,185,211 T→C 33.3% F1323L (TTC→CTC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,2110TC33.3% 156.9 / 137.2 153F1323L (TTC→CTC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (52/50);  new base C (26/25);  total (78/75)
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.43e-01

Only 100 of 156 total aligned reads displayed.

TGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAC  >  NC_000913/4185113‑4185310
                                                                                                  |                                                                                                   
tGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCtt                                                                                                    <  1:2393955/100‑1 (MQ=255)
 gCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTc                                                                                                   <  1:3364422/100‑1 (MQ=255)
  ntCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCa                                                                                                  >  1:348912/2‑100 (MQ=255)
  cTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCa                                                                                                  >  2:808533/1‑100 (MQ=255)
  cTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCa                                                                                                  >  1:3344652/1‑100 (MQ=255)
   tCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCaa                                                                                                 >  1:1219086/1‑100 (MQ=255)
    cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATACAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                                                                                                >  2:774955/1‑100 (MQ=255)
    cACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAAc                                                                                                <  1:798357/100‑1 (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)
           aaGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCATTCAACg                                                                                               <  2:2007026/94‑1 (MQ=255)
           aaGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                                                                                         >  2:1858746/1‑100 (MQ=255)
           aaGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgt                                                                                         >  1:2494490/1‑100 (MQ=255)
            aGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTAttgtt                                                                                        <  1:2429377/100‑1 (MQ=255)
             gTCTGATGACGTGAACGGTCTTACCAAGGTGTATAAAAACATCGTGGACGGAAACCATCAGATGGATCCGGGCATGCTAGAATCCTGGCACGTAtt                                                                                           >  2:3355118/1‑96 (MQ=255)
             gTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTAttgttg                                                                                       >  1:842346/1‑100 (MQ=255)
               cTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGaa                                                                                     >  2:129782/1‑100 (MQ=255)
                tgatgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGaca                                                                                    >  2:2512148/1‑98 (MQ=255)
                 gatgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAag                                                                                   <  1:3374139/100‑1 (MQ=255)
                  atgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAaga                                                                                  <  2:2196202/100‑1 (MQ=255)
                   tgaCGTGAACGTTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAagag                                                                                 <  1:3409158/100‑1 (MQ=255)
                   tgaCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAagag                                                                                 <  1:382716/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)
                      cGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAtt                                                                              <  2:347667/100‑1 (MQ=255)
                       gTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttc                                                                             <  2:2769010/100‑1 (MQ=255)
                       gTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttc                                                                             <  1:1368574/100‑1 (MQ=255)
                        tGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcg                                                                            >  2:104844/1‑100 (MQ=255)
                        tGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttcg                                                                            <  1:3407430/100‑1 (MQ=255)
                          aaCGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCCTGCCAGAATCCCTCAACGTATTGTTGAAAGAGAttcgtt                                                                          >  2:2686113/1‑100 (MQ=255)
                            cGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGAttcgttcg                                                                        >  2:3352474/1‑100 (MQ=255)
                              gTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCt                                                                      >  2:834216/1‑100 (MQ=255)
                                 gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg                                                                   >  1:814074/1‑100 (MQ=255)
                                 gTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATGCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTggg                                                                   <  1:2695093/100‑1 (MQ=255)
                                   aCCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTa                                                                 <  1:1822121/100‑1 (MQ=255)
                                      aaGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCa                                                              >  2:509370/1‑100 (MQ=255)
                                      aaGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCa                                                              <  2:1425249/100‑1 (MQ=255)
                                      aaGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCa                                                              >  2:1713939/1‑100 (MQ=255)
                                        gATGTATAAAAACATTGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAAc                                                            >  1:1290476/1‑100 (MQ=255)
                                        gATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAAc                                                            <  2:2399694/100‑1 (MQ=255)
                                         aTGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAGAGAGATTCGTTCGCTGGGTATCAACa                                                           <  1:2330419/100‑1 (MQ=255)
                                         aTGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACa                                                           <  2:455802/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)
                                              taAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGaa                                                      >  1:2416017/1‑100 (MQ=255)
                                                aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACt                                                    <  1:2225766/100‑1 (MQ=255)
                                                aaaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACt                                                    >  1:2177987/1‑100 (MQ=255)
                                                 aaaCATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTg                                                   <  1:3661101/100‑1 (MQ=255)
                                                  aaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                                  <  2:719308/100‑1 (MQ=255)
                                                  aaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                                  <  2:2805826/100‑1 (MQ=255)
                                                  aaCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                                  <  1:142867/100‑1 (MQ=255)
                                                  aaCAGCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTgg                                                  <  1:2786621/100‑1 (MQ=255)
                                                   aCATCGTGGACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGa                                                 <  1:2149911/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)
                                                          ggACGGCAACCATCAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAg                                          >  2:874742/1‑100 (MQ=255)
                                                            aCGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTa                                        >  1:2987174/1‑100 (MQ=255)
                                                            aCGGCAACCATCAGATCGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGTGTATCAACATCGAACTGGAAGACGAGTa                                        <  2:1558635/100‑1 (MQ=255)
                                                              ggCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAt                                      <  1:348267/100‑1 (MQ=255)
                                                               gCAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAAtt                                     >  1:1215374/1‑100 (MQ=255)
                                                                cAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtc                                    >  2:241153/1‑100 (MQ=255)
                                                                cAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtc                                    <  2:162599/100‑1 (MQ=255)
                                                                cAACCATCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtc                                    >  1:2498352/1‑100 (MQ=255)
                                                                 aaCCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATtc                                    <  2:1032625/99‑1 (MQ=255)
                                                                    cATCAGATTGAGGAGGGCATGCCATAATGCTGCAACGTATTGTTGAAAGGGATTCGTTCGCTGGGTATCAACATCGAGCTGGAAGACGAGTAGTTCTCGc                                <  2:1118914/100‑1 (MQ=255)
                                                                    cATCAGATGGAGCCGGGCATGCCAGACTCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCACCATCGAACTGGAAGACGAGTAATTCTCGc                                <  2:1512343/100‑1 (MQ=255)
                                                                     aTCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCt                               >  1:2703192/1‑100 (MQ=255)
                                                                     aTCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCt                               >  1:1902346/1‑100 (MQ=255)
                                                                     aTCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCt                               >  2:2729856/1‑100 (MQ=255)
                                                                      tCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTc                              <  1:1520645/100‑1 (MQ=255)
                                                                      tCAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTc                              >  1:929035/1‑100 (MQ=255)
                                                                       cAGATGGAGTCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCa                             <  1:3155125/100‑1 (MQ=255)
                                                                       cAGATGGAGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCa                             <  1:252645/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)
                                                                             gAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAgg                       <  2:3546508/100‑1 (MQ=255)
                                                                              aGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGt                      >  2:2064568/1‑100 (MQ=255)
                                                                              aGCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGt                      >  1:1869753/1‑100 (MQ=255)
                                                                               gtcGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc                     >  2:3071197/3‑100 (MQ=255)
                                                                               gCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc                     >  1:516868/1‑100 (MQ=255)
                                                                               gCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc                     <  2:2247383/100‑1 (MQ=255)
                                                                               gCCGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc                     >  1:2267316/1‑100 (MQ=255)
                                                                                 cGGGCATGCCAGAATCCCTCAACGTATTGTTGAAAGATATTCGTTCGCTGTGTATCAACACCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTc                     <  2:2268598/98‑1 (MQ=255)
                                                                                        gCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTc            <  1:1156373/100‑1 (MQ=255)
                                                                                         ccAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCg           >  2:126843/1‑100 (MQ=255)
                                                                                          cAGAATCCTTCAGCGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCgg          >  2:3422607/1‑100 (MQ=255)
                                                                                          cAGAATCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCgg          >  1:1282374/1‑100 (MQ=255)
                                                                                             aaTCCCTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGtt       <  2:289648/100‑1 (MQ=255)
                                                                                                 cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaaa   >  2:3808802/1‑100 (MQ=255)
                                                                                                 cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaaa   >  2:288447/1‑100 (MQ=255)
                                                                                                 cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTaaaa   <  2:886503/100‑1 (MQ=255)
                                                                                                 cTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGCTaaaa   >  1:542979/1‑100 (MQ=255)
                                                                                                  ttCAACGTATTGTTGAAAGAGATTCGTTCGCTGAGCAGCACCATCGAACTGGAAGACGAGTAATTCTCGATCAAACAGGTCACTGCTGCCGGGCTACAAc  >  2:383280/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGAATCCTTCAACGTATTGTTGAAAGAGATTCGTTCGCTGGGTATCAACATCGAACTGGAAGACGAGTAATTCTCGCTCAAACAGGTCACTGCTGTCGGGTTAAAAC  >  NC_000913/4185113‑4185310

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: