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.1% 710.5 / NA 226I1017T (ATC→ACC) mfdtranscription‑repair coupling mutation frequency decline
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (108/115);  minor base C (0/2);  total (108/118)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

ATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGATGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCT  >  NC_000913/1170048‑1170227
                                                                                          |                                                                                         
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        >  1:4911118/1‑93 (MQ=255)
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        >  2:1632738/1‑93 (MQ=255)
 tAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGa                                                                                        <  2:2056406/93‑1 (MQ=255)
  aaGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGaa                                                                                       >  1:52489/1‑93 (MQ=255)
    gCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAat                                                                                     >  2:5053378/1‑93 (MQ=255)
     cTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatc                                                                                    <  2:5077101/93‑1 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   >  2:4119991/1‑93 (MQ=255)
      tCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatca                                                                                   >  2:1424321/1‑93 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  1:45470/93‑1 (MQ=255)
        gACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAatcatc                                                                                 <  1:3135926/93‑1 (MQ=255)
         aCTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCt                                                                                <  2:3064366/93‑1 (MQ=255)
          cTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTg                                                                               >  1:4105508/1‑93 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGGAATCATCTgg                                                                              >  2:3789483/1‑93 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTgg                                                                              <  2:5924218/93‑1 (MQ=255)
           tttGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTgg                                                                              >  1:558618/1‑93 (MQ=255)
             tGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCa                                                                            >  2:733202/1‑93 (MQ=255)
               aTCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAAt                                                                          <  2:4734586/93‑1 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  2:6033153/1‑93 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATa                                                                         >  2:4015714/1‑93 (MQ=255)
                tctcTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTACACGTCAGGGGTGAGATCATCTGGCAATa                                                                         >  1:4392480/1‑93 (MQ=255)
                 ctctTCCAGTTGGTTATCCGTTGTCGCGCTGGCAATGCGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAg                                                                        <  1:5149413/93‑1 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACTTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  1:3471039/93‑1 (MQ=255)
                  tctTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGc                                                                       <  2:2347890/93‑1 (MQ=255)
                    ttCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGa                                                                     <  1:6145256/93‑1 (MQ=255)
                     tCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGAc                                                                    >  2:3545808/1‑93 (MQ=255)
                     tCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGAc                                                                    <  2:6116106/93‑1 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   <  1:3281612/93‑1 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  2:3862910/1‑93 (MQ=255)
                      ccAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACg                                                                   >  2:6090955/1‑93 (MQ=255)
                      ccAGTTCGGTGTCCGTTTTGGCGCAGGCAATACGTTTGAAGAACGGCAGACGCGTGTTTACGTCTGGGGTGAAATCATCTGGCAATAGCGAcc                                                                   <  2:4197873/93‑2 (MQ=255)
                       cAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGGGTTCACGTCAGGGGGGAAATAATCTGGCAATATTGACgt                                                                  >  2:4903646/1‑92 (MQ=255)
                       cAGTTCGTTTTCCGTTTTGGCGCTAGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACgg                                                                  <  1:3909065/93‑1 (MQ=255)
                        aGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGc                                                                 <  1:226318/93‑1 (MQ=255)
                        aGTTCGTTTTCCGATTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAACAGCGACGGc                                                                 >  1:1879740/1‑93 (MQ=255)
                         gTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                <  2:2507818/93‑1 (MQ=255)
                         gTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCa                                                                >  1:1080904/1‑93 (MQ=255)
                           tCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATc                                                              >  1:3653263/1‑93 (MQ=255)
                            cGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATcc                                                             >  2:5446533/1‑93 (MQ=255)
                              ttttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGc                                                           >  2:4580603/1‑93 (MQ=255)
                              ttttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGGGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGc                                                           >  2:5131013/1‑93 (MQ=255)
                                ttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGAGGGCATCCGCAg                                                         >  2:5599443/1‑93 (MQ=255)
                                ttCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAg                                                         >  2:4247078/1‑93 (MQ=255)
                                 tCCGTTTTGGCGCTGGCAATACGTTTGTGGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGc                                                        <  2:1032015/93‑1 (MQ=255)
                                  ccGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCt                                                       <  1:2596498/93‑1 (MQ=255)
                                     ttttGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGa                                                    >  2:1248640/1‑93 (MQ=255)
                                         ggCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTc                                                <  2:6204585/93‑1 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            <  1:5326336/93‑1 (MQ=255)
                                             cTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAACCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGtt                                            >  2:3609331/1‑93 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           >  2:2694905/1‑93 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           <  1:3701659/93‑1 (MQ=255)
                                              tGGCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGttt                                           >  1:2633561/1‑93 (MQ=255)
                                               ggCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTg                                          >  2:2867549/1‑93 (MQ=255)
                                                gCAATACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgc                                         <  2:5985560/93‑1 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      >  2:4816018/1‑93 (MQ=255)
                                                   aTACGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctgg                                      <  2:1565849/93‑1 (MQ=255)
                                                    tACGTTTGTAGAACGACAGACACGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggc                                     >  2:1314809/1‑93 (MQ=255)
                                                     aCGTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTgctggct                                    <  2:3516943/93‑1 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  <  1:375955/93‑1 (MQ=255)
                                                       gTTTGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctgg                                  >  2:3750036/1‑93 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGctggctggg                                 >  1:2068382/1‑92 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 <  2:1772245/93‑1 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 >  1:1730272/1‑93 (MQ=255)
                                                        tttGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGt                                 >  1:3019359/1‑93 (MQ=255)
                                                          tGTAGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTga                               <  1:1275842/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            <  1:1814506/93‑1 (MQ=255)
                                                             aGAACGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGAt                            >  2:4440726/1‑93 (MQ=255)
                                                               aaCGACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCt                          <  2:3291273/93‑1 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  2:164531/93‑1 (MQ=255)
                                                                  gACAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTcc                       <  1:356988/93‑1 (MQ=255)
                                                                   aCAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      >  2:809667/1‑93 (MQ=255)
                                                                   aCAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      <  2:2161995/93‑1 (MQ=255)
                                                                   aCAGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCa                      <  1:714633/93‑1 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    <  1:1147380/93‑1 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    <  2:1166563/93‑1 (MQ=255)
                                                                     aGACGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGc                    <  2:2030767/93‑1 (MQ=255)
                                                                       aCGCGTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGa                  <  1:2013089/93‑1 (MQ=255)
                                                                        cgcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGAc                 >  2:3893234/1‑93 (MQ=255)
                                                                         gcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACg                <  1:5495746/93‑1 (MQ=255)
                                                                         gcgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACg                <  1:2239160/93‑1 (MQ=255)
                                                                          cgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACgg               >  2:2210024/1‑93 (MQ=255)
                                                                          cgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACgg               >  1:2437197/1‑93 (MQ=255)
                                                                          cgTGTTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACgg               <  2:570736/93‑1 (MQ=255)
                                                                           gtgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGc              >  1:5896965/1‑93 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             >  1:5025250/1‑93 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             <  1:4802428/93‑1 (MQ=255)
                                                                            tgtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCAACCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCt             >  1:3384000/1‑93 (MQ=255)
                                                                             gtTCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTc            <  1:3999882/93‑1 (MQ=255)
                                                                               tCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTcgc          <  1:3599548/93‑1 (MQ=255)
                                                                               tCACGTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTcgc          >  1:6242820/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  1:4138040/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  2:4981890/1‑93 (MQ=255)
                                                                                   gTCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTcc      >  1:134932/1‑93 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  2:5233249/1‑93 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     <  1:1632738/93‑1 (MQ=255)
                                                                                    tCAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     >  1:1458465/1‑93 (MQ=255)
                                                                                    gcAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCg     <  1:448390/92‑1 (MQ=255)
                                                                                     cAGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTNCCAGCGACGGCTCGCGTCCgg    >  1:2192413/1‑93 (MQ=255)
                                                                                      aGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCgg  >  2:1510962/1‑92 (MQ=255)
                                                                                      aGGGGTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGc   <  1:3045613/93‑1 (MQ=255)
                                                                                       ggggTGAAATCATCTGGCAATAGCGACGGCATCCGCAGCTCGACTTCTGTTTGCTGGCTGGTGAGATCTTCCAGCGACGGCTCGCGTCCGGCt  >  1:528160/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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: