Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,213,859 C→T N271N (AAC→AAT rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,213,8590CT100.0% 1799.9 / NA 415N271N (AAC→AATrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (220/195);  total (220/195)

Only 100 of 415 total aligned reads displayed.

GGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAACAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTCATTA  >  NC_000913/3213761‑3213956
                                                                                                  |                                                                                                 
ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAg                                                                                                 >  1:5840441/1‑101 (MQ=255)
  tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCa                                                                                               <  1:5840027/101‑1 (MQ=255)
   cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCAt                                                                                              >  2:5840444/1‑101 (MQ=255)
    gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATg                                                                                             <  1:5839739/101‑1 (MQ=255)
     cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATgc                                                                                            <  1:5839798/101‑1 (MQ=255)
     cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATgc                                                                                            <  1:5839844/101‑1 (MQ=255)
     cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATgc                                                                                            >  1:5840450/1‑101 (MQ=255)
      aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATgcg                                                                                           <  2:5839807/101‑1 (MQ=255)
        tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGcgcg                                                                                         >  2:5840453/1‑101 (MQ=255)
         cACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGt                                                                                        >  1:5840455/1‑101 (MQ=255)
             cTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCatg                                                                                    >  2:5840459/1‑101 (MQ=255)
             cTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCatg                                                                                    <  1:5839892/101‑1 (MQ=255)
              tACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCatga                                                                                   <  2:5840112/101‑1 (MQ=255)
                 cGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATgg                                                                                >  1:5840466/1‑101 (MQ=255)
                  gCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGa                                                                               >  2:5840468/1‑101 (MQ=255)
                    tCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGAcc                                                                             >  2:5840473/1‑101 (MQ=255)
                    tCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGAcc                                                                             >  1:5840474/1‑101 (MQ=255)
                       ggAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACcgcg                                                                          >  2:5840484/1‑101 (MQ=255)
                       ggAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACcgcg                                                                          >  1:5840480/1‑101 (MQ=255)
                         aaGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGtt                                                                        >  1:5840490/1‑101 (MQ=255)
                         aaGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGtt                                                                        <  2:5839878/101‑1 (MQ=255)
                          agagATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTc                                                                       <  1:5839997/101‑1 (MQ=255)
                           gagaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCg                                                                      >  2:5840495/1‑101 (MQ=255)
                            agaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGt                                                                     <  1:5840115/101‑1 (MQ=255)
                            agaTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGt                                                                     >  1:5840496/1‑101 (MQ=255)
                              aTCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACcgc                                                                           <  1:5839918/93‑1 (MQ=255)
                               tCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACg                                                                  >  2:5840502/1‑101 (MQ=255)
                                 cTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCa                                                                >  1:5840507/1‑101 (MQ=255)
                                  tGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCaa                                                               >  2:5840509/1‑101 (MQ=255)
                                   gAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAg                                                              <  1:5840101/101‑1 (MQ=255)
                                    aaaCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTGCGCAAGa                                                             >  2:5840516/1‑101 (MQ=255)
                                    aaaCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCaa                                                               <  2:5840031/99‑1 (MQ=255)
                                    aaaCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGa                                                             <  2:5839722/101‑1 (MQ=255)
                                      aCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGAcc                                                                             <  2:5839965/83‑1 (MQ=255)
                                      aCTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAAc                                                           <  2:5840037/101‑1 (MQ=255)
                                       cTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACg                                                          >  2:5840521/1‑101 (MQ=255)
                                       cTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACg                                                          >  1:5840522/1‑101 (MQ=255)
                                         gTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTc                                                        <  2:5840181/101‑1 (MQ=255)
                                           cTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTg                                                      >  2:5840530/1‑101 (MQ=255)
                                            tGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGa                                                     >  1:5840534/1‑101 (MQ=255)
                                            tGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGa                                                     >  2:5840536/1‑101 (MQ=255)
                                              aaGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATc                                                   <  2:5839935/101‑1 (MQ=255)
                                              aaGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATc                                                   >  1:5840543/1‑101 (MQ=255)
                                              aaGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATc                                                   >  2:5840544/1‑101 (MQ=255)
                                                 tATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATg                                                <  1:5839882/101‑1 (MQ=255)
                                                 tATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATg                                                <  1:5840059/101‑1 (MQ=255)
                                                 tATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATg                                                <  1:5839533/101‑1 (MQ=255)
                                                  aTTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGa                                               >  1:5840552/1‑101 (MQ=255)
                                                   ttCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGaa                                              <  1:5840164/101‑1 (MQ=255)
                                                   ttCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGaa                                              <  2:5839912/101‑1 (MQ=255)
                                                   ttCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGaa                                              <  1:5840207/101‑1 (MQ=255)
                                                   ttCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGaa                                              >  2:5840555/1‑101 (MQ=255)
                                                    tCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAg                                             >  1:5840558/1‑101 (MQ=255)
                                                     cAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGc                                            >  1:5840562/1‑101 (MQ=255)
                                                      aaaCAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCAt                                                 >  2:5840566/1‑95 (MQ=255)
                                                        aCAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGctct                                         >  1:5840571/1‑101 (MQ=255)
                                                        aCAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGctct                                         <  2:5839982/101‑1 (MQ=255)
                                                         cAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTg                                        >  1:5840575/1‑101 (MQ=255)
                                                           gTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCg                                      >  1:5840581/1‑101 (MQ=255)
                                                           gTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCg                                      >  2:5840580/1‑101 (MQ=255)
                                                             tCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCAt                                                                                              >  1:5840587/1‑43 (MQ=255)
                                                             tCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGtt                                    <  1:5840083/101‑1 (MQ=255)
                                                               cGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGa                                  >  2:5840590/1‑101 (MQ=255)
                                                                  cTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCa                               >  1:5840597/1‑101 (MQ=255)
                                                                  cTGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCa                               >  2:5840598/1‑101 (MQ=255)
                                                                   tGGTGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAg                              <  1:5839975/101‑1 (MQ=255)
                                                                      tGCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGc                           <  2:5840205/101‑1 (MQ=255)
                                                                       gCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCa                          <  1:5839665/101‑1 (MQ=255)
                                                                       gCCGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCa                          >  1:5840606/1‑101 (MQ=255)
                                                                        ccGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaa                         <  1:5839953/101‑1 (MQ=255)
                                                                         cGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCGTGAAGCTCTGCGTTGAGCAGTGCaaa                        <  2:5840048/101‑1 (MQ=255)
                                                                         cGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaa                        <  1:5840069/101‑1 (MQ=255)
                                                                         cGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaa                        <  2:5840006/101‑1 (MQ=255)
                                                                         cGAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaa                        >  1:5840613/1‑101 (MQ=255)
                                                                          gAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaaa                       <  1:5839938/101‑1 (MQ=255)
                                                                          gAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaaa                       >  2:5840614/1‑101 (MQ=255)
                                                                          gAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCaaaa                       <  1:5840257/101‑1 (MQ=255)
                                                                          gAAGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCCCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCa                          >  1:5840619/1‑98 (MQ=255)
                                                                            aGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATg                     <  1:5840341/101‑1 (MQ=255)
                                                                            aGCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATg                     <  1:5840133/101‑1 (MQ=255)
                                                                             gCAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGc                    >  1:5840624/1‑101 (MQ=255)
                                                                              cAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGcc                   <  2:5840231/101‑1 (MQ=255)
                                                                              cAGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGcc                   >  2:5840627/1‑101 (MQ=255)
                                                                               aGTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCg                  <  2:5840145/101‑1 (MQ=255)
                                                                                gTTTGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCga                 <  2:5840096/101‑1 (MQ=255)
                                                                                  ttGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCgaag               >  1:5840632/1‑101 (MQ=255)
                                                                                  ttGACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCgaag               <  1:5840014/101‑1 (MQ=255)
                                                                                    gACTACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCgaagaa             >  1:5840637/1‑101 (MQ=255)
                                                                                       tACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGaaaaa          <  2:5840071/101‑1 (MQ=255)
                                                                                       tACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGaaaaa          >  2:5840643/1‑101 (MQ=255)
                                                                                       tACCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGaaaaa          >  1:5840642/1‑101 (MQ=255)
                                                                                        aCCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAAc         <  1:5840046/101‑1 (MQ=255)
                                                                                        aCCTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAAc         <  1:5839987/101‑1 (MQ=255)
                                                                                          cTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGaaaaa          <  1:5840219/98‑1 (MQ=255)
                                                                                          cTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACtt       <  2:5840606/101‑1 (MQ=255)
                                                                                          cTGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACtt       <  2:5840313/101‑1 (MQ=255)
                                                                                           tGGTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTc      >  2:5840649/1‑101 (MQ=255)
                                                                                             gTCAATAGCATGCGCGTCATGATGGATCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTCAt    <  1:5839875/101‑1 (MQ=255)
                                                                                             gTCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTCAt    >  2:5840654/1‑101 (MQ=255)
                                                                                              tCAATAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTCAtt   >  1:5840655/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTACCTGGTCAACAGCATGCGCGTCATGATGGACCGCGTTCGTACGCAAGAACGTCTGATCATGAAGCTCTGCGTTGAGCAGTGCAAAATGCCGAAGAAAAACTTCATTA  >  NC_000913/3213761‑3213956

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