Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,983,890 T→C T88A (ACT→GCT)  araF ← L‑arabinose‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,8900TC98.2% 332.0 / ‑2.6 108T88A (ACT→GCT) araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (53/53);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

Only 100 of 118 total aligned reads displayed.

GTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983801‑1983979
                                                                                         |                                                                                         
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  1:3205534/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  1:3205530/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        >  2:2298294/1‑93 (MQ=255)
gTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGcaca                                                                                        <  2:1402644/93‑1 (MQ=255)
 ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:4692767/1‑93 (MQ=255)
 ttAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACaa                                                                                       >  2:2650109/1‑93 (MQ=255)
  tAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACaaa                                                                                      <  2:4033835/93‑1 (MQ=255)
    aCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATa                                                                                    <  2:3951241/93‑1 (MQ=255)
    aCAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGATTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATa                                                                                    >  2:2855989/1‑93 (MQ=255)
     cAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGCTTGGGGTCCGGAGCACAAATaa                                                                                   >  1:4425494/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  2:3492267/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  2:4377632/1‑93 (MQ=255)
       aaCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACg                                                                                 >  1:4185626/1‑93 (MQ=255)
        aCTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGa                                                                                >  2:152064/1‑93 (MQ=255)
           ggTCATCCACGGCAATGACTTTCCTATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             <  2:2591969/93‑1 (MQ=255)
           ggTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAAc                                                                             <  2:2141841/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:2978843/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  1:477784/93‑1 (MQ=255)
             tCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCt                                                                           <  2:1109441/93‑1 (MQ=255)
              cATCCGCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtt                                                                          <  1:2146545/93‑1 (MQ=255)
              cATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCtt                                                                          <  2:3486055/93‑1 (MQ=255)
                tCCAAGGCAGTGACGTTTATATCGTAGCCACGCGCGTACGCGACGATGGCAGAGCCGAGTTTGGGGGCCGGAGCACAAATAACGAAACCtttt                                                                        <  2:3838531/93‑1 (MQ=255)
                 ccACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTg                                                                       <  1:2280234/93‑1 (MQ=255)
                 ccACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTg                                                                       >  2:3643963/1‑93 (MQ=255)
                  cACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgc                                                                      <  1:2316966/93‑1 (MQ=255)
                  cACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgc                                                                      <  2:2720536/93‑1 (MQ=255)
                   aCGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTgcg                                                                     <  2:1196968/93‑1 (MQ=255)
                     ggCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGcc                                                                   <  2:1945172/93‑1 (MQ=255)
                      gCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGccc                                                                  >  1:3664619/1‑92 (MQ=255)
                       cAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCAc                                                                 <  2:4658685/93‑1 (MQ=255)
                         aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               <  1:2115598/93‑1 (MQ=255)
                         aTGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTg                                                               >  2:3422356/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:2806657/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:504196/1‑93 (MQ=255)
                          tGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTgg                                                              >  1:7747/1‑93 (MQ=255)
                           gACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGc                                                             >  2:202569/1‑93 (MQ=255)
                           gACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGc                                                             <  1:1069391/93‑1 (MQ=255)
                            aCTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCa                                                            >  2:4728277/1‑93 (MQ=255)
                               ttCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGTAACCTTTTGCGCCACTGGCAGcc                                                         >  2:2583916/1‑93 (MQ=255)
                                  atatCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAgg                                                      >  1:1691377/1‑93 (MQ=255)
                                   aatCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGc                                                     <  2:3362525/92‑1 (MQ=255)
                                    atCGTAGGGACGCGCTTTCGAGTTGATGGCAGAGACGAGTTTGGGGTCCGGAGCACAAATAAAGCAACCTTTTGCGCCACTGGAAGCAAGGCt                                                    >  2:3067417/1‑93 (MQ=255)
                                    atCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCt                                                    >  1:3034109/1‑93 (MQ=255)
                                     tCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTg                                                   >  2:4134381/1‑93 (MQ=255)
                                     tCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTg                                                   <  1:4352849/93‑1 (MQ=255)
                                      cGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGt                                                  <  1:845114/93‑1 (MQ=255)
                                       gTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTc                                                 <  1:1580964/93‑1 (MQ=255)
                                          gCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGAt                                              <  2:4022986/93‑1 (MQ=255)
                                          gCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGAt                                              <  2:1155172/93‑1 (MQ=255)
                                            cACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcg                                            >  2:1016163/1‑93 (MQ=255)
                                             ccgcgCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcgc                                           <  1:1904828/92‑1 (MQ=255)
                                             aCGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGGCGATcgc                                           >  2:4307275/1‑93 (MQ=255)
                                              cgcgCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATcgcg                                          >  2:2100639/1‑93 (MQ=255)
                                                cgcTTTCGCGACGACGGCAGAGCCGAGGTTGGGGTCCGGCGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGtt                                        >  1:884595/1‑93 (MQ=255)
                                                    ttCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAAt                                    >  2:1597793/1‑93 (MQ=255)
                                                    ttCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAAt                                    <  2:3993571/93‑1 (MQ=255)
                                                      cgcgACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGt                                  <  2:2993955/93‑1 (MQ=255)
                                                       gcgACGATGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtt                                 <  2:79800/93‑1 (MQ=255)
                                                        cgacgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttt                                <  1:4743991/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  1:4595111/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  2:2252333/93‑1 (MQ=255)
                                                          acgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGttttt                              <  2:2252469/93‑1 (MQ=255)
                                                           cgaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGtttttt                             <  2:1145817/93‑1 (MQ=255)
                                                            gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            <  2:3312246/93‑1 (MQ=255)
                                                            gaTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTc                            >  2:1805699/1‑93 (MQ=255)
                                                             aTGGCAGAGCCGAGTTTTGTGTCCGTAGCACAAATAACGAAACCTTTTGCGTCTCTTGTAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCg                           >  2:4337597/1‑93 (MQ=255)
                                                             aTGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCg                           >  2:2011182/1‑93 (MQ=255)
                                                              tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          <  2:1613403/93‑1 (MQ=255)
                                                              tGGCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGc                          >  2:2354484/1‑93 (MQ=255)
                                                                gCAGAGCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCa                        <  1:3931553/93‑1 (MQ=255)
                                                                  agagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATc                      >  1:2468267/1‑93 (MQ=255)
                                                                  agagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATc                      >  2:4703321/1‑93 (MQ=255)
                                                                   gagCCGAGTTTGGGGTCCGGCGCACAAATAACGAAGCCTTTTGCGCCTCGGGCAGCCAGGCTGTCGCTCGCGTTCAATGTTTTTTCGCCCTcc                     <  2:244098/93‑1 (MQ=255)
                                                                   gagCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATcc                     >  2:1971147/1‑93 (MQ=255)
                                                                     gCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgg                   >  1:204496/1‑93 (MQ=255)
                                                                     gCCGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgg                   >  1:411457/1‑93 (MQ=255)
                                                                       cGAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggca                 <  2:1523397/93‑1 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcgggac                >  1:443283/1‑93 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcac                <  2:1656784/93‑1 (MQ=255)
                                                                        gAGTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcac                <  1:749748/93‑1 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              >  1:2166167/1‑93 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              >  2:188935/1‑93 (MQ=255)
                                                                          gTTTGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacgg              <  1:1449610/93‑1 (MQ=255)
                                                                           tttGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacggc             <  2:856603/93‑1 (MQ=255)
                                                                           tttGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCcggcacggc             >  1:769485/1‑93 (MQ=255)
                                                                             tGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCaa           <  2:104476/93‑1 (MQ=255)
                                                                             tGGGGTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCaa           <  1:4586716/93‑1 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        >  2:3262300/1‑93 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        >  1:2977106/1‑93 (MQ=255)
                                                                                ggTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCt        <  1:1143838/93‑1 (MQ=255)
                                                                                 gTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCtt       >  2:2877566/1‑93 (MQ=255)
                                                                                 gTCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCtt       >  1:1434119/1‑93 (MQ=255)
                                                                                  tCCGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTta      >  2:4565518/1‑93 (MQ=255)
                                                                                   ccGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaa     <  2:376038/93‑1 (MQ=255)
                                                                                    cGGAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  1:1144038/93‑1 (MQ=255)
                                                                                    cGGAGCACAAATAACGAAACCTGTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaat    <  1:4556031/93‑1 (MQ=255)
                                                                                     cgAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaata   <  1:2712025/92‑1 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  >  2:2158838/1‑93 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  <  2:2073150/93‑1 (MQ=255)
                                                                                      gAGCACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTtaataa  >  2:1663260/1‑93 (MQ=255)
                                                                                         |                                                                                         
GTTAACAAACTGGTCATCCACGGCAATGACTTTCATATCGTAGCCACGCGCTTTCGCGACGATGGCAGAGCCGAGTTTGGGGTCCGGAGTACAAATAACGAAACCTTTTGCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAA  >  NC_000913/1983801‑1983979

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: