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,1380AG99.6% 721.6 / ‑5.1 229I1017T (ATC→ACC) mfdtranscription‑repair coupling mutation frequency decline
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (112/116);  total (112/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 246 total aligned reads displayed.

ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227
                                                                                          |                                                                                         
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        >  1:4949152/1‑93 (MQ=255)
   aGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaaa                                                                                      <  2:5096094/93‑1 (MQ=255)
     ctctACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAAAACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    <  2:1701548/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   <  1:7841288/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   <  2:5625246/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   >  2:6055342/1‑93 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   <  1:2856858/93‑1 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  1:30479/93‑1 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  1:5043336/93‑1 (MQ=255)
         ccTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                <  2:3332038/92‑1 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTgg                                                                              >  1:6810439/1‑93 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTgg                                                                              >  1:7577046/1‑93 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGGGAAATCATCTgg                                                                              >  2:1017563/1‑93 (MQ=255)
            ttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGc                                                                             >  2:4973124/1‑93 (MQ=255)
            ttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGc                                                                             <  1:3587845/93‑1 (MQ=255)
             tGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCa                                                                            <  1:1702162/93‑1 (MQ=255)
             tGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCa                                                                            <  2:7891651/93‑1 (MQ=255)
              gATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCaa                                                                           <  2:12335/93‑1 (MQ=255)
              gATCTCTGCCAGTTCGGTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCaa                                                                           <  1:4106906/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  1:3663256/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  1:5419776/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       >  1:3347283/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  1:5239319/93‑1 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  2:3405532/93‑1 (MQ=255)
                  tctTCCAGTGCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  2:4503179/93‑1 (MQ=255)
                   ctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCg                                                                      <  1:2764054/93‑1 (MQ=255)
                     tCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGAc                                                                    <  1:7282188/93‑1 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  1:957471/1‑93 (MQ=255)
                         gTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                <  1:441817/93‑1 (MQ=255)
                         gTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                >  2:2303637/1‑93 (MQ=255)
                          ttCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCAt                                                               >  2:4806485/1‑93 (MQ=255)
                           tctttttCTTTTTGGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              <  1:7408841/90‑1 (MQ=255)
                           tCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              <  2:6476855/93‑1 (MQ=255)
                           tCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGGCAGGGGTGAAACCATCTGGCAATAGCGACGGCATc                                                              >  1:4517002/1‑93 (MQ=255)
                            tgTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATcc                                                             <  2:6678428/92‑1 (MQ=255)
                              ttttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGc                                                           >  2:4153737/1‑93 (MQ=255)
                               tttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCa                                                          <  1:5394655/93‑1 (MQ=255)
                               tttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCa                                                          >  2:1894608/1‑93 (MQ=255)
                                ttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGGGAAATCATCTGGCAATAGCGACGGCATCCGCac                                                         >  2:6511131/1‑92 (MQ=255)
                                 tCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGc                                                        <  2:5770139/93‑1 (MQ=255)
                                 tCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGc                                                        >  1:2316398/1‑93 (MQ=255)
                                   cGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTc                                                      <  1:847165/93‑1 (MQ=255)
                                     ttttGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGa                                                    <  2:4321016/93‑1 (MQ=255)
                                        tGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACtt                                                 >  2:2343051/1‑93 (MQ=255)
                                        tGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACtt                                                 <  2:7810221/93‑1 (MQ=255)
                                           cgcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTg                                              <  1:3477665/93‑1 (MQ=255)
                                            gcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGt                                             >  1:186698/1‑93 (MQ=255)
                                            gcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGt                                             >  1:7752868/1‑93 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            >  1:887202/1‑93 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            >  2:7573223/1‑93 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            <  1:2567260/93‑1 (MQ=255)
                                              tGGCAATACGTTTTTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCACCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTCttt                                           >  2:3887196/1‑93 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           >  1:640247/1‑93 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           <  2:3987415/93‑1 (MQ=255)
                                               ggCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTg                                          <  2:6967697/93‑1 (MQ=255)
                                               ggCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTg                                          <  1:6142590/93‑1 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      <  1:3766778/93‑1 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:3509070/1‑93 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:2920116/1‑93 (MQ=255)
                                                    tACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggc                                     <  2:2122920/93‑1 (MQ=255)
                                                       gattGTAGAACGCCAGACGCGTGTTCACGTCGGGGGTGAAATCATCTAGCAATACCGACGCCATTCGCAGCTCGACTTCTGTTTGctggctgg                                  >  2:710860/3‑93 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  >  1:3856293/1‑93 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  >  1:750522/1‑93 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  2:4877110/93‑1 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  1:4364858/93‑1 (MQ=255)
                                                           gTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTgag                              >  1:7144542/1‑93 (MQ=255)
                                                            tAGAACGACAGACGCGTGTTCGCGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTgaga                             <  1:7718886/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            <  2:2156326/93‑1 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCtt                         <  1:6988944/93‑1 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCtt                         >  1:4059790/1‑93 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGGGAGATCtt                         >  2:239644/1‑93 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        <  2:2588640/93‑1 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        <  2:742608/93‑1 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        <  2:1247285/93‑1 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  1:1811422/93‑1 (MQ=255)
                                                                   aCAGACGCTTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      <  1:6943464/93‑1 (MQ=255)
                                                                   aCAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      <  1:6383626/93‑1 (MQ=255)
                                                                    cAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAg                     >  1:6566121/1‑93 (MQ=255)
                                                                    cAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAg                     >  1:2210250/1‑93 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    >  1:1204120/1‑93 (MQ=255)
                                                                      gACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCg                   >  2:4433408/1‑93 (MQ=255)
                                                                      gACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCg                   >  1:581390/1‑93 (MQ=255)
                                                                        cgcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGAc                 <  1:5036138/93‑1 (MQ=255)
                                                                        cgcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGAc                 <  1:3154757/93‑1 (MQ=255)
                                                                           gtgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGc              >  2:4507523/1‑93 (MQ=255)
                                                                           gtgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGc              >  1:476932/1‑93 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             <  2:5273864/93‑1 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             <  2:7592435/93‑1 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             >  2:1317076/1‑93 (MQ=255)
                                                                            tgtTCACGTCAGGGGCGTACTCATCTGGCACTACCGACGGCATTCGCTGCTCGTCTTCTGTTTGCTGGCTGGTGAAATATTTCAGCAACgggt             >  2:6587107/1‑91 (MQ=255)
                                                                             gtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTc            >  1:4299223/1‑93 (MQ=255)
                                                                             gtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTc            <  2:2067392/93‑1 (MQ=255)
                                                                              ttCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTcc           >  2:1470132/1‑92 (MQ=255)
                                                                                  cGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTc       >  2:7324359/1‑93 (MQ=255)
                                                                                  cGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGCGAGATCTTCCAGCGACGGCTCGCGTc       >  2:285638/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      <  2:3604214/93‑1 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  1:2394386/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  2:5283147/1‑93 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  1:13365/1‑93 (MQ=255)
                                                                                       ggggTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCt  >  1:1086524/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227

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