Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,170,138 A→G I1017T (ATC→ACC)  mfd ← transcription‑repair coupling mutation frequency decline

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,1380AG100.0% 353.7 / NA 112I1017T (ATC→ACC) mfdtranscription‑repair coupling mutation frequency decline
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (54/58);  total (54/58)

Only 100 of 118 total aligned reads displayed.

ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227
                                                                                          |                                                                                         
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        >  1:1084586/1‑93 (MQ=255)
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        >  1:2632145/1‑93 (MQ=255)
  aaGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaa                                                                                       >  2:4124765/1‑93 (MQ=255)
   aGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaaa                                                                                      >  1:1124492/1‑93 (MQ=255)
   aGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaaa                                                                                      >  2:2264584/1‑93 (MQ=255)
     cTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGTCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    >  1:3854672/1‑93 (MQ=255)
     cTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    <  1:3769631/93‑1 (MQ=255)
     cTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    <  1:3433311/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   <  2:4283834/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTCTGGCGCTGTCAATACGTTTGTAGAACGACAGACGCGTTTTCACGTCAGGGGTGAAatca                                                                                   <  2:3514840/93‑1 (MQ=255)
       cGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcat                                                                                  >  1:3849184/1‑93 (MQ=255)
       cGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcat                                                                                  >  2:436456/1‑93 (MQ=255)
       cGACTTTGATCTCTTCCAATTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcat                                                                                  <  2:3470813/93‑1 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                >  1:598172/1‑93 (MQ=255)
          cTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTg                                                                               >  1:489339/1‑93 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTgg                                                                              >  1:1315539/1‑93 (MQ=255)
            ttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGc                                                                             <  1:3094698/93‑1 (MQ=255)
              gATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCaa                                                                           <  2:302924/93‑1 (MQ=255)
              gATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCaa                                                                           <  1:236517/93‑1 (MQ=255)
                tctcTTCGAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACGGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  2:1968279/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  1:3323971/1‑93 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  2:2928490/1‑93 (MQ=255)
                 ctctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAg                                                                        >  2:3883023/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  1:2776195/93‑1 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  1:111655/93‑1 (MQ=255)
                   ctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCg                                                                      >  1:3744340/1‑93 (MQ=255)
                   ctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAAAACGACAGACGCGTGTTCACGTCATGGGTGAACTCATCTCGCAATAGCg                                                                      >  2:1165075/1‑93 (MQ=255)
                    ttCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGa                                                                     >  2:1572680/1‑93 (MQ=255)
                     tCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGAc                                                                    <  2:723055/93‑1 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   <  2:2188207/93‑1 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  2:1285542/1‑93 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  2:4675199/1‑93 (MQ=255)
                       cAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACgg                                                                  >  2:2794108/1‑93 (MQ=255)
                        aGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGc                                                                 >  1:4323219/1‑93 (MQ=255)
                         gTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                >  2:2741160/1‑93 (MQ=255)
                           tCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              >  2:3151508/1‑93 (MQ=255)
                           tCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              <  2:2689435/93‑1 (MQ=255)
                            cGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATcc                                                             <  2:1870631/93‑1 (MQ=255)
                             gTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCAGCTGGCAATAGCGACGGCATCCg                                                            <  2:2374145/93‑1 (MQ=255)
                              ttttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGc                                                           >  1:420254/1‑93 (MQ=255)
                              ttttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGc                                                           >  1:2535458/1‑93 (MQ=255)
                                ttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAg                                                         <  2:1927206/93‑1 (MQ=255)
                                ttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAAGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAg                                                         <  1:2510657/93‑1 (MQ=255)
                                 tCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGc                                                        <  2:360952/93‑1 (MQ=255)
                                    gTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGGAATAGCGACGGCATCCGCAGCTCg                                                     >  1:1623414/1‑93 (MQ=255)
                                    gTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCg                                                     >  2:879461/1‑93 (MQ=255)
                                     ttttGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGGTCACGTCAGGGGTGGAATCATCTGGCAATAGCGACGGCATCCGCAGCTCgc                                                    >  2:3233660/1‑92 (MQ=255)
                                      tttGGCGCTGGGAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGAc                                                   <  2:1237010/93‑1 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            <  1:4445748/93‑1 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           >  2:2474305/1‑93 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           <  2:1601116/93‑1 (MQ=255)
                                               ggCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTg                                          >  2:2517756/1‑93 (MQ=255)
                                               ggCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTg                                          <  2:1853716/93‑1 (MQ=255)
                                                gCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgc                                         <  1:1016386/93‑1 (MQ=255)
                                                gCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgc                                         <  1:2763402/93‑1 (MQ=255)
                                                 cAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgct                                        >  1:3236009/1‑93 (MQ=255)
                                                 cAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgct                                        <  2:144117/93‑1 (MQ=255)
                                                 cAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgct                                        <  1:795780/93‑1 (MQ=255)
                                                  aaTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctg                                       >  2:3519950/1‑93 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:3422176/1‑93 (MQ=255)
                                                    tACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggc                                     >  1:3248678/1‑93 (MQ=255)
                                                      cGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggctg                                   <  1:1019213/93‑1 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctggg                                 >  1:1529984/1‑92 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 >  2:4100425/1‑93 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 <  1:3965546/93‑1 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 >  2:3849405/1‑93 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 <  1:1300845/93‑1 (MQ=255)
                                                         ttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTg                                >  2:627676/1‑93 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  2:912574/93‑1 (MQ=255)
                                                           gTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTgag                              <  1:4251110/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            <  1:2404449/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            >  1:320296/1‑93 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            <  1:4241307/93‑1 (MQ=255)
                                                              gAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATc                           <  1:782477/93‑1 (MQ=255)
                                                              gAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATc                           >  1:4673585/1‑93 (MQ=255)
                                                               aaCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCt                          <  2:3159240/93‑1 (MQ=255)
                                                               aaCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCt                          <  2:4084718/93‑1 (MQ=255)
                                                               aaCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGGGAGATCt                          >  1:3841836/1‑93 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCtt                         <  1:2844748/93‑1 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGCCATCCGCAGCTCGGCTTCTGATTGCTGGCTGGTGAAATcgt                         >  2:3818983/1‑91 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  2:1682962/93‑1 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  2:4038936/93‑1 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  1:1565108/93‑1 (MQ=255)
                                                                   aCAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      <  1:1847044/93‑1 (MQ=255)
                                                                    cAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAg                     >  2:1204445/1‑93 (MQ=255)
                                                                       aCGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGa                  <  1:2952272/93‑1 (MQ=255)
                                                                         gcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACg                <  2:1432384/93‑1 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             <  1:665401/93‑1 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             >  2:114503/1‑93 (MQ=255)
                                                                             gtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTc            <  2:3033706/93‑1 (MQ=255)
                                                                              ttCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTcg           >  2:1644875/1‑93 (MQ=255)
                                                                                  cGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTc       <  2:3263338/93‑1 (MQ=255)
                                                                                  cGTCAGGGGTAAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTc       <  2:3552838/93‑1 (MQ=255)
                                                                                   gTCAGGTGTTAACTCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTGCTGTTTGCGCGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      <  2:2339063/93‑1 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  1:1557425/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  2:2472777/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      <  2:4220310/93‑1 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  1:2179292/1‑93 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  2:4000091/1‑93 (MQ=255)
                                                                                       ggggTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCt  >  2:287173/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227

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