Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,356,775 C→T A279T (GCC→ACC)  cadB ← transport of lysine/cadaverine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,7750CT98.2% 554.4 / ‑2.4 168A279T (GCC→ACC) cadBtransport of lysine/cadaverine
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (81/82);  total (84/84)
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.71e-01

Only 100 of 180 total aligned reads displayed.

CTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAG  >  NC_000913/4356685‑4356864
                                                                                          |                                                                                         
cTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTgg                                                                                         <  2:1804744/93‑1 (MQ=255)
 tttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTggg                                                                                        >  2:2833113/1‑92 (MQ=255)
  ttCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGgtggtg                                                                                       >  1:2393470/1‑93 (MQ=255)
   tCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGa                                                                                      >  2:4284272/1‑93 (MQ=255)
   tCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGa                                                                                      <  2:2430503/93‑1 (MQ=255)
      ggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATg                                                                                   >  1:907948/1‑93 (MQ=255)
      ggAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATg                                                                                   >  2:2133519/1‑93 (MQ=255)
       gAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGc                                                                                  <  1:2920102/93‑1 (MQ=255)
       gAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGc                                                                                  <  1:3545272/93‑1 (MQ=255)
       gAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGCGAATGc                                                                                  >  2:4230388/1‑93 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               >  2:559042/1‑93 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               <  2:437219/93‑1 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               >  2:3234934/1‑93 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               >  1:3660171/1‑93 (MQ=255)
          gTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGa                                                                               >  1:3974871/1‑93 (MQ=255)
            tACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGGCCAGAGAAGTCAGGCACGCAAAGGTGGTGCATGCAGaaa                                                                             >  1:450073/1‑93 (MQ=255)
              ccGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAAcc                                                                           >  1:4340039/1‑93 (MQ=255)
              ccGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAAcc                                                                           >  1:2788699/1‑93 (MQ=255)
               cgtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCa                                                                          <  1:2061773/93‑1 (MQ=255)
               cgtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCa                                                                          <  2:2364478/93‑1 (MQ=255)
                gtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAg                                                                         >  2:1604751/1‑93 (MQ=255)
                gtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAg                                                                         >  1:3099724/1‑93 (MQ=255)
                gtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAg                                                                         <  1:3558979/93‑1 (MQ=255)
                gtcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAg                                                                         >  1:1766953/1‑93 (MQ=255)
                 tcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAACCCAgc                                                                        >  1:3014354/1‑93 (MQ=255)
                 tcgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgc                                                                        >  1:3016725/1‑93 (MQ=255)
                   gtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAgcgg                                                                      <  2:3258231/93‑1 (MQ=255)
                       gCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCa                                                                  >  2:735401/1‑93 (MQ=255)
                           cACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGccc                                                              <  1:4188299/93‑1 (MQ=255)
                           cACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGccc                                                              <  1:3978431/93‑1 (MQ=255)
                           cACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGccc                                                              <  2:3006848/93‑1 (MQ=255)
                            ccGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  1:1062197/92‑1 (MQ=255)
                            ccGTACACCTGCCTGAACTACCATCATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  1:2955012/92‑1 (MQ=255)
                            aCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  2:1821931/93‑1 (MQ=255)
                            aCGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCa                                                             <  1:769021/93‑1 (MQ=255)
                             cGTACACCTGCCTGGCCTACCAACATCATCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCCGt                                                           >  1:1776175/1‑93 (MQ=255)
                               tACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGtt                                                          <  2:1392920/93‑1 (MQ=255)
                                agacCTGCCTTGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTa                                                         <  2:2186629/91‑1 (MQ=255)
                                acacCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTa                                                         >  2:1882311/1‑93 (MQ=255)
                                acacCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTa                                                         <  2:2691588/93‑1 (MQ=255)
                                  acCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTAcc                                                       >  1:3345352/1‑93 (MQ=255)
                                   cctgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCg                                                      >  2:1585096/1‑93 (MQ=255)
                                    ctgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGa                                                     <  2:1911774/93‑1 (MQ=255)
                                     tgcctgGCCTACCAACATCATCCAGGATCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAg                                                    <  2:1501450/93‑1 (MQ=255)
                                     cgcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAg                                                    <  1:3186590/92‑1 (MQ=255)
                                      gcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAgg                                                   <  1:3183479/93‑1 (MQ=255)
                                      gcctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGGGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAgg                                                   >  1:575281/1‑93 (MQ=255)
                                       cctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGa                                                  <  1:2858857/93‑1 (MQ=255)
                                        ctgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAt                                                 <  1:4070271/93‑1 (MQ=255)
                                        ctgGCCTACCAACATCATCCAGGAGCCCAGACAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAt                                                 <  1:433949/93‑1 (MQ=255)
                                         tgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATa                                                <  1:527801/93‑1 (MQ=255)
                                         tgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATa                                                >  2:2951740/1‑93 (MQ=255)
                                         tgGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATa                                                >  1:1304688/1‑93 (MQ=255)
                                          ggCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAg                                               <  2:2537132/93‑1 (MQ=255)
                                           gCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGt                                              <  2:2262494/93‑1 (MQ=255)
                                            ccTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAgct                                             >  2:2901600/1‑91 (MQ=255)
                                            ccTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGtt                                             >  1:1712068/1‑93 (MQ=255)
                                             cTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTg                                            >  1:3903677/1‑93 (MQ=255)
                                              tACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGa                                           >  1:3439692/1‑93 (MQ=255)
                                               aCCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGaa                                          <  2:1331695/93‑1 (MQ=255)
                                                   aCATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAc                                      <  1:1421063/93‑1 (MQ=255)
                                                     atcatcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTg                                    >  2:3199120/1‑93 (MQ=255)
                                                      tcatcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGa                                   <  1:1574532/93‑1 (MQ=255)
                                                       catcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGAt                                  <  1:2259808/93‑1 (MQ=255)
                                                        atcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGAtt                                 >  1:1022560/1‑93 (MQ=255)
                                                        atcCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGAtt                                 <  1:631944/93‑1 (MQ=255)
                                                          ccAGGAGCCCAGAGACGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGc                               <  2:4212428/93‑1 (MQ=255)
                                                          ccAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGc                               <  1:890172/93‑1 (MQ=255)
                                                          ccAGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGc                               <  1:2229038/93‑1 (MQ=255)
                                                            aGGAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaa                             >  2:1008123/1‑93 (MQ=255)
                                                             ggAGCCCATAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTCATTGCaaa                            >  1:4609844/1‑93 (MQ=255)
                                                             ggAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaaa                            <  2:934981/93‑1 (MQ=255)
                                                             ggAGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCaaa                            >  2:4645904/1‑93 (MQ=255)
                                                               aGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACg                          <  1:1191893/93‑1 (MQ=255)
                                                               aGCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACg                          <  2:4128853/93‑1 (MQ=255)
                                                                gCCCAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACgg                         <  2:3893885/93‑1 (MQ=255)
                                                                  ccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAg                       >  2:4031528/1‑93 (MQ=255)
                                                                  ccAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAg                       <  2:3396872/93‑1 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      >  2:4594944/1‑93 (MQ=255)
                                                                   cAGAGAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGc                      >  2:2516246/1‑93 (MQ=255)
                                                                    agagAAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCCCTGATTGCAAACGGAGCa                     <  1:4001281/93‑1 (MQ=255)
                                                                     gagaAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCAc                    >  2:4456570/1‑93 (MQ=255)
                                                                     gagaAGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCAc                    >  2:205393/1‑93 (MQ=255)
                                                                        aaGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCgg                 <  2:307237/93‑1 (MQ=255)
                                                                         aGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGa                >  1:1149937/1‑93 (MQ=255)
                                                                         aGTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGa                <  2:2078792/93‑1 (MQ=255)
                                                                          gTCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGaa               >  2:3026443/1‑93 (MQ=255)
                                                                           tCAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAg              >  1:1611192/1‑93 (MQ=255)
                                                                            cAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgc             >  2:3549232/1‑93 (MQ=255)
                                                                            cAGGCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgc             <  2:193498/93‑1 (MQ=255)
                                                                              ggCACGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAgccg           >  1:1260316/1‑93 (MQ=255)
                                                                                  cGCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCAt       <  1:1949487/93‑1 (MQ=255)
                                                                                   gCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCAtt      >  1:4430544/1‑93 (MQ=255)
                                                                                   gCAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCAtt      >  2:379008/1‑93 (MQ=255)
                                                                                    cAAAGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTa     >  1:2017582/1‑93 (MQ=255)
                                                                                     aaaGGTGCTAAATGCAGAAACCAGCGGCGCAGCCCAGGTACCGAGGAGAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTAc    <  1:3531117/93‑1 (MQ=255)
                                                                                      aaGGTGGTGAATGCAGAAACCAGCGGCGCATCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACa   <  1:4277613/93‑1 (MQ=255)
                                                                                      aaGGTGGGGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAAACGAGCACCGGAAGCCGCCATTACa   >  2:2589998/1‑93 (MQ=255)
                                                                                       aGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCAGTGATTGGAAACGGAGCACCGGCAGGCGGCATTACag  >  2:389448/1‑93 (MQ=255)
                                                                                       aGGTGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACag  >  2:4043289/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAG  >  NC_000913/4356685‑4356864

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