Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,209,048 A→G intergenic (+84/‑185) b2122 → / → yehR hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,209,0480AG88.1% 604.9 / 69.2 252intergenic (+84/‑185)b2122/yehRhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/19);  new base G (102/120);  total (113/139)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 257 total aligned reads displayed.

CAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGA  >  NC_000913/2208956‑2209140
                                                                                             |                                                                                            
  ggCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGcgcg                                                                                            >  2:3528576/1‑93 (MQ=255)
   gCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCgg                                                                                           >  2:914644/1‑93 (MQ=255)
   gCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCgg                                                                                           >  2:4525079/1‑93 (MQ=255)
   gCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACgg                                                                                           <  1:1744306/93‑1 (MQ=255)
    cATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGGACATGTGTTGGGGCAATTTTGTAATACAGAAAGCGAGGCGCggg                                                                                          >  1:2633953/1‑93 (MQ=255)
    cATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAACGAGGCGCggg                                                                                          >  1:6729300/1‑93 (MQ=255)
     aTAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGt                                                                                         <  2:6561806/93‑1 (MQ=255)
      tAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTa                                                                                        >  2:5029687/1‑93 (MQ=255)
        gTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGGCAATTTCGTAATACAGAAAGCGAGGCGCGGGgaga                                                                                      >  1:7185400/1‑93 (MQ=255)
        gTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGa                                                                                      <  2:681814/93‑1 (MQ=255)
         ttGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGaa                                                                                     <  2:3870927/93‑1 (MQ=255)
         ttGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGaa                                                                                     <  1:5045827/93‑1 (MQ=255)
          tGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAg                                                                                    <  1:1291148/93‑1 (MQ=255)
          tGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAg                                                                                    <  1:3297476/93‑1 (MQ=255)
           gCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGc                                                                                   <  1:5084571/93‑1 (MQ=255)
             tGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCaa                                                                                 <  2:7246022/93‑1 (MQ=255)
             tGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCaa                                                                                 <  1:6131104/93‑1 (MQ=255)
             tGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAACGAGGCGCGGGGAGAAGCaa                                                                                 >  1:6590681/1‑93 (MQ=255)
              ggCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCCGAGGCGCGGGTAGAAGCaa                                                                                 <  1:2378196/93‑1 (MQ=255)
              ggCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAAt                                                                                <  2:5231247/93‑1 (MQ=255)
                cGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGa                                                                              <  2:2737529/93‑1 (MQ=255)
                 gAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGaa                                                                             >  2:7104998/1‑93 (MQ=255)
                 gAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAGCAATGaa                                                                             >  1:3213113/1‑93 (MQ=255)
                 gAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAGCAATGaa                                                                             >  1:7455482/1‑93 (MQ=255)
                  aaTGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAAt                                                                            <  2:1413538/93‑1 (MQ=255)
                   aTGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACCTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAACAAataata                                                                           >  1:5460582/1‑93 (MQ=255)
                    tGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATaa                                                                          <  2:5633599/93‑1 (MQ=255)
                     gACCTCACCGCGGTGGGTTTTTATCGGCTGGGACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAGCAATGAATAAc                                                                         >  1:7064897/1‑93 (MQ=255)
                     gACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAAc                                                                         <  1:2831640/93‑1 (MQ=255)
                      aCCTCACCGCGGTGGGTTTTTATCGGCTGGGACATGTGTTGGGGCAATTTTGTAATACAGAAAGCGAGGGGCGGGGAGAAGCAAAGAATAACg                                                                        >  1:5136423/1‑93 (MQ=255)
                      aCCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAGCAATGAATAACg                                                                        >  1:1665864/1‑93 (MQ=255)
                        cTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgt                                                                      <  1:5886001/93‑1 (MQ=255)
                        cTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGgt                                                                      <  2:6449723/93‑1 (MQ=255)
                         tCACCGCGGTGGGTTTTTATCGGCTGGCCCATGTGTTGGGGTAATTTCGTAATACAGAAAGCGAGGGGCGGGGAGAAGCAAAGAATAACgggg                                                                     >  1:4406540/1‑91 (MQ=255)
                         tCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgtg                                                                     <  2:4907513/93‑1 (MQ=255)
                         tCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgtg                                                                     <  2:3797087/93‑1 (MQ=255)
                         tCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgtg                                                                     >  2:769087/1‑93 (MQ=255)
                          cACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgtgt                                                                    <  1:2149386/93‑1 (MQ=255)
                          cACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGgtgt                                                                    <  2:2987747/93‑1 (MQ=255)
                            ccGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGttt                                                                  <  1:3938179/93‑1 (MQ=255)
                            ccGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGttt                                                                  <  2:699222/93‑1 (MQ=255)
                             cgcgGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGtttt                                                                 <  2:4441171/93‑1 (MQ=255)
                               cgGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTtttgc                                                               <  1:3835253/93‑1 (MQ=255)
                                ggTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgct                                                              >  2:4187889/1‑93 (MQ=255)
                                  tGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGGAGAAGCAATGAATAACGGTGTtttgcttt                                                            >  1:3111155/1‑93 (MQ=255)
                                  tGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                                                            <  2:3548495/93‑1 (MQ=255)
                                     gTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg                                                         <  1:6477704/93‑1 (MQ=255)
                                     gTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg                                                         >  2:1989715/1‑93 (MQ=255)
                                          tATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACaa                                                    <  1:5980709/93‑1 (MQ=255)
                                            tCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatt                                                  <  2:2094508/93‑1 (MQ=255)
                                            tCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatt                                                  >  2:2458956/1‑93 (MQ=255)
                                             cGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatta                                                 <  2:568768/93‑1 (MQ=255)
                                               gCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAattatt                                               <  2:2914739/93‑1 (MQ=255)
                                                 tGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTcc                                             >  1:5575707/1‑93 (MQ=255)
                                                  ggCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTccc                                            >  1:1443376/1‑93 (MQ=255)
                                                    cacaTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCa                                          <  1:3712838/93‑1 (MQ=255)
                                                    cacaTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCa                                          <  2:1505978/93‑1 (MQ=255)
                                                    cacaTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCa                                          >  1:3262714/1‑93 (MQ=255)
                                                     acaTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCAt                                         <  2:4671179/93‑1 (MQ=255)
                                                       aTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGc                                       <  2:6195010/93‑1 (MQ=255)
                                                       aTGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGc                                       <  1:6264352/93‑1 (MQ=255)
                                                          tgtTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTc                                    <  1:1183952/93‑1 (MQ=255)
                                                          tgtTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTc                                    <  1:292107/93‑1 (MQ=255)
                                                          tgtTGGGACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTc                                    <  1:6538768/93‑1 (MQ=255)
                                                               ggACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaa                               <  2:143044/93‑1 (MQ=255)
                                                               ggACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaa                               <  2:1294769/93‑1 (MQ=255)
                                                               ggACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaa                               <  1:5280303/93‑1 (MQ=255)
                                                                gACAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaac                              <  1:955833/93‑1 (MQ=255)
                                                                 aCAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaaca                             >  2:3054535/1‑93 (MQ=255)
                                                                 aCAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaaca                             >  2:3963170/1‑93 (MQ=255)
                                                                  cAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaacaa                            <  2:5482487/93‑1 (MQ=255)
                                                                  cAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaacaa                            <  1:4282644/93‑1 (MQ=255)
                                                                  cAATTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGcaacaa                            <  2:4363863/93‑1 (MQ=255)
                                                                   aaTTTCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAg                           <  1:722557/93‑1 (MQ=255)
                                                                     tttCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGa                         <  1:5721026/93‑1 (MQ=255)
                                                                      ttCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGAt                        <  2:2229508/93‑1 (MQ=255)
                                                                      ttCGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGAt                        >  1:1831480/1‑93 (MQ=255)
                                                                        cGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGa                      <  2:3343274/93‑1 (MQ=255)
                                                                        cGTAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGa                      <  1:3994944/93‑1 (MQ=255)
                                                                          tAATACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGaaa                    <  1:6832177/93‑1 (MQ=255)
                                                                           aaTACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAAt                   >  2:1170357/1‑93 (MQ=255)
                                                                            aTACAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATa                  >  2:1625287/1‑93 (MQ=255)
                                                                               cAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTg               >  1:2118286/1‑93 (MQ=255)
                                                                               cAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTg               <  2:5266943/93‑1 (MQ=255)
                                                                               cAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTg               >  1:3739377/1‑93 (MQ=255)
                                                                               cAGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTg               <  1:3333055/93‑1 (MQ=255)
                                                                                aGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGc              <  2:3668198/93‑1 (MQ=255)
                                                                                aGAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGc              <  2:670980/93‑1 (MQ=255)
                                                                                 gAAAGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCt             >  2:4606798/1‑93 (MQ=255)
                                                                                  aaaGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCtt            <  1:2493797/93‑1 (MQ=255)
                                                                                    aGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCtttt          <  1:1927669/93‑1 (MQ=255)
                                                                                    aGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCtttt          <  2:2507356/93‑1 (MQ=255)
                                                                                    aGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCtttt          <  2:5913946/93‑1 (MQ=255)
                                                                                    aGCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCtttt          <  1:730726/93‑1 (MQ=255)
                                                                                     gCGAGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCttttt         >  2:6077901/1‑93 (MQ=255)
                                                                                     gCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCttttt         <  1:2507995/93‑1 (MQ=255)
                                                                                         aGGCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACt      >  1:1119391/1‑93 (MQ=255)
                                                                                          ggCGCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTa     <  1:7007913/93‑1 (MQ=255)
                                                                                            cACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGTGAACAATTATTCCCCAGGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAag   <  1:4816170/93‑1 (MQ=255)
                                                                                             aCGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAaga  <  1:2992540/93‑1 (MQ=255)
                                                                                             |                                                                                            
CAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCCCATGCTTCAGCAACAAGGATGAAATAGTGCTTTTTACTAAGA  >  NC_000913/2208956‑2209140

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