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.3% 909.5 / ‑3.7 270T385T (ACA→ACGfrvBPTS system, fructose‑like enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (114/154);  total (115/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 288 total aligned reads displayed.

CACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGC  >  NC_000913/4088703‑4088882
                                                                                          |                                                                                         
cACGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgt                                                                                         <  2:4039858/93‑1 (MQ=255)
 ccGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgtg                                                                                        <  2:211840/92‑1 (MQ=255)
 aCGCCGGTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACcgtg                                                                                        <  2:3063193/93‑1 (MQ=255)
       gTAATGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCg                                                                                  >  2:4120265/1‑93 (MQ=255)
         aaTGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGAc                                                                                <  2:5463031/93‑1 (MQ=255)
         aaTGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGAc                                                                                >  1:5926662/1‑93 (MQ=255)
          aTGCCGCCCGCCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCCTGGCACCGAcc                                                                               >  1:2589345/1‑93 (MQ=255)
             ccgcccgcCAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACcaca                                                                            >  2:8253289/1‑93 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACaa                                                                           >  2:5573953/1‑93 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACaa                                                                           >  2:3564798/1‑93 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACaa                                                                           <  1:4197978/93‑1 (MQ=255)
              cgcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACaa                                                                           <  1:666063/93‑1 (MQ=255)
               gcccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAAt                                                                          <  1:3800120/93‑1 (MQ=255)
                cccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATg                                                                         <  1:654552/93‑1 (MQ=255)
                 ccgccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATgg                                                                        <  2:5882207/93‑1 (MQ=255)
                   gccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCAACCACAATGGCa                                                                      <  2:4421921/93‑1 (MQ=255)
                    ccAGCGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCac                                                                     <  2:1792987/93‑1 (MQ=255)
                        cGTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCactact                                                                 >  1:1355585/1‑93 (MQ=255)
                         gtgttGGCGGTGATCATCGGCAGCGGAGCTGCCAGAGCGTACGGAATAGCCGGTTCGGTTGCGGCCGTGGCACCGACCACAATGGCACTACTg                                                                <  1:6684830/90‑1 (MQ=255)
                         gTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTg                                                                <  1:4558353/93‑1 (MQ=255)
                         gTATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTg                                                                <  1:5091967/93‑1 (MQ=255)
                          tATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTgg                                                               <  2:7138528/93‑1 (MQ=255)
                          tATTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTgg                                                               >  2:6952745/1‑93 (MQ=255)
                           aTTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTgtc                                                              >  1:6837943/1‑91 (MQ=255)
                           aTTGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGc                                                              >  1:237130/1‑93 (MQ=255)
                            ttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCa                                                             >  1:6757412/1‑93 (MQ=255)
                            ttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCa                                                             <  1:1097599/93‑1 (MQ=255)
                            ttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCa                                                             <  2:419865/93‑1 (MQ=255)
                            ttGGCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCa                                                             <  1:2218675/93‑1 (MQ=255)
                              ggCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGc                                                           <  2:3138020/93‑1 (MQ=255)
                               gCGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGcc                                                          >  2:1276705/1‑93 (MQ=255)
                                cGGTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCct                                                         <  1:7513732/93‑1 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                        >  2:6567678/1‑93 (MQ=255)
                                 ggTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctc                                                        <  2:3452234/93‑1 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                       <  1:6471437/93‑1 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                       >  1:1941208/1‑93 (MQ=255)
                                  gTAATCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCctct                                                       <  2:1848331/93‑1 (MQ=255)
                                    aaTCATCGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCttt                                                     <  1:8290289/93‑1 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                    <  1:8202706/93‑1 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                    <  1:101793/93‑1 (MQ=255)
                                     atcatcGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCtttt                                                    <  1:1608446/93‑1 (MQ=255)
                                      tcatcaGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCttttt                                                   >  1:7366043/1‑93 (MQ=255)
                                          cGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTg                                               <  2:4315186/93‑1 (MQ=255)
                                          cGGCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTg                                               <  2:5423798/93‑1 (MQ=255)
                                           ggCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGt                                              <  1:4120019/93‑1 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                             <  2:6118964/93‑1 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                             <  2:2016244/93‑1 (MQ=255)
                                            gCAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTc                                             <  1:8047814/93‑1 (MQ=255)
                                             cAACGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCg                                            <  2:1476529/93‑1 (MQ=255)
                                               aCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCggg                                          >  2:6844076/1‑93 (MQ=255)
                                               aCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCggg                                          <  1:5676680/93‑1 (MQ=255)
                                               aCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCggg                                          >  1:3685904/1‑93 (MQ=255)
                                               aCGGAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCggg                                          >  1:2468296/1‑93 (MQ=255)
                                                cGGAGCTGCGAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGt                                         <  2:3276519/93‑1 (MQ=255)
                                                 ggAGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTa                                        <  1:8731687/93‑1 (MQ=255)
                                                   aGCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGaaaa                                      >  1:939383/1‑93 (MQ=255)
                                                    gCTGCGAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                     <  1:1684708/93‑1 (MQ=255)
                                                    gCTGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAAc                                     >  2:7472844/1‑93 (MQ=255)
                                                      tGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAg                                   <  1:7018371/93‑1 (MQ=255)
                                                      tGCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAg                                   <  1:4912704/93‑1 (MQ=255)
                                                       gCCAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGt                                  >  1:473320/1‑93 (MQ=255)
                                                        ccAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCGCCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGCCGGGTAAACAGtt                                 <  2:1116779/93‑1 (MQ=255)
                                                         cAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGttt                                >  2:2526592/1‑93 (MQ=255)
                                                         cAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGttt                                >  1:7983619/1‑93 (MQ=255)
                                                         cAGAGCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGttt                                <  1:374603/93‑1 (MQ=255)
                                                          agagCGTAATGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                               <  1:7410592/93‑1 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTTGGGTAAACAGTTTg                               <  1:2895583/93‑1 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                               >  1:7579126/1‑93 (MQ=255)
                                                          agagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTg                               <  1:2032320/93‑1 (MQ=255)
                                                           gagCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTgg                              >  2:6295731/1‑93 (MQ=255)
                                                            agCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTggg                             <  2:2340577/93‑1 (MQ=255)
                                                            agCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCGTTTTCCTGTCGGGTAAACATGTTggg                             >  1:5462888/1‑93 (MQ=255)
                                                             gCGTAAGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGgc                            <  2:2277537/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       >  1:1513512/1‑93 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       <  1:5138120/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAAt                       <  1:2703740/93‑1 (MQ=255)
                                                                  aGGAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGGCGGGTAAACAGTTTGGGCGCAAt                       >  2:4740705/1‑93 (MQ=255)
                                                                    gAATAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGa                     >  1:910550/1‑93 (MQ=255)
                                                                     aaTAGCCGGTTGGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGaa                    <  2:5012332/93‑1 (MQ=255)
                                                                      aTAGCTGGTTCGGTTGCGCCCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAc                   <  1:8489917/93‑1 (MQ=255)
                                                                      aTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAc                   <  2:8412908/93‑1 (MQ=255)
                                                                      aTAGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAc                   >  2:8001368/1‑93 (MQ=255)
                                                                        aGCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACgc                 >  2:6862239/1‑92 (MQ=255)
                                                                         gCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTc                >  2:2594148/1‑93 (MQ=255)
                                                                         gCCGGTTCGGTTGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTc                <  1:6080863/93‑1 (MQ=255)
                                                                           cggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGc              >  1:3974932/1‑93 (MQ=255)
                                                                           cggttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGc              >  2:762197/1‑93 (MQ=255)
                                                                             gtttggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa            >  2:6882216/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa            >  1:5811576/1‑93 (MQ=255)
                                                                             gttcggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCaa            <  2:4806006/93‑1 (MQ=255)
                                                                                 ggttGCTACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAcc        <  2:925010/93‑1 (MQ=255)
                                                                                 ggttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAcc        <  2:3902908/93‑1 (MQ=255)
                                                                                  gttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAccc       <  2:1619592/93‑1 (MQ=255)
                                                                                  gttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAccc       >  2:2840607/1‑93 (MQ=255)
                                                                                  gttGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGAccc       >  2:2459684/1‑93 (MQ=255)
                                                                                    tGCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGCCccgc     <  2:7744258/93‑1 (MQ=255)
                                                                                     gCGACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACccgcc    <  2:6766391/93‑1 (MQ=255)
                                                                                       gACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  >  1:3472328/1‑93 (MQ=255)
                                                                                       gACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  <  1:4501433/93‑1 (MQ=255)
                                                                                       gACCGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCcgccgc  <  1:2906925/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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: