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.7% 529.0 / ‑4.4 154E418D (GAA→GAC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (70/81);  total (71/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

ACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGC  >  NC_000913/4186507‑4186697
                                                                                                |                                                                                              
aCTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAccac                                                                                             >  1:489121/1‑100 (MQ=255)
  tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCAccc                                                                                           <  1:3453682/100‑1 (MQ=255)
  tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCAccc                                                                                           >  1:1312054/1‑100 (MQ=255)
   gCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCg                                                                                          >  2:560894/1‑100 (MQ=255)
    cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCgg                                                                                         <  2:215217/100‑1 (MQ=255)
     gTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGt                                                                                        <  2:2261774/100‑1 (MQ=255)
     gTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGt                                                                                        >  1:2929506/1‑100 (MQ=255)
      tGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTa                                                                                       >  2:503785/1‑100 (MQ=255)
        gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGCAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTAct                                                                                     >  2:2923320/1‑100 (MQ=255)
         tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActg                                                                                    >  2:3562626/1‑100 (MQ=255)
          cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                   >  1:2418634/1‑100 (MQ=255)
          cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                   >  2:2284373/1‑100 (MQ=255)
          cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                   >  2:2850993/1‑100 (MQ=255)
          cTTGCTACCAACATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgc                                                                                   <  1:1012478/100‑1 (MQ=255)
           ttGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTActgct                                                                                  <  2:749462/100‑1 (MQ=255)
              cTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGaa                                                                               <  2:2685872/100‑1 (MQ=255)
               tACCACCATTAAAGCTGCGCAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTCGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAAc                                                                              <  1:3626697/100‑1 (MQ=255)
                 ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCg                                                                            <  2:1144626/100‑1 (MQ=255)
                 ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCg                                                                            <  1:775917/100‑1 (MQ=255)
                  caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                                           >  2:2503993/1‑100 (MQ=255)
                     caTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCa                                                                        >  1:1588665/1‑100 (MQ=255)
                       ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCAcc                                                                      <  2:441544/100‑1 (MQ=255)
                         aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGa                                                                    >  1:1515098/1‑100 (MQ=255)
                            gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTGGTCTGGACGAAGTGATCAGCGATCAACCGGAACTGCTGAACCGTGCACCGCctc                                                                 <  2:3849520/100‑1 (MQ=255)
                            gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActc                                                                 >  1:929310/1‑100 (MQ=255)
                            gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGActc                                                                 <  2:728646/100‑1 (MQ=255)
                                cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                             <  1:2263487/100‑1 (MQ=255)
                                cGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                             >  1:1342011/1‑100 (MQ=255)
                                 gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                            <  1:3341906/100‑1 (MQ=255)
                                 gaagaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTCCAc                                                            >  2:611052/1‑100 (MQ=255)
                                     aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                        <  1:1455632/100‑1 (MQ=255)
                                      aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATGCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                       >  2:1796278/1‑100 (MQ=255)
                                        aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCTGCGTCCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                     <  1:3581513/100‑1 (MQ=255)
                                         tGGTTGAGCGTGAAGGATCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGAGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                    <  1:1906571/100‑1 (MQ=255)
                                         tGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                    <  1:2141430/100‑1 (MQ=255)
                                         tGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                    <  2:780079/100‑1 (MQ=255)
                                           gTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTa                                                  >  1:2499251/1‑100 (MQ=255)
                                            ttGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTAt                                                 <  2:1513166/100‑1 (MQ=255)
                                              gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                               >  1:664324/1‑100 (MQ=255)
                                               aGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                                              <  1:823743/100‑1 (MQ=255)
                                                 cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                            >  2:1708439/1‑100 (MQ=255)
                                                 cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAgg                                            <  1:1333121/100‑1 (MQ=255)
                                                 cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACGCTGCACCGTCTGGGTATCCAgg                                            >  1:2025312/1‑100 (MQ=255)
                                                  gcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                           >  1:313836/1‑100 (MQ=255)
                                                  gcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGc                                           >  2:1983333/1‑100 (MQ=255)
                                                   cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCa                                          <  2:1393720/100‑1 (MQ=255)
                                                   cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCa                                          >  2:902543/1‑100 (MQ=255)
                                                   cgAAGAAGCTGTCGTTTGGGATAGCCTTGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGGGCACCGACTTTGCACCGTCTGGGTATCCAGGCa                                          <  2:1830000/100‑1 (MQ=255)
                                                    gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                         >  1:2163653/1‑100 (MQ=255)
                                                      agaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                       <  1:1846596/100‑1 (MQ=255)
                                                       gaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                      >  1:2062838/1‑100 (MQ=255)
                                                         agCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                    >  1:3057666/1‑100 (MQ=255)
                                                          gCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                   >  1:3444157/1‑100 (MQ=255)
                                                          gCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                   <  1:1530356/100‑1 (MQ=255)
                                                           cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAcc                                  <  1:670063/100‑1 (MQ=255)
                                                            tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                 >  1:3323444/1‑100 (MQ=255)
                                                            tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                 >  1:4067006/1‑100 (MQ=255)
                                                             gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                <  1:3370700/100‑1 (MQ=255)
                                                             gTCGTTCGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                <  1:3370715/100‑1 (MQ=255)
                                                              tCGTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                               <  1:2999591/100‑1 (MQ=255)
                                                               cGTTTGGGGTATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCTGGAATTTGATGCGGTa                              >  2:2507157/1‑100 (MQ=255)
                                                                gTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                             <  1:1093548/100‑1 (MQ=255)
                                                                gTTTGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                             >  2:3743824/1‑100 (MQ=255)
                                                                 tttGGGATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACt                            <  1:3672558/100‑1 (MQ=255)
                                                                    gggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                         >  2:2080701/1‑100 (MQ=255)
                                                                    gggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                         <  1:2319791/100‑1 (MQ=255)
                                                                    gggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                         <  1:2319829/100‑1 (MQ=255)
                                                                    gggATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTG‑ACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                        <  2:1824610/100‑1 (MQ=255)
                                                                     ctgtatCCTGGACGAAGCTATCCGCGTCCACCCGGTACTGCTGAACCGCGCACCGACTCTTCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATc                        <  2:3622346/97‑1 (MQ=255)
                                                                      gATATCCTGGACGAAGTTATCCGCGACCCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCCGGGTATCCAGGCATTTGAACCGGTACTGATCg                       >  1:1226220/1‑100 (MQ=255)
                                                                      gATATCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTATGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCg                       >  2:1113710/1‑100 (MQ=255)
                                                                        tatCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                     >  2:2397297/1‑100 (MQ=255)
                                                                         atCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                    <  1:2654987/100‑1 (MQ=255)
                                                                          tCCTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAgg                   <  1:1271276/100‑1 (MQ=255)
                                                                           ccTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTCCTGATCGAAGGt                  >  1:3229259/1‑100 (MQ=255)
                                                                           ccTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                  >  1:2594826/1‑100 (MQ=255)
                                                                           ccTGGACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                  >  1:2624756/1‑100 (MQ=255)
                                                                               gACGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg              <  1:2887010/100‑1 (MQ=255)
                                                                                aCGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc             <  2:3480179/100‑1 (MQ=255)
                                                                                aCGAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGc             <  1:2808749/100‑1 (MQ=255)
                                                                                 cagaGTTAGCCGCGTCCACCCGGTACTGCTGAAGCGTTCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCt            <  1:1018441/97‑1 (MQ=255)
                                                                                  gAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa           <  1:2261028/100‑1 (MQ=255)
                                                                                  gAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa           <  2:2381504/100‑1 (MQ=255)
                                                                                  gAAGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTa           >  2:2614056/1‑100 (MQ=255)
                                                                                   aaGTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTAt          >  1:603012/1‑100 (MQ=255)
                                                                                     gTTATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATcc        <  1:1161642/100‑1 (MQ=255)
                                                                                      ttATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCa       >  2:2704161/1‑100 (MQ=255)
                                                                                      ttATCCACGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCa       <  1:3108969/100‑1 (MQ=255)
                                                                                       tATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg      >  1:2751748/1‑100 (MQ=255)
                                                                                       tATCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAg      <  2:3830915/100‑1 (MQ=255)
                                                                                        tctcgcgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc     <  1:1614557/97‑1 (MQ=255)
                                                                                        aTCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc     <  2:1362281/100‑1 (MQ=255)
                                                                                        aTCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc     <  2:41802/100‑1 (MQ=255)
                                                                                         tCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCTAGCt    >  2:1543732/1‑100 (MQ=255)
                                                                                         tCCGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt    <  1:105292/100‑1 (MQ=255)
                                                                                          ccGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg   <  1:2099576/100‑1 (MQ=255)
                                                                                          ccGCGACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTg   <  2:2885601/100‑1 (MQ=255)
                                                                                           cgcgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTCTCCAGCTGc  >  1:610970/1‑100 (MQ=255)
                                                                                           cgcgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGc  >  1:898063/1‑100 (MQ=255)
                                                                                           cgcgACCACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATGTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGc  <  2:1192214/100‑1 (MQ=255)
                                                                                                |                                                                                              
ACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGC  >  NC_000913/4186507‑4186697

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