Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,198,118 A→G D87G (GAC→GGC)  hupA → DNA‑binding protein HU‑alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,1180AG99.6% 822.7 / ‑5.1 261D87G (GAC→GGC) hupADNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (114/145);  total (114/147)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

Only 100 of 265 total aligned reads displayed.

AGCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198028‑4198207
                                                                                          |                                                                                         
 gCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGgcgc                                                                                        >  2:3394485/1‑93 (MQ=255)
    tACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGt                                                                                     >  1:2619694/1‑93 (MQ=255)
      cTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTta                                                                                   >  2:2929065/1‑93 (MQ=255)
      cTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTta                                                                                   >  1:1675368/1‑93 (MQ=255)
         gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagt                                                                                <  2:464444/93‑1 (MQ=255)
         gCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaaga                                                                                >  1:1294065/1‑92 (MQ=255)
           cGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaa                                                                              <  1:4598806/93‑1 (MQ=255)
           cGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaa                                                                              <  2:3774256/93‑1 (MQ=255)
           cGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTtaagtaa                                                                              <  2:3090799/93‑1 (MQ=255)
            gCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTtaagtaag                                                                             >  2:2594949/1‑93 (MQ=255)
             cAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGa                                                                            <  1:1111551/93‑1 (MQ=255)
              aaCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAt                                                                           <  2:3224348/93‑1 (MQ=255)
              aaCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAt                                                                           <  2:4657305/93‑1 (MQ=255)
              aaCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAt                                                                           >  1:3246477/1‑93 (MQ=255)
               aCCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGAtt                                                                          <  1:26963/93‑1 (MQ=255)
                 ccGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGc                                                                        >  2:1710503/1‑93 (MQ=255)
                  cGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCg                                                                       <  2:3001847/93‑1 (MQ=255)
                   gCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGt                                                                      <  2:3629680/93‑1 (MQ=255)
                     aGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTgg                                                                    >  2:3551532/1‑93 (MQ=255)
                     aGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTgg                                                                    >  2:2229933/1‑93 (MQ=255)
                        ccGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAg                                                                 <  2:4495750/93‑1 (MQ=255)
                         cGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGt                                                                <  2:3709589/93‑1 (MQ=255)
                         cGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGt                                                                >  2:3243299/1‑93 (MQ=255)
                          ggTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTg                                                               >  1:2317343/1‑93 (MQ=255)
                            tAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGaa                                                             >  1:580892/1‑93 (MQ=255)
                              aaGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACa                                                           >  1:2930394/1‑93 (MQ=255)
                               aGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAg                                                          >  2:343397/1‑93 (MQ=255)
                               aGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAg                                                          <  1:2954701/93‑1 (MQ=255)
                                gAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGt                                                         >  2:1893737/1‑93 (MQ=255)
                                gAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGt                                                         <  2:4037458/93‑1 (MQ=255)
                                 aaaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtt                                                        >  1:601893/1‑93 (MQ=255)
                                  aaTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGttt                                                       <  2:117833/93‑1 (MQ=255)
                                   aTCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGtttt                                                      >  1:2205184/1‑93 (MQ=255)
                                      aaaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAAc                                                   >  1:522130/1‑93 (MQ=255)
                                       aaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACg                                                  >  1:4558772/1‑93 (MQ=255)
                                       aaaTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACg                                                  <  1:3873824/93‑1 (MQ=255)
                                         aTTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGaa                                                <  1:3025824/93‑1 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               <  2:585144/93‑1 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               <  1:862460/93‑1 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               <  2:2481863/93‑1 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               >  1:3057457/1‑93 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               >  1:4136688/1‑93 (MQ=255)
                                          ttGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTGAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAg                                               <  1:1790115/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  1:2783709/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:3194890/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:1636862/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  1:924643/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:4257570/93‑1 (MQ=255)
                                           tGCCGCAGCTAACGTACCGGCATGTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAgg                                              <  2:2110424/93‑1 (MQ=255)
                                            gCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAggg                                             <  1:502854/93‑1 (MQ=255)
                                                 aGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtggt                                        <  1:2646767/93‑1 (MQ=255)
                                                  gCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGggtgggt                                       >  1:2104697/1‑91 (MQ=255)
                                                  gCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGtt                                       >  1:3739199/1‑93 (MQ=255)
                                                  gCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGtt                                       <  2:2000971/93‑1 (MQ=255)
                                                   cTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGttt                                      <  2:518030/93‑1 (MQ=255)
                                                   cTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGGGGttt                                      >  1:4532868/1‑93 (MQ=255)
                                                    tAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTc                                     >  1:197646/1‑93 (MQ=255)
                                                     aaCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCa                                    >  1:1326092/1‑93 (MQ=255)
                                                      aCGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAc                                   >  1:1190039/1‑93 (MQ=255)
                                                        gTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAccc                                 <  2:2536889/93‑1 (MQ=255)
                                                        gTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAccc                                 <  2:2405148/93‑1 (MQ=255)
                                                         tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                <  1:4241610/93‑1 (MQ=255)
                                                         tACCGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCAcccc                                <  2:359702/93‑1 (MQ=255)
                                                           ccGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtt                              >  2:3327539/1‑93 (MQ=255)
                                                            cGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCttt                             <  2:2197552/93‑1 (MQ=255)
                                                            cGGCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCCTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCttt                             <  1:2505505/93‑1 (MQ=255)
                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  2:4280943/93‑1 (MQ=255)
                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            <  1:3402750/93‑1 (MQ=255)
                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            >  2:1489663/1‑93 (MQ=255)
                                                             ggCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCtttt                            >  1:1458339/1‑93 (MQ=255)
                                                              gCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTg                           <  1:3284469/93‑1 (MQ=255)
                                                              gCATTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTg                           >  1:3380097/1‑93 (MQ=255)
                                                                aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACATTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTc                         <  2:1877497/93‑1 (MQ=255)
                                                                aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGTGGAGGGTTCACCCCTTgtgtt                         >  1:1404698/1‑92 (MQ=255)
                                                                aTTTGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGtt                         >  2:1405593/1‑92 (MQ=255)
                                                                 tttGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCt                        >  1:4429746/1‑93 (MQ=255)
                                                                   tGTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCttt                      <  2:3421400/93‑1 (MQ=255)
                                                                    gTTTCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTc                     <  1:2190805/93‑1 (MQ=255)
                                                                     tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  2:835201/93‑1 (MQ=255)
                                                                     tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    >  1:2721440/1‑93 (MQ=255)
                                                                     tttCTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCt                    <  2:1219267/93‑1 (MQ=255)
                                                                        cTGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGc                 >  1:3445375/1‑93 (MQ=255)
                                                                        cTGGCAAGACACAGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTCTTAACGAAGGGGTGTTTTCACCCCTTTTGTCTTTCTGGc                 <  2:16076/93‑1 (MQ=255)
                                                                         tGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCg                <  1:1017368/93‑1 (MQ=255)
                                                                         tGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCg                >  2:3960372/1‑93 (MQ=255)
                                                                         tGGCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCg                <  2:1732827/93‑1 (MQ=255)
                                                                          ggCAAGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGt               <  1:657617/93‑1 (MQ=255)
                                                                             aaGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGa            >  2:1480010/1‑93 (MQ=255)
                                                                              aGGCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGAt           <  1:684552/93‑1 (MQ=255)
                                                                                gCACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCa         <  2:2641504/93‑1 (MQ=255)
                                                                                 cACTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCAt        <  1:780663/93‑1 (MQ=255)
                                                                                   cTGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTg      <  2:1982281/93‑1 (MQ=255)
                                                                                    tGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa     >  2:2745032/1‑93 (MQ=255)
                                                                                    tGAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGa     <  1:963535/93‑1 (MQ=255)
                                                                                     gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    <  2:700423/93‑1 (MQ=255)
                                                                                     gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    <  1:2872688/93‑1 (MQ=255)
                                                                                     gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    >  2:125162/1‑93 (MQ=255)
                                                                                     gAAAGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGAt    <  2:2449372/93‑1 (MQ=255)
                                                                                      aaaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATg   >  1:3907228/1‑93 (MQ=255)
                                                                                       aaGGCGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATgc  <  1:4091260/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCGTACTGGCCGCAACCCGCAGACCGGTAAAGAAATCAAAATTGCCGCAGCTAACGTACCGGCATTTGTTTCTGGCAAGGCACTGAAAGACGCAGTTAAGTAAGATTGCGTGGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGC  >  NC_000913/4198028‑4198207

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: