Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 367,140 C→T A207T (GCC→ACC)  mhpR ← transcriptional regulator for mhp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913367,1400CT94.7% 408.4 / NA 133A207T (GCC→ACC) mhpRtranscriptional regulator for mhp operon
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (42/83);  minor base G (7/0);  total (49/84)
Fisher's exact test for biased strand distribution p-value = 7.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGT  >  NC_000913/367050‑367228
                                                                                          |                                                                                        
gTCCGTAACCCTCTTTGCGCGCGCGCTCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTa                                                                                        <  1:2365671/93‑1 (MQ=255)
gTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTa                                                                                        <  2:2558140/93‑1 (MQ=255)
 tCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTaa                                                                                       <  2:3800399/93‑1 (MQ=255)
  ccGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAAc                                                                                      <  1:2037895/93‑1 (MQ=255)
   cGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACa                                                                                     <  1:3930197/93‑1 (MQ=255)
     tagcccGCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACAtt                                                                                   <  2:1607710/90‑1 (MQ=255)
     taccccTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACAtt                                                                                   <  2:584093/90‑1 (MQ=255)
     tAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACAtt                                                                                   >  1:615971/1‑93 (MQ=255)
     tAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACAtt                                                                                   <  2:3871853/93‑1 (MQ=255)
        cccTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCg                                                                                <  1:3311541/93‑1 (MQ=255)
         ccTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGa                                                                               >  1:693816/1‑93 (MQ=255)
         ccTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGa                                                                               >  2:1004135/1‑93 (MQ=255)
         ccTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGa                                                                               <  1:2343214/93‑1 (MQ=255)
         ccTCTTTGCGCGCGCGCGCCAGAATGGCTCCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGa                                                                               <  2:4043801/93‑1 (MQ=255)
          ctctTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGAt                                                                              <  2:871588/93‑1 (MQ=255)
           tctTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGGTAACATTTCGAtt                                                                             >  2:2443201/1‑93 (MQ=255)
           tctTTGCGCGCGCGCGCCAGAATGGCTGCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGAtt                                                                             <  1:1391771/93‑1 (MQ=255)
            ctTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTa                                                                            >  1:3727243/1‑93 (MQ=255)
            ctTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTa                                                                            <  2:1213702/93‑1 (MQ=255)
             tttGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTaa                                                                           <  2:1255596/93‑1 (MQ=255)
              ttGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAAt                                                                          <  2:1119665/93‑1 (MQ=255)
               tGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAAtt                                                                         <  1:3034086/93‑1 (MQ=255)
                 cgcgcgcgcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTcc                                                                       <  2:366295/93‑1 (MQ=255)
                  gcgcgcgcgcCAGAATGGCTTCCAGCTTTAACTGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCt                                                                      <  2:640992/93‑1 (MQ=255)
                  gcgcgcgcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCt                                                                      <  1:4141070/93‑1 (MQ=255)
                     cgcgcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGc                                                                   <  1:1340200/93‑1 (MQ=255)
                       cgcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCgg                                                                 <  1:1612863/93‑1 (MQ=255)
                       cgcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCgg                                                                 <  1:800015/93‑1 (MQ=255)
                        gcgcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGGTAACATTTCGATTAATTCCTTGCGGt                                                                >  1:3801819/1‑93 (MQ=255)
                          gcCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCt                                                              <  2:353500/93‑1 (MQ=255)
                             aGAATGTCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGt                                                           <  1:4107269/93‑1 (MQ=255)
                             aGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGt                                                           <  2:1726908/93‑1 (MQ=255)
                              gAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGtt                                                          >  2:1516658/1‑93 (MQ=255)
                              gAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGGTAACATTTCGATTAATTTCTTGCGGTgttgtt                                                          >  1:4079264/1‑93 (MQ=255)
                              gAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGGTAACATTTCGATTAATTCCTTGCGGTCTTGtt                                                          >  2:4436228/1‑93 (MQ=255)
                               aaTGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTc                                                         <  1:2014871/93‑1 (MQ=255)
                                aTGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTcc                                                        >  2:1691199/1‑93 (MQ=255)
                                aTGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTcc                                                        >  1:239326/1‑93 (MQ=255)
                                 tGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCg                                                       <  2:478462/93‑1 (MQ=255)
                                  ggCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTTGATTAATTCCTTGCGGTCTTGTTCCgg                                                      <  1:3562696/93‑1 (MQ=255)
                                  ggCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCgg                                                      <  1:572053/93‑1 (MQ=255)
                                     ttCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCa                                                   <  2:3633572/93‑1 (MQ=255)
                                     ttCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCa                                                   >  2:3114899/1‑93 (MQ=255)
                                      tCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCaa                                                  >  2:145417/1‑93 (MQ=255)
                                       ccAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCaaa                                                 >  2:3550024/1‑93 (MQ=255)
                                        cAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCaaaa                                                <  1:4258418/93‑1 (MQ=255)
                                        cAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCaaaa                                                >  2:2579008/1‑93 (MQ=255)
                                         aGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTGCTTGCGGTATTGTTCCGGGCaaacg                                               >  1:1475075/1‑91 (MQ=255)
                                         aGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAg                                               >  1:3927257/1‑93 (MQ=255)
                                         aGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAg                                               <  1:1140193/93‑1 (MQ=255)
                                            tttAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAggac                                            >  1:961475/1‑91 (MQ=255)
                                            tttAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgcc                                            >  1:2274963/1‑93 (MQ=255)
                                            tttAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgcc                                            <  1:2871971/93‑1 (MQ=255)
                                            tttAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgcc                                            >  2:679579/1‑93 (MQ=255)
                                            tttAACGGTGCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgcc                                            <  2:4503734/93‑1 (MQ=255)
                                              tAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAACAGgccag                                          >  2:2454775/1‑93 (MQ=255)
                                              tAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGggcag                                          >  2:855381/1‑93 (MQ=255)
                                               aaCGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGTGGTCTTGTTCCGGGCAAAAGgccagc                                         <  1:27504/93‑1 (MQ=255)
                                                aCGGTGCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGgccagcc                                        <  1:2617770/93‑1 (MQ=255)
                                                aCGGTCCCCGGGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATCCCTTGCGGTCTTGTTCCGGGCAAATGgccagcc                                        <  2:4549317/93‑1 (MQ=255)
                                                  ggTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGccagccag                                      <  1:817994/93‑1 (MQ=255)
                                                    ttCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGt                                    <  2:1546113/93‑1 (MQ=255)
                                                    ttCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGt                                    >  1:1768228/1‑93 (MQ=255)
                                                     tCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCACGTc                                   <  1:1292485/93‑1 (MQ=255)
                                                       ccGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAg                                 <  2:709705/93‑1 (MQ=255)
                                                         gTGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGc                               >  2:3244750/1‑93 (MQ=255)
                                                          tGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGcc                              <  2:136594/93‑1 (MQ=255)
                                                          tGCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGcc                              >  1:4644369/1‑93 (MQ=255)
                                                           gCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGccc                             <  1:349497/93‑1 (MQ=255)
                                                           gCCAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGccc                             >  2:2248041/1‑93 (MQ=255)
                                                             cAGTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGa                           <  1:1728207/93‑1 (MQ=255)
                                                               gTTGATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAgg                         <  1:4288189/93‑1 (MQ=255)
                                                                  gATAGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCgg                      >  2:4532879/1‑93 (MQ=255)
                                                                     aGTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGttt                   <  1:3295612/93‑1 (MQ=255)
                                                                      gTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGtttt                  >  1:437719/1‑93 (MQ=255)
                                                                      gTCATCACCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGtttt                  <  1:3320351/93‑1 (MQ=255)
                                                                       tcatcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTc                 >  1:1688864/1‑93 (MQ=255)
                                                                       tcatcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTc                 >  1:3644871/1‑93 (MQ=255)
                                                                       tcatcaCCGGGGCGGGAGGGTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTc                 >  2:4269332/1‑93 (MQ=255)
                                                                        catcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCa                <  1:4368710/93‑1 (MQ=255)
                                                                        catcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCa                <  1:4638823/93‑1 (MQ=255)
                                                                          tcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGa              <  1:1630459/93‑1 (MQ=255)
                                                                          tcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGa              >  1:3904349/1‑93 (MQ=255)
                                                                          tcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGa              >  1:4453648/1‑93 (MQ=255)
                                                                          tcaCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGGTTTCAGa              >  1:812147/1‑93 (MQ=255)
                                                                           caCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGaa             <  1:2290374/93‑1 (MQ=255)
                                                                           caCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGaa             >  1:1374927/1‑93 (MQ=255)
                                                                            aCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAg            >  1:3796179/1‑93 (MQ=255)
                                                                            aCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAg            <  1:2475302/93‑1 (MQ=255)
                                                                            aCCGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAg            >  2:2480107/1‑93 (MQ=255)
                                                                             ccGGGGCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGc           <  1:3157241/93‑1 (MQ=255)
                                                                               ggggCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCgg         >  2:4481237/1‑93 (MQ=255)
                                                                                  gCGGGAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCaa      >  2:3577399/1‑93 (MQ=255)
                                                                                   cGGTAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCaaa     <  2:2411339/93‑1 (MQ=255)
                                                                                    gggAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAAc    <  2:2151624/93‑1 (MQ=255)
                                                                                     ggAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACg   >  2:2576316/1‑93 (MQ=255)
                                                                                      gAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGt  <  1:3098025/93‑1 (MQ=255)
                                                                                      gAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGt  <  1:3473896/93‑1 (MQ=255)
                                                                                      gAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGt  <  1:2495040/93‑1 (MQ=255)
                                                                                      gAGGTTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGt  >  1:1943555/1‑93 (MQ=255)
                                                                                          |                                                                                        
GTCCGTAACCCTCTTTGCGCGCGCGCGCCAGAATGGCTTCCAGCTTTAACGGTTCCCGTGCCAGTTGATAGTCATCACCGGGGCGGGAGGCTAACATTTCGATTAATTCCTTGCGGTCTTGTTCCGGGCAAAAGGCCAGCCAGGTCAGGCCCGAGGCGGTTTTCAGAAGCGGCAAACGT  >  NC_000913/367050‑367228

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