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,7930TC99.6% 793.7 / ‑4.6 230T385T (ACA→ACGfrvBPTS system, fructose‑like enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (105/124);  total (105/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

Only 100 of 241 total aligned reads displayed.

CACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGC  >  NC_000913/4088703‑4088882
                                                                                          |                                                                                         
   gCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGc                                                                                      <  2:6243751/93‑1 (MQ=255)
   cccGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGc                                                                                      <  1:6135343/92‑1 (MQ=255)
    ccGTTAATCCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCa                                                                                     <  1:2583075/93‑1 (MQ=255)
     cGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCAc                                                                                    <  1:5234124/93‑1 (MQ=255)
     cGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCAc                                                                                    >  1:2903303/1‑93 (MQ=255)
      ggTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCAcc                                                                                   >  1:4908078/1‑93 (MQ=255)
        tAAGGCGGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGTCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGa                                                                                 <  2:1525133/93‑1 (MQ=255)
           tGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACca                                                                              >  2:5064559/1‑93 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACaa                                                                           <  2:3141026/93‑1 (MQ=255)
               gcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAAt                                                                          >  2:4010183/1‑93 (MQ=255)
               gcccgccAGCATATTTGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAAt                                                                          >  2:5429903/1‑93 (MQ=255)
               ccccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAAt                                                                          <  1:1926945/92‑1 (MQ=255)
                cccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATg                                                                         <  1:914755/93‑1 (MQ=255)
                 ccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAACGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGGACATGCCACAATgg                                                                        >  1:6435603/1‑93 (MQ=255)
                  cgccAGCGTATTGGGTGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGc                                                                       <  1:1744992/93‑1 (MQ=255)
                  cgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGc                                                                       <  1:5970106/93‑1 (MQ=255)
                  cgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGc                                                                       <  1:3269825/93‑1 (MQ=255)
                     cAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACGATGGCact                                                                    >  2:1220307/1‑93 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                  >  2:3977355/1‑93 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                  <  2:372847/93‑1 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                  <  2:4944016/93‑1 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactac                                                                  <  1:2040785/93‑1 (MQ=255)
                       gCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCaatac                                                                  >  1:3583322/1‑93 (MQ=255)
                        cGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactact                                                                 <  2:3344324/93‑1 (MQ=255)
                           aTTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGc                                                              <  2:39253/93‑1 (MQ=255)
                           aTTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGc                                                              <  2:5547022/93‑1 (MQ=255)
                                cGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCct                                                         <  2:245031/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                        <  2:2792454/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                        <  2:873386/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                        <  1:5654352/93‑1 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                       <  1:3424134/93‑1 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                       <  2:2653298/93‑1 (MQ=255)
                                   tAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtt                                                      >  1:2615003/1‑93 (MQ=255)
                                   tAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtt                                                      <  1:2209641/93‑1 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                    <  1:1216189/93‑1 (MQ=255)
                                      tcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCttttt                                                   >  2:489716/1‑93 (MQ=255)
                                      gcatcGGGAACGGAGCGGCCCGAGCGTAAGGAATAGCCGGTTAGGTTGCGAACGTGGCACCGACCACAATGGCACTACTGGCAGCCTCttttt                                                   <  1:2604902/92‑1 (MQ=255)
                                        atcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTcc                                                 >  2:2867789/1‑93 (MQ=255)
                                         tcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCt                                                <  1:5876443/93‑1 (MQ=255)
                                         tcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCt                                                <  2:1510476/93‑1 (MQ=255)
                                           ggCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGt                                              >  2:1795270/1‑93 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                             <  2:262171/93‑1 (MQ=255)
                                              aaCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCgg                                           <  1:3770984/93‑1 (MQ=255)
                                              aaCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCgg                                           <  1:348405/93‑1 (MQ=255)
                                              aaCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCgg                                           <  1:5621187/93‑1 (MQ=255)
                                              aaCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCgg                                           >  1:1353790/1‑93 (MQ=255)
                                                cGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGt                                         <  1:3308684/93‑1 (MQ=255)
                                                cGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGt                                         >  1:2771468/1‑93 (MQ=255)
                                                cGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGt                                         <  1:1035368/93‑1 (MQ=255)
                                                 ggAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGGACCGACCACAATGGCAATAATGGCAGCCTCCTTTTCCTGTCGGGTa                                        >  2:4440703/1‑93 (MQ=255)
                                                 ggAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTTCTGTCGGGTa                                        >  2:1167848/1‑93 (MQ=255)
                                                  gAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTaa                                       <  2:2512641/93‑1 (MQ=255)
                                                   aGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTaaa                                      >  1:1613661/1‑93 (MQ=255)
                                                     cTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACa                                    <  2:6413778/93‑1 (MQ=255)
                                                      tGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAg                                   >  1:3610789/1‑93 (MQ=255)
                                                       gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGt                                  <  1:6472734/93‑1 (MQ=255)
                                                       gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGt                                  <  1:4289268/93‑1 (MQ=255)
                                                        ccAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGtt                                 <  2:2939603/93‑1 (MQ=255)
                                                         cAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGttt                                <  2:6106395/93‑1 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                              <  1:1998547/93‑1 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                              >  1:1251432/1‑93 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                              <  1:4635798/93‑1 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                              >  2:2040217/1‑93 (MQ=255)
                                                             gCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgc                            >  1:1153862/1‑93 (MQ=255)
                                                             gCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgc                            <  1:2318951/93‑1 (MQ=255)
                                                             gCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgc                            >  1:3564113/1‑93 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCTGGTAAACAGTTTGGgcgc                          >  1:4674327/1‑93 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgcgc                          >  1:3817175/1‑93 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgcgc                          >  2:2723280/1‑93 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgcgc                          <  2:3837404/93‑1 (MQ=255)
                                                               gTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgcgc                          >  1:42739/1‑93 (MQ=255)
                                                                 aaGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCaa                        >  2:1595975/1‑93 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGTGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       <  2:3894758/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       <  2:6062641/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       >  1:5408605/1‑93 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGCAAACAGTTTGGGCGCAAt                       >  1:1114370/1‑93 (MQ=255)
                                                                   ggAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATg                      >  1:3534618/1‑93 (MQ=255)
                                                                   ggAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATg                      <  1:637145/93‑1 (MQ=255)
                                                                   ggAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATg                      >  2:3840555/1‑93 (MQ=255)
                                                                    gAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGa                     >  2:3396473/1‑93 (MQ=255)
                                                                    gAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGa                     >  1:18404/1‑93 (MQ=255)
                                                                     aaTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGaa                    <  1:740517/93‑1 (MQ=255)
                                                                      aTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAc                   >  1:157681/1‑93 (MQ=255)
                                                                       tAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACg                  <  2:5898381/93‑1 (MQ=255)
                                                                           cggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGc              >  1:4061615/1‑93 (MQ=255)
                                                                            ggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCa             <  1:5106887/93‑1 (MQ=255)
                                                                            ggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCa             >  1:6092105/1‑93 (MQ=255)
                                                                            ggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCa             >  2:2228753/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa            >  2:4199896/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa            >  2:1657303/1‑93 (MQ=255)
                                                                              ttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAg           >  1:3167640/1‑93 (MQ=255)
                                                                               tcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGa          >  1:5027177/1‑93 (MQ=255)
                                                                                cggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAc         >  1:3926651/1‑93 (MQ=255)
                                                                                 ggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAcc        <  2:5619366/93‑1 (MQ=255)
                                                                                 ggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAcc        >  2:4766210/1‑93 (MQ=255)
                                                                                   ttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACccg      >  1:5747314/1‑93 (MQ=255)
                                                                                     gCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACccgcc    <  2:2746951/93‑1 (MQ=255)
                                                                                       gACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  <  2:4889116/93‑1 (MQ=255)
                                                                                       gACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  >  2:2271494/1‑93 (MQ=255)
                                                                                       gACCGTGCCACCGACCACAATGGCACTACTGTCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  <  2:743577/93‑1 (MQ=255)
                                                                                          |                                                                                         
CACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGC  >  NC_000913/4088703‑4088882

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