Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,902,722 G→A G168G (GGC→GGT bglF ← PTS system beta‑glucosides‑specific enzyme IIABC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,902,7220GA99.4% 1262.8 / ‑5.3 358G168G (GGC→GGTbglFPTS system beta‑glucosides‑specific enzyme IIABC component
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (184/172);  total (186/172)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 382 total aligned reads displayed.

CGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTACC  >  NC_000913/3902632‑3902811
                                                                                          |                                                                                         
 gcATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTaccacca                                                                                        >  1:707094/1‑92 (MQ=255)
 gcATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccaccg                                                                                        >  1:5198645/1‑93 (MQ=255)
 gcATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccaccg                                                                                        <  1:1734263/93‑1 (MQ=255)
 gcATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccaccg                                                                                        >  2:2258008/1‑93 (MQ=255)
    tCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAg                                                                                     >  2:3305917/1‑93 (MQ=255)
     ccGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGc                                                                                    <  2:2440692/93‑1 (MQ=255)
        cTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGtt                                                                                 <  1:1630699/93‑1 (MQ=255)
        cTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGtt                                                                                 <  2:1224931/93‑1 (MQ=255)
         tttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGttt                                                                                <  2:3447511/93‑1 (MQ=255)
          ttttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTc                                                                               <  1:792544/93‑1 (MQ=255)
           tttGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTcc                                                                              <  1:2266192/93‑1 (MQ=255)
             tGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTcccc                                                                            >  1:6034047/1‑93 (MQ=255)
              gCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTccccc                                                                           >  2:2028765/1‑92 (MQ=255)
              gCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCcg                                                                           >  2:692467/1‑93 (MQ=255)
              gCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCcg                                                                           <  2:666599/93‑1 (MQ=255)
               cccGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCcgc                                                                          <  1:102894/93‑1 (MQ=255)
                 cGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCgg                                                                        <  2:6744430/93‑1 (MQ=255)
                 cGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCgg                                                                        >  1:3545892/1‑93 (MQ=255)
                  gTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGgt                                                                       >  2:2558671/1‑93 (MQ=255)
                      tcGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTAt                                                                   >  2:3598033/1‑93 (MQ=255)
                      tcGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTAt                                                                   >  2:5297496/1‑93 (MQ=255)
                        gAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATcc                                                                 <  1:6363634/93‑1 (MQ=255)
                        gAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATcc                                                                 >  1:4047974/1‑93 (MQ=255)
                         aaaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATccc                                                                >  1:7883068/1‑93 (MQ=255)
                         aaaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATccc                                                                >  2:7293664/1‑93 (MQ=255)
                         aaaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATccc                                                                <  1:450649/93‑1 (MQ=255)
                          aaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCa                                                               >  2:3117565/1‑93 (MQ=255)
                          aaGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCa                                                               >  2:6431087/1‑93 (MQ=255)
                           aGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAg                                                              <  1:6965450/93‑1 (MQ=255)
                           aGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAg                                                              <  2:4728934/93‑1 (MQ=255)
                             cAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGa                                                            <  1:6816987/93‑1 (MQ=255)
                              aGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGAt                                                           >  2:6629020/1‑93 (MQ=255)
                              aGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGAt                                                           <  1:2907650/93‑1 (MQ=255)
                               gTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGAtt                                                          >  2:6513296/1‑93 (MQ=255)
                                tCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTa                                                         >  1:535883/1‑93 (MQ=255)
                                     aTTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg                                                    >  2:5181393/1‑93 (MQ=255)
                                      ttAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGga                                                   >  1:4274958/1‑93 (MQ=255)
                                      ttAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGga                                                   >  2:3221510/1‑93 (MQ=255)
                                        aaTGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaag                                                 >  1:3882660/1‑93 (MQ=255)
                                          tGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaagaa                                               >  1:6180023/1‑93 (MQ=255)
                                          tGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaagaa                                               >  1:5399162/1‑93 (MQ=255)
                                          tGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaagaa                                               <  2:2761996/93‑1 (MQ=255)
                                           ggATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAAc                                              <  1:2776039/93‑1 (MQ=255)
                                            gATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAAcc                                             >  2:1787551/1‑93 (MQ=255)
                                             aTGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCa                                            >  1:2172597/1‑93 (MQ=255)
                                               gCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCaaa                                          >  1:2507569/1‑93 (MQ=255)
                                                 aCTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAAc                                        >  2:4658883/1‑93 (MQ=255)
                                                  cTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACa                                       >  1:2350374/1‑93 (MQ=255)
                                                  cTAACGCTCCACCAATCACCATGGCAGTAAATGGATCACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACa                                       >  1:2417000/1‑93 (MQ=255)
                                                   tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACaa                                      >  1:4965365/1‑93 (MQ=255)
                                                   tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACaa                                      >  1:6479080/1‑93 (MQ=255)
                                                   tAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACaa                                      <  2:1437561/93‑1 (MQ=255)
                                                    aaCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAg                                     >  1:70221/1‑93 (MQ=255)
                                                     aCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAgg                                    <  1:3725763/93‑1 (MQ=255)
                                                     aCGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTAACCCAGGATTATCGGGAAGAACCAAAACAAgg                                    <  1:3693097/93‑1 (MQ=255)
                                                      cGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGc                                   <  2:4463722/93‑1 (MQ=255)
                                                      cGCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGc                                   <  1:525315/93‑1 (MQ=255)
                                                       gCTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCa                                  >  1:56381/1‑93 (MQ=255)
                                                        cTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCAt                                 >  1:5901020/1‑93 (MQ=255)
                                                        cTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCAt                                 <  1:4099981/93‑1 (MQ=255)
                                                        cTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCAt                                 <  2:1099977/93‑1 (MQ=255)
                                                        cTCCACCAATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGAAt                                 <  2:7017484/93‑1 (MQ=255)
                                                          ccaccaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCa                               <  1:1349533/93‑1 (MQ=255)
                                                           caccaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAGCAAGGCATCAc                              <  1:7420350/93‑1 (MQ=255)
                                                           caccaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCAc                              >  2:42815/1‑93 (MQ=255)
                                                            accaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACt                             >  2:846858/1‑93 (MQ=255)
                                                            accaATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCACAGGATTATCGGGAAGAACCAAAACAAGGCATCACt                             >  2:4595766/1‑93 (MQ=255)
                                                              caATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTgg                           <  2:4328784/93‑1 (MQ=255)
                                                              caATCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTgg                           >  1:628442/1‑93 (MQ=255)
                                                               aaTCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgc                          >  2:1906311/1‑93 (MQ=255)
                                                                aTCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgcg                         >  2:6204465/1‑93 (MQ=255)
                                                                 tCACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGgcgc                        >  1:7666575/1‑93 (MQ=255)
                                                                  cACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCt                       <  2:6034046/93‑1 (MQ=255)
                                                                  cACCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCt                       >  1:2025741/1‑93 (MQ=255)
                                                                   aCCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTg                      <  2:7874143/93‑1 (MQ=255)
                                                                   aCCATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTg                      >  1:4609171/1‑93 (MQ=255)
                                                                    ccATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGa                     <  1:1896887/93‑1 (MQ=255)
                                                                     cATGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGaa                    <  1:7295104/93‑1 (MQ=255)
                                                                      aTGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAAt                   >  1:2702344/1‑93 (MQ=255)
                                                                       tGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATa                  <  2:3869147/93‑1 (MQ=255)
                                                                       tGGCAGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGAGCTGAATa                  >  2:2948/1‑93 (MQ=255)
                                                                           aGTAAATGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAAt              >  1:929939/1‑93 (MQ=255)
                                                                              aaaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaa           <  1:3237010/93‑1 (MQ=255)
                                                                               aaTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTaaa          >  2:5001449/1‑93 (MQ=255)
                                                                                aTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAAt         >  2:4872471/1‑93 (MQ=255)
                                                                                aTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAAt         <  1:7568940/93‑1 (MQ=255)
                                                                                aTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAAt         >  2:6856618/1‑93 (MQ=255)
                                                                                aTGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAAt         >  1:2399835/1‑93 (MQ=255)
                                                                                 tGGATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCAATGGCGCTGAATAAAATTAAATa        <  1:1392917/93‑1 (MQ=255)
                                                                                  tgATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACAAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa       <  2:5837293/92‑1 (MQ=255)
                                                                                  ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa       >  2:5440915/1‑93 (MQ=255)
                                                                                  ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa       >  1:671061/1‑93 (MQ=255)
                                                                                  ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa       >  2:4128824/1‑93 (MQ=255)
                                                                                  ggATTACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATaa       <  2:7641623/93‑1 (MQ=255)
                                                                                     ttACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTa    >  1:308284/1‑93 (MQ=255)
                                                                                     ttACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTa    <  1:4621150/93‑1 (MQ=255)
                                                                                     ttACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTa    >  2:1610467/1‑93 (MQ=255)
                                                                                     ttACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTa    <  1:116489/93‑1 (MQ=255)
                                                                                      tACCACCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTAc   <  2:7141953/93‑1 (MQ=255)
                                                                                       accaccGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTAcc  <  2:7473670/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGCATCACTGGCGCTGAATAAAATTAAATAAGTACC  >  NC_000913/3902632‑3902811

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