Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,246,188 A→G E213E (GAA→GAG lamB → phage lambda receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,246,1880AG98.0% 648.5 / ‑0.1 205E213E (GAA→GAGlamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (97/104);  total (100/105)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 215 total aligned reads displayed.

CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGAACA  >  NC_000913/4246098‑4246280
                                                                                          |                                                                                            
 aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                           <  1:2512438/93‑1 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          <  1:3539807/93‑1 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          >  2:3469372/1‑93 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          >  1:3973098/1‑93 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          <  2:4396182/93‑1 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          <  2:411868/93‑1 (MQ=255)
  atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                          <  2:2530962/93‑1 (MQ=255)
   tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                         <  1:4360345/93‑1 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                       <  2:2900931/93‑1 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                       >  1:2914307/1‑93 (MQ=255)
     tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                       <  1:2222430/93‑1 (MQ=255)
        aTGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCg                                                                                    >  2:3509793/1‑93 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                   <  2:3117794/93‑1 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                   >  2:1156028/1‑93 (MQ=255)
         tGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGa                                                                                   <  2:1716091/93‑1 (MQ=255)
          gACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAc                                                                                  >  2:1292839/1‑93 (MQ=255)
           aCTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACt                                                                                 >  1:3328286/1‑93 (MQ=255)
           aCTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACt                                                                                 <  2:2099089/93‑1 (MQ=255)
            cTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTa                                                                                <  2:2301932/93‑1 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                               >  1:4254111/1‑93 (MQ=255)
             tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                               >  2:1099466/1‑93 (MQ=255)
              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCACCCCGGGCGGCCCATTAGAGCTGGGTGTCGACTACg                                                                              >  1:706995/1‑93 (MQ=255)
              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                              >  1:4380719/1‑93 (MQ=255)
              ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                              >  2:2476516/1‑93 (MQ=255)
               taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                             <  1:3894979/93‑1 (MQ=255)
               taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                             >  2:2846755/1‑93 (MQ=255)
               taCCAACGAAACCGCGAACGACGTTTTCGATGGGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                             <  1:1586788/93‑1 (MQ=255)
                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                            >  2:2027469/1‑93 (MQ=255)
                aCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGt                                                                            <  1:3495772/93‑1 (MQ=255)
                  cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGCCg                                                                          >  2:872055/1‑93 (MQ=255)
                   aaCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGt                                                                         <  2:380995/93‑1 (MQ=255)
                      gAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                      >  1:3734979/1‑93 (MQ=255)
                        aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCaa                                                                    >  1:3147298/1‑93 (MQ=255)
                        aaCCGCGAACGACGCTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGCCGGCACATTAGAGCTGGGTGTCGACTACGGACGTCCCac                                                                    >  2:2509339/1‑92 (MQ=255)
                           cgcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACtt                                                                 <  2:3316762/93‑1 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                                <  2:4538533/93‑1 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                                >  1:1524078/1‑93 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                                >  2:3536405/1‑93 (MQ=255)
                            gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCCGGCGGCACATTAGAGATGGGCGTCGACACCGGTCGTGCCAACTTg                                                                >  2:650237/1‑93 (MQ=255)
                             cgAACGACGTTTTCGATTTGCGTTTAGCGCAGATGGAAATCAACCCGTGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGc                                                               >  2:3627420/1‑93 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCg                                                              >  2:1532037/1‑93 (MQ=255)
                              gAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCg                                                              >  2:1791441/1‑93 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             <  2:3417947/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             >  1:688032/1‑93 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             <  1:3181438/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             <  2:2950539/93‑1 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             >  2:2996641/1‑93 (MQ=255)
                               aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                             <  1:1299858/93‑1 (MQ=255)
                                acgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTg                                                            >  1:1625670/1‑93 (MQ=255)
                                 cgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCCACTTGCGTGa                                                           >  2:4315/1‑93 (MQ=255)
                                 cgacgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGa                                                           >  1:2224731/1‑93 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                         <  1:2369002/93‑1 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                         >  1:3124381/1‑93 (MQ=255)
                                   acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                         <  1:3580651/93‑1 (MQ=255)
                                    cgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATaa                                                        <  1:1675197/93‑1 (MQ=255)
                                      ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                      >  1:4474322/1‑93 (MQ=255)
                                      ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                      >  2:420313/1‑93 (MQ=255)
                                       tttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTa                                                     <  1:3249551/93‑1 (MQ=255)
                                            aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                                >  1:552286/1‑93 (MQ=255)
                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                               <  1:3191242/93‑1 (MQ=255)
                                             tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                               >  2:2385981/1‑93 (MQ=255)
                                              gtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTg                                              >  1:4324529/1‑93 (MQ=255)
                                              gtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTg                                              <  1:4014789/93‑1 (MQ=255)
                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                             <  2:2814276/93‑1 (MQ=255)
                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                             <  1:483569/93‑1 (MQ=255)
                                               tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                             <  1:2490627/93‑1 (MQ=255)
                                                 cGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGtt                                           <  1:79844/93‑1 (MQ=255)
                                                 cGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGtt                                           >  1:2813640/1‑93 (MQ=255)
                                                   tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                         <  1:2344213/93‑1 (MQ=255)
                                                     tAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGTGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                       <  1:3064702/93‑1 (MQ=255)
                                                     tAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                       <  2:1668660/93‑1 (MQ=255)
                                                     tAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                       >  2:4494084/1‑93 (MQ=255)
                                                         gcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcgc                                   >  1:3488907/1‑93 (MQ=255)
                                                          cAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCa                                  <  1:1089282/93‑1 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                                >  1:1967571/1‑93 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                                >  1:182430/1‑93 (MQ=255)
                                                            gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                                <  1:1696824/93‑1 (MQ=255)
                                                              tGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGa                              >  1:115676/1‑93 (MQ=255)
                                                                 aaaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAg                           >  2:3051115/1‑93 (MQ=255)
                                                                  aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                          >  2:404735/1‑93 (MQ=255)
                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                         <  2:3824265/93‑1 (MQ=255)
                                                                   aTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGAc                         <  1:4644250/93‑1 (MQ=255)
                                                                    tCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACg                        >  1:3724227/1‑93 (MQ=255)
                                                                    tCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACg                        >  2:2616275/1‑93 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                      <  1:4537000/93‑1 (MQ=255)
                                                                      aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGc                      >  2:1427861/1‑93 (MQ=255)
                                                                       aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                     >  1:2434211/1‑93 (MQ=255)
                                                                       aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                     <  2:4618685/93‑1 (MQ=255)
                                                                          cGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGt                  <  1:2048444/93‑1 (MQ=255)
                                                                            ggcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTa                <  1:470965/93‑1 (MQ=255)
                                                                            ggcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTa                >  2:4726534/1‑93 (MQ=255)
                                                                                gcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCa            >  2:3694891/1‑93 (MQ=255)
                                                                                gcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCa            >  1:3384161/1‑93 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct          <  1:3802931/93‑1 (MQ=255)
                                                                                  acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct          <  1:970504/93‑1 (MQ=255)
                                                                                   caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg         >  1:769278/1‑93 (MQ=255)
                                                                                    aTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgc        <  1:3407244/93‑1 (MQ=255)
                                                                                     ttAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgct       >  1:53159/1‑93 (MQ=255)
                                                                                     ttAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgct       >  2:1682248/1‑93 (MQ=255)
                                                                                          aCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGAACa  <  1:2763909/93‑1 (MQ=255)
                                                                                          |                                                                                            
CAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGAACA  >  NC_000913/4246098‑4246280

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