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,1380AG97.0% 412.4 / 0.9 134I1017T (ATC→ACC) mfdtranscription‑repair coupling mutation frequency decline
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (58/71);  total (61/73)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 148 total aligned reads displayed.

ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227
                                                                                          |                                                                                         
atAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTg                                                                                         >  1:2628644/1‑93 (MQ=255)
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        <  1:3158396/93‑1 (MQ=255)
   aGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaaa                                                                                      >  2:753474/1‑93 (MQ=255)
   aGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaaa                                                                                      <  2:2290264/93‑1 (MQ=255)
    gCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAat                                                                                     <  2:1559472/93‑1 (MQ=255)
     cTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    <  1:1196472/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   <  2:3275341/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   >  2:2642979/1‑93 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGCTGAAatca                                                                                   >  1:2714327/1‑93 (MQ=255)
       cGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcat                                                                                  >  2:2444679/1‑93 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  2:1687370/93‑1 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  2:1998746/93‑1 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                >  1:2851001/1‑93 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                >  2:4283106/1‑93 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                <  2:3407515/93‑1 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTGTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                <  2:2271888/93‑1 (MQ=255)
          cTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTg                                                                               <  1:1511647/93‑1 (MQ=255)
             tGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCa                                                                            <  1:908711/93‑1 (MQ=255)
               aTCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAAt                                                                          <  2:3700650/93‑1 (MQ=255)
               aTCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAAt                                                                          <  2:1847677/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  2:3587831/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  1:3958766/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  2:618919/1‑93 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  1:1344587/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGCCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  1:749304/1‑93 (MQ=255)
                tctcGTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         <  2:2318855/93‑1 (MQ=255)
                 ctctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGGGTTCACGTCAGGGGGGAAATCATCTGGCAATAg                                                                        >  2:2225868/1‑93 (MQ=255)
                 ctctTCCAGTTCGTTTTCCGTTTTCGCGCTGGCAATACGTTTGCAGTACGACAGACGCGTGTTCACGTCATGGGTGAAATCATCCTGCAACAg                                                                        >  2:1271852/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  2:1087154/93‑1 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       >  2:4619290/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       >  1:1835490/1‑93 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       >  1:1908945/1‑93 (MQ=255)
                    ttCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGa                                                                     <  1:3281985/93‑1 (MQ=255)
                    ttCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGa                                                                     >  2:4641937/1‑93 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  1:3295955/1‑93 (MQ=255)
                       cAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACgg                                                                  >  2:1590598/1‑93 (MQ=255)
                        aGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAAGGGTGAAATCATCTGGCAATAGCGACGGc                                                                 >  1:2833815/1‑93 (MQ=255)
                         gTTCGTTTTCTGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                >  1:3831364/1‑93 (MQ=255)
                          ttCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCAt                                                               <  2:1501267/93‑1 (MQ=255)
                           tggTTTTCCTTGTTGGCGCTGGCATTAGTTTTGTAGAACGACAGGCGCGTTTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              <  2:4134129/91‑1 (MQ=255)
                            cGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATcc                                                             <  1:4623727/93‑1 (MQ=255)
                             gTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCg                                                            >  2:1973467/1‑93 (MQ=255)
                             gTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCg                                                            <  2:1534745/93‑1 (MQ=255)
                                 tCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGc                                                        <  2:3672880/93‑1 (MQ=255)
                                  ccGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCt                                                       <  2:3827052/93‑1 (MQ=255)
                                  ccGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCt                                                       >  2:958619/1‑93 (MQ=255)
                                   cGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTc                                                      <  2:1054037/93‑1 (MQ=255)
                                    gTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCg                                                     >  2:1918634/1‑93 (MQ=255)
                                    gTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCg                                                     <  1:4264915/93‑1 (MQ=255)
                                     ttttGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGa                                                    >  2:1461385/1‑93 (MQ=255)
                                        tGGCGCTTGCAATTCATTTTTAGAACGACGGACTGGTATTCAGTCCAAGCTTGAAATCATCTGGTAATAGCGACGGCATCCGCAGCTCTACtt                                                 <  2:2478924/93‑1 (MQ=255)
                                         ggCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTc                                                <  1:4460494/93‑1 (MQ=255)
                                          gcgcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCt                                               >  1:3318208/1‑93 (MQ=255)
                                           cgcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTg                                              <  2:531457/93‑1 (MQ=255)
                                           cgcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTg                                              <  2:4277696/93‑1 (MQ=255)
                                            gcTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGt                                             >  1:2633093/1‑93 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            <  1:2087928/93‑1 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            <  2:3903186/93‑1 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            >  1:2362258/1‑93 (MQ=255)
                                               ggCAACACGTTTGTAGAACGACAGACGCGGGTGCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCGCGACTTCTGTTTg                                          >  1:4212645/1‑93 (MQ=255)
                                                 cAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGTTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgct                                        <  1:1195749/93‑1 (MQ=255)
                                                 cAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgct                                        <  1:4743797/93‑1 (MQ=255)
                                                  aaTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctg                                       >  1:2485765/1‑93 (MQ=255)
                                                  aaTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctg                                       >  2:4532069/1‑93 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCCGGGGGGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgcttt                                      >  2:4257514/1‑91 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:1228259/1‑93 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:2431037/1‑93 (MQ=255)
                                                     aCGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggct                                    <  1:447599/93‑1 (MQ=255)
                                                      cGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggctg                                   <  2:2449535/93‑1 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  >  1:3505461/1‑93 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGTCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  <  2:2794026/93‑1 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 >  1:347978/1‑93 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 <  2:3434749/93‑1 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGTGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  1:610365/93‑1 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  2:1711481/93‑1 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  2:3677828/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            >  2:4477889/1‑93 (MQ=255)
                                                              gAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATc                           >  1:118072/1‑93 (MQ=255)
                                                              gAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATc                           >  1:2679788/1‑93 (MQ=255)
                                                                aCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGGGAGATCtt                         >  1:1718833/1‑93 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        >  2:3361836/1‑93 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        >  1:3996676/1‑93 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        <  2:2098501/93‑1 (MQ=255)
                                                                 cGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTc                        <  2:2145847/93‑1 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    >  2:2941001/1‑93 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    <  1:1613461/93‑1 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    <  1:1544073/93‑1 (MQ=255)
                                                                       aCGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGa                  <  1:670751/93‑1 (MQ=255)
                                                                           gtgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGc              <  1:1351413/93‑1 (MQ=255)
                                                                           gtgtCCACGTCAGGAGTGAAATCATCTGGCCATAGCGACGGCATCCGCAGCTCGCCTTCTGTTTGCTGGCTGGCGAGATCTTCCAGCGACGGc              <  2:1665969/93‑1 (MQ=255)
                                                                             gtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTc            >  2:135384/1‑93 (MQ=255)
                                                                                 aCGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTcgcgc        >  2:2630284/1‑92 (MQ=255)
                                                                                  cGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTc       <  1:3147381/93‑1 (MQ=255)
                                                                                  cGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTc       <  1:3825081/93‑1 (MQ=255)
                                                                                   gtgaGGGGTGAAATCATCTGGCAATAGCGACGGTATGCGCAGCTCGACTTCTGTTTGCTGGGTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      <  1:3668594/90‑1 (MQ=255)
                                                                                   gtaaGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      <  1:3668595/90‑1 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  2:3262777/1‑93 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACCGCTCGCGTCCg     <  2:4676954/93‑1 (MQ=255)
                                                                                      aGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGc   <  1:2250395/93‑1 (MQ=255)
                                                                                       ggggTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCt  <  2:803653/93‑1 (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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑