Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,088,793 T→C T385T (ACA→ACG frvB ← PTS system, fructose‑like enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,7930TC85.7% 646.9 / 100.6 279T385T (ACA→ACGfrvBPTS system, fructose‑like enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (16/24);  new base C (122/117);  total (138/141)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 288 total aligned reads displayed.

CACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCG  >  NC_000913/4088703‑4088883
                                                                                          |                                                                                          
cACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgt                                                                                          <  1:1773322/93‑1 (MQ=255)
cACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgt                                                                                          <  2:6469811/93‑1 (MQ=255)
 aCGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgtg                                                                                         >  2:4653461/1‑93 (MQ=255)
  cGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTgg                                                                                        <  2:5255535/93‑1 (MQ=255)
   gCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGc                                                                                       >  1:544541/1‑93 (MQ=255)
      ggTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCAcc                                                                                    <  1:5106348/93‑1 (MQ=255)
      ggTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCAcc                                                                                    <  1:2158577/93‑1 (MQ=255)
        tAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGa                                                                                  >  1:6529648/1‑93 (MQ=255)
        tAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGa                                                                                  <  1:1043262/93‑1 (MQ=255)
           tGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACca                                                                               >  2:2540990/1‑93 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACcacat                                                                            >  2:3868639/1‑92 (MQ=255)
               gcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAAt                                                                           >  2:4397008/1‑93 (MQ=255)
                cccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATg                                                                          <  2:701829/93‑1 (MQ=255)
                  cgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGc                                                                        <  1:6901543/93‑1 (MQ=255)
                   gccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCa                                                                       <  1:7384054/93‑1 (MQ=255)
                    ccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCac                                                                      <  2:6952291/93‑1 (MQ=255)
                    ccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCac                                                                      >  1:4578574/1‑93 (MQ=255)
                    ccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCac                                                                      <  2:4263937/93‑1 (MQ=255)
                     cAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCact                                                                     >  1:6086656/1‑93 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                   >  1:682806/1‑93 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                   >  2:225963/1‑93 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                   <  1:5786436/93‑1 (MQ=255)
                        cGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCactact                                                                  >  2:1561081/1‑93 (MQ=255)
                        cGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactact                                                                  <  1:975948/93‑1 (MQ=255)
                        cGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactact                                                                  <  2:7131624/93‑1 (MQ=255)
                         gTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTg                                                                 >  2:4046297/1‑93 (MQ=255)
                         gTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTg                                                                 <  1:2517193/93‑1 (MQ=255)
                         gTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTg                                                                 >  2:550734/1‑93 (MQ=255)
                            ttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCa                                                              >  2:5933899/1‑93 (MQ=255)
                              ggCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGc                                                            <  1:6023289/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                         <  1:4018450/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                         <  2:4756687/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                         >  1:2986695/1‑93 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                        <  1:114543/93‑1 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                        >  2:4982161/1‑93 (MQ=255)
                                  gTAATCATCGGCACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtt                                                       >  1:620960/1‑93 (MQ=255)
                                    aaTCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCttt                                                      >  1:4213668/1‑93 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                     <  1:6680850/93‑1 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                     >  2:2675271/1‑93 (MQ=255)
                                        atcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTcc                                                  <  1:641353/93‑1 (MQ=255)
                                        atcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTcc                                                  >  1:2304310/1‑93 (MQ=255)
                                        atcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTcc                                                  >  1:1212709/1‑93 (MQ=255)
                                        atcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTcc                                                  >  2:444884/1‑93 (MQ=255)
                                          cGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTg                                                <  1:5563676/93‑1 (MQ=255)
                                           ggCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGt                                               <  2:4329549/93‑1 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                              >  1:5871025/1‑93 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                              <  2:2190143/93‑1 (MQ=255)
                                             cAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTcc                                             <  1:4874999/93‑2 (MQ=255)
                                              aaCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCgg                                            <  2:1333690/93‑1 (MQ=255)
                                                cGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGt                                          >  1:1361820/1‑93 (MQ=255)
                                                 ggAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTa                                         >  2:3611239/1‑93 (MQ=255)
                                                  gAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTaa                                        >  1:3854643/1‑93 (MQ=255)
                                                  gAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTaa                                        >  2:6218919/1‑93 (MQ=255)
                                                   aGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTaaa                                       <  1:3106024/93‑1 (MQ=255)
                                                    gCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                      >  2:1131526/1‑93 (MQ=255)
                                                    gCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                      <  2:3433626/93‑1 (MQ=255)
                                                    gCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                      <  2:6278616/93‑1 (MQ=255)
                                                    gCTGCCAGAGCGCAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                      <  1:5135408/93‑1 (MQ=255)
                                                     cTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACa                                     <  1:5047033/93‑1 (MQ=255)
                                                      tGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAg                                    >  1:2061698/1‑93 (MQ=255)
                                                       gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGt                                   <  2:5705524/93‑1 (MQ=255)
                                                       gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGt                                   <  2:3839352/93‑1 (MQ=255)
                                                        ccAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGtt                                  <  1:4827129/93‑1 (MQ=255)
                                                         cAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGttt                                 >  1:4136966/1‑93 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                                >  2:7372747/1‑93 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                                <  2:1373590/93‑1 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                                >  1:4297016/1‑93 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                               <  2:614590/93‑1 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                               <  2:751314/93‑1 (MQ=255)
                                                            agCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTggg                              <  1:3366776/93‑1 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgcgc                           <  2:2576525/93‑1 (MQ=255)
                                                                tAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCa                          <  2:2057157/93‑1 (MQ=255)
                                                                tAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCa                          >  2:3736709/1‑93 (MQ=255)
                                                                 aaGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCCCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCaa                         <  2:1820048/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                        <  1:4434518/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                        >  1:3903434/1‑93 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                        <  1:5309414/93‑1 (MQ=255)
                                                                     aaTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGaa                     <  2:4904354/93‑1 (MQ=255)
                                                                     aaTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGaa                     <  2:5459336/93‑1 (MQ=255)
                                                                       tAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACg                   >  1:2186720/1‑93 (MQ=255)
                                                                       tAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACg                   <  1:1957143/93‑1 (MQ=255)
                                                                       tAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACg                   >  2:659251/1‑93 (MQ=255)
                                                                         gCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTc                 <  2:87779/93‑1 (MQ=255)
                                                                         gCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTc                 <  2:3309557/93‑1 (MQ=255)
                                                                          ccGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCg                >  1:7046594/1‑93 (MQ=255)
                                                                            ggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCa              >  2:7248775/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa             <  2:4700202/93‑1 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa             <  1:2374795/93‑1 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa             >  1:1845471/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa             >  2:5859242/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa             >  2:2805408/1‑93 (MQ=255)
                                                                              ttcggttGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAg            >  1:3436914/1‑93 (MQ=255)
                                                                              ttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAg            <  1:7133371/93‑1 (MQ=255)
                                                                                cggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAc          <  1:6756213/93‑1 (MQ=255)
                                                                                 ggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAcc         <  2:5083435/93‑1 (MQ=255)
                                                                                  gttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAccc        <  2:5378717/93‑1 (MQ=255)
                                                                                   ttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACccg       >  2:3018647/1‑93 (MQ=255)
                                                                                        aCTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCcgcg  >  2:1407443/1‑93 (MQ=255)
                                                                                        aCTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCcgcg  <  2:1207827/93‑1 (MQ=255)
                                                                                        aCTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCCTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCcgcg  >  2:3198424/1‑93 (MQ=255)
                                                                                          |                                                                                          
CACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCG  >  NC_000913/4088703‑4088883

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