Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,603 A→C E418D (GAA→GAC rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,6030AC98.4% 680.1 / 0.2 192E418D (GAA→GACrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base C (97/92);  total (97/95)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 196 total aligned reads displayed.

AACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAC  >  NC_000913/4186506‑4186699
                                                                                                 |                                                                                                
aaCTGCGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAcca                                                                                                >  1:5442629/1‑100 (MQ=255)
aaCTGCGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGCCGTTTGGGATATCCTGGACGAAGTTATCCGCGAcca                                                                                                >  2:4933831/1‑100 (MQ=255)
  cTGCGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGaccacc                                                                                              <  1:3573742/100‑1 (MQ=255)
    gCGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCg                                                                                            <  1:4681200/100‑1 (MQ=255)
    gCGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGCTGAGCGCGAAGAAGCTGTCCCTTGGGATATCCTCGACGAAGTTATCCGCGACCACCCg                                                                                            >  2:3401706/1‑100 (MQ=255)
     cGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCgg                                                                                           <  1:3138382/100‑1 (MQ=255)
     cGTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCgg                                                                                           <  1:521080/100‑1 (MQ=255)
      gTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGt                                                                                          >  1:4499299/1‑100 (MQ=255)
      gTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGt                                                                                          <  2:1439033/100‑1 (MQ=255)
      gTGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGt                                                                                          <  1:3852384/100‑1 (MQ=255)
       tGGTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGCTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTa                                                                                         >  1:2495008/1‑100 (MQ=255)
        ggTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTAc                                                                                        <  2:3736825/100‑1 (MQ=255)
        ggTCTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTAc                                                                                        <  2:289426/100‑1 (MQ=255)
           cTTGCTACCACCATGAAAGCTGCGCAGAAAATGTTTGAGCGCGAAGACGCCGTCGTTTGGGATATCCTGGACGCAGTTATCCGCGACCACCCGGTActgc                                                                                     >  2:4101188/1‑100 (MQ=255)
           cTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                     <  2:3352985/100‑1 (MQ=255)
           cTTGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                     >  1:1493465/1‑100 (MQ=255)
             tGCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgctg                                                                                   <  1:5516357/100‑1 (MQ=255)
              gCTACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGa                                                                                  >  2:2080205/1‑100 (MQ=255)
                tACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAc                                                                                <  2:1266132/100‑1 (MQ=255)
                tACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAc                                                                                <  1:3815474/100‑1 (MQ=255)
                tACCACCATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAc                                                                                >  1:4929293/1‑100 (MQ=255)
                 accaccATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAcc                                                                               >  1:3619927/1‑100 (MQ=255)
                 accaccATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAcc                                                                               >  1:1659237/1‑100 (MQ=255)
                 accaccATGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAcc                                                                               >  1:4642229/1‑100 (MQ=255)
                  ccaccaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCg                                                                              <  1:798144/100‑1 (MQ=255)
                  ccaccaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCg                                                                              >  2:579575/1‑100 (MQ=255)
                  ccaccaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCg                                                                              >  1:1019794/1‑100 (MQ=255)
                   caccaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGCAGAAGCTGTCGTTTGGGATATCCTGGACGACGTTATCCGCGACCACCCGGTACTGCTGAACCGt                                                                             >  2:4275271/1‑100 (MQ=255)
                    accaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTg                                                                            >  1:4108126/1‑100 (MQ=255)
                    accaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTg                                                                            <  1:4599877/100‑1 (MQ=255)
                     ccaTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGc                                                                           <  1:5137601/100‑1 (MQ=255)
                      caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCa                                                                          <  1:875447/100‑1 (MQ=255)
                      caTGAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCa                                                                          >  1:2153614/1‑100 (MQ=255)
                       atgaaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCAc                                                                         >  2:148872/4‑100 (MQ=255)
                           aaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTCGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGAc                                                                     <  1:112668/100‑1 (MQ=255)
                           aaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGGTTGGGATATCCTGGCCGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGAc                                                                     <  1:122130/100‑1 (MQ=255)
                            aGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGAct                                                                    <  2:102999/100‑1 (MQ=255)
                             gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActc                                                                   <  1:1341620/100‑1 (MQ=255)
                             gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActc                                                                   <  2:1214187/100‑1 (MQ=255)
                              cTGCTAAGAAAATGGTTGAGCTCGAAGAGGCTGTCGTTTGGGAGATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActat                                                                  >  2:2229691/1‑98 (MQ=255)
                              cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActct                                                                  >  2:4026651/1‑100 (MQ=255)
                              cTGCGAAGAAAATCGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActct                                                                  >  2:4754228/1‑100 (MQ=255)
                               tGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                 <  2:3315729/100‑1 (MQ=255)
                                 cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGAACTGCTGAACCGTGCACCGACTCTGCa                                                               >  1:4827340/1‑100 (MQ=255)
                                  gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                              <  1:2550737/100‑1 (MQ=255)
                                  gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATAACCTGGACCAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                              >  1:5277363/1‑100 (MQ=255)
                                    agaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                            >  2:159776/1‑100 (MQ=255)
                                    agaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                            >  2:362855/1‑100 (MQ=255)
                                    agaaAATCGTTGAGCGCTAAGAAGCTGTCGTTCGGGATATCCTGTACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                            <  1:5115831/100‑1 (MQ=255)
                                     gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                           >  1:731149/1‑100 (MQ=255)
                                     gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                           >  2:3360841/1‑100 (MQ=255)
                                     gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                           >  1:4898078/1‑100 (MQ=255)
                                     gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGt                                                           >  1:2565757/1‑100 (MQ=255)
                                       aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                         >  1:533157/1‑100 (MQ=255)
                                        aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                        >  2:646365/1‑100 (MQ=255)
                                        aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                        <  2:859634/100‑1 (MQ=255)
                                          tGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                      <  1:2095523/100‑1 (MQ=255)
                                          tGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                      <  2:5369585/100‑1 (MQ=255)
                                          tGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                      <  1:2879396/100‑1 (MQ=255)
                                            gTTTAGCGCGAAGTAGCTGTCGTTTGGGATCCCCTGGACGAAGTTATTCGCGACTACGCGGTACTGCCGAACCGTGCACCGACTCTGCGCCGTCGGGGTa                                                    >  2:532761/1‑100 (MQ=255)
                                                aGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                                <  1:654183/100‑1 (MQ=255)
                                                  cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                              <  2:5191426/100‑1 (MQ=255)
                                                  cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                              <  2:925129/100‑1 (MQ=255)
                                                    cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCa                                            >  2:1589672/1‑100 (MQ=255)
                                                     gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                           <  1:1735433/100‑1 (MQ=255)
                                                     gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                           <  1:3157616/100‑1 (MQ=255)
                                                        gaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                        >  1:4778344/1‑100 (MQ=255)
                                                        gaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTCCTGAAGCGTGCGCCGACTCTGCACCGTCTGGGTATCCAGGCAtttc                                        >  1:3934118/1‑99 (MQ=255)
                                                            cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAcc                                    >  2:2532924/1‑100 (MQ=255)
                                                             tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                   >  2:234251/1‑100 (MQ=255)
                                                              gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGCACCgg                                  <  2:4840805/100‑1 (MQ=255)
                                                              gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                  >  1:4480378/1‑100 (MQ=255)
                                                              gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                  >  2:283777/1‑100 (MQ=255)
                                                                 gTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                               <  1:4766867/100‑1 (MQ=255)
                                                                    tGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                            >  2:5049194/1‑100 (MQ=255)
                                                                     gggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                           <  2:3146350/100‑1 (MQ=255)
                                                                      ggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                          <  2:2131069/100‑1 (MQ=255)
                                                                      ggATATCCTGGACGAAGTTATCCGCGACCACCAGGTACTGCTGAACCGTGCCCCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTTATc                          <  1:3058144/100‑1 (MQ=255)
                                                                        atatCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                        >  2:2180616/1‑100 (MQ=255)
                                                                         tatCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                       <  2:1875194/100‑1 (MQ=255)
                                                                           tCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                     >  2:3826792/1‑100 (MQ=255)
                                                                           tCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                     >  1:4990052/1‑100 (MQ=255)
                                                                           tCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                     <  2:2199554/100‑1 (MQ=255)
                                                                            ccTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                    >  2:2727306/1‑100 (MQ=255)
                                                                             cTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTa                   >  2:3675545/1‑100 (MQ=255)
                                                                             cTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTa                   >  1:1460926/1‑100 (MQ=255)
                                                                                gACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                >  2:1638172/1‑100 (MQ=255)
                                                                                gACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                <  2:2582347/100‑1 (MQ=255)
                                                                                 aCGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc               <  1:1864836/100‑1 (MQ=255)
                                                                                 aCGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc               <  1:4569458/100‑1 (MQ=255)
                                                                                  cGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt              <  2:4587851/100‑1 (MQ=255)
                                                                                  cGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt              >  1:4136646/1‑100 (MQ=255)
                                                                                  cGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt              <  1:1674163/100‑1 (MQ=255)
                                                                                        tATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg        >  1:1166761/1‑100 (MQ=255)
                                                                                         aTCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc       >  2:4215077/1‑100 (MQ=255)
                                                                                          tCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt      <  2:510197/100‑1 (MQ=255)
                                                                                           ccccgACCACCCGGTACTGCTGAACAGTGCACAGCCTCTCCACCATCTCAGTATCCAGGAATTTCAACCGGTACTACTCGAAGCTAACGCTATCCAGCTg     <  2:1362541/97‑1 (MQ=255)
                                                                                           ccGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg     <  1:2745817/100‑1 (MQ=255)
                                                                                            cgcgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGc    >  2:3907549/1‑100 (MQ=255)
                                                                                              cgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAc  >  1:3208563/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAC  >  NC_000913/4186506‑4186699

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