Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,665,195:1 (T)7→8 100% intergenic (+75/‑120) dmsD → / → clcB twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,1881.T98.1% 356.9 / ‑0.4 107intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (48/57);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 127 total aligned reads displayed.

TCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665090‑1665282
                                                                                                   |                                                                                                
tCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                   >  1:1951244/1‑99 (MQ=255)
 gaaacagcgGTTTAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tt                                                                                              >  1:1703128/7‑100 (MQ=255)
  cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttt                                                                                             >  2:2620853/1‑100 (MQ=255)
   tGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttt                                                                                            >  2:3574100/1‑100 (MQ=255)
    gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttt                                                                                           >  2:1428853/1‑100 (MQ=255)
      cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                         <  1:3653776/100‑1 (MQ=255)
      cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                         <  2:1512681/100‑1 (MQ=255)
       gctgTTGAATCGCGGTTTCGATAGTGTCACCGTAATGTTGCCGGATCGTGTCGATTGCGCGCTATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                         <  1:1008507/96‑1 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttt                                                                                        >  2:3415831/1‑99 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttt                                                                                        >  2:3502197/1‑99 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                         <  1:3733107/99‑1 (MQ=255)
        cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttta                                                                                       <  2:345740/100‑3 (MQ=255)
          gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                     >  2:2443043/1‑96 (MQ=255)
            tAAACGGCTGTTTCGATAGGATCAGCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                     <  1:2714855/98‑5 (MQ=255)
            tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTAtctc                                                                                    >  1:429986/1‑100 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                  <  2:3722783/100‑1 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                  <  1:1802639/100‑1 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                  >  1:8328/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                                >  2:3599404/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                                >  2:3138245/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                                >  1:2313700/1‑100 (MQ=255)
                 cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaa                                                                               <  2:283451/100‑1 (MQ=255)
                 cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaa                                                                               >  1:2287931/1‑100 (MQ=255)
                  gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                              <  2:2434246/100‑1 (MQ=255)
                   cTGGTTCGATAGGATCACCGTAACATTTTCGGATTTTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                              <  1:1208057/99‑1 (MQ=255)
                      tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacg                                                                          <  2:1133078/100‑1 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                       >  2:3938313/1‑100 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                       <  1:2708607/100‑1 (MQ=255)
                          ggtAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGTAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:2889926/3‑100 (MQ=255)
                          gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:1050954/1‑100 (MQ=255)
                          gATAGGATCACCGTAATATTGCCGGATGGTGACGAGTACGCGCAATCCGGCAATAAGAGGTTACAGTGTCAGGTTTTTTTTATCTCTTAAAGCACGCACt                                                                      >  1:844805/1‑100 (MQ=255)
                           aTAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTg                                                                     <  1:3471403/100‑1 (MQ=255)
                            tAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGc                                                                    >  2:2259184/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCCATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCAAGCACTGCtt                                                                  >  1:1773914/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATTTCTTAAAGAACGCACtac                                                                    >  1:431593/1‑96 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  2:2489232/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  1:2797808/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  <  2:3806278/100‑1 (MQ=255)
                               gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAATCACGCACTGCttt                                                                 <  2:2904236/100‑1 (MQ=255)
                                 tCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                               <  2:3154422/100‑1 (MQ=255)
                                   aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:3990957/100‑1 (MQ=255)
                                   aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:3595723/100‑1 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            >  2:2173936/1‑100 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            <  2:3514041/100‑1 (MQ=255)
                                      gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                          >  1:2752718/1‑100 (MQ=255)
                                       tAATATTGTAGGATGGTGACGATTACGCGCAAGCCGTCAAGAATAGGTTACAGTGTCTCGTTTTTTTTATCTCTGAAAGCACGCACTGCTTTTGCGGCTg                                                         <  1:1477972/100‑1 (MQ=255)
                                        aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  1:2166925/100‑1 (MQ=255)
                                            ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                     <  1:4207781/100‑1 (MQ=255)
                                            ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TGTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                    <  1:1290116/100‑1 (MQ=255)
                                             tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                   <  2:3271109/100‑1 (MQ=255)
                                              gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                  >  1:2823809/1‑100 (MQ=255)
                                               ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 >  1:4089258/1‑100 (MQ=255)
                                               ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 >  2:3062686/1‑100 (MQ=255)
                                                cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 <  1:931537/99‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                               <  2:303231/100‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCAATCCGGCAAGAATAGGTTACAGTGTCACGTTTTTTTTATGTCTTAAAGCACGCACGGCTTTTggtgctagc                                                       >  1:2447972/1‑84 (MQ=255)
                                                  gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                              <  1:3000537/100‑1 (MQ=255)
                                                    tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                              <  2:3364622/98‑1 (MQ=255)
                                                       tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                         >  2:543901/1‑100 (MQ=255)
                                                        gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                        <  1:3574087/100‑1 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                      <  1:3256350/100‑1 (MQ=255)
                                                            aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                             >  1:2408759/1‑91 (MQ=255)
                                                            aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTACGTTTTAAAGCACCCCCTGCTTTTGCGGTTGGCCTCTTTTGCGGCAAAATAg                                    >  2:2987836/1‑100 (MQ=255)
                                                             ttACGCGGAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCTCTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                   <  1:2289146/100‑1 (MQ=255)
                                                              tACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc                                  <  1:4177109/100‑1 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 <  1:1645582/100‑1 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 >  2:1407568/1‑100 (MQ=255)
                                                               aCGCGCAATCCGGCAACAATAGGTTACAGTGACACGTTTTTTTTATCTCTTAAAGCAAGCAATGCTTGTGTGGATAGCCTCTTTTGCCCCAAAAGAGt                                   >  1:108840/1‑98 (MQ=255)
                                                                cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  2:1903707/1‑100 (MQ=255)
                                                                   gcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                             <  1:1714644/100‑1 (MQ=255)
                                                                    cAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                            <  1:1575408/100‑1 (MQ=255)
                                                                       tCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                         <  2:2910498/100‑1 (MQ=255)
                                                                        ccGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                         <  2:1887767/99‑1 (MQ=255)
                                                                         cGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                       >  2:3251750/1‑100 (MQ=255)
                                                                          ggCAATAATAGGTTACGGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                      <  1:128651/100‑1 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGTTTTTTTTTTCTCTTAAAGCACGCACTGCTTGTGCGGCTGGCCTCGTTTTCCGCGAAATAATCGCCCgtgt                          >  1:1603025/1‑96 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCTCCCGTGTTTCa                      >  1:1535798/1‑100 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTGTTGCCGCAAAATAGTCGCCCGTGTTTc                       <  1:2436796/99‑1 (MQ=255)
                                                                            cAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                    >  2:280761/1‑100 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGTTTTTTTTTTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                    <  1:1280359/99‑1 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                   <  2:2938713/100‑1 (MQ=255)
                                                                              ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                  <  2:3376095/100‑1 (MQ=255)
                                                                               taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                 <  2:2735439/100‑1 (MQ=255)
                                                                                aataGGTTACAGTGTCACGTTTTTTTTTTCTTTTAAAGCACGTACTGCTTTTGCGGCTGGCCTCTTTTTCCGCAAAATAGTCGCCCTTGTTTCAtt                    >  1:3841817/1‑96 (MQ=255)
                                                                                aataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                <  1:3958157/100‑1 (MQ=255)
                                                                                 ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                >  1:4214411/1‑99 (MQ=255)
                                                                                 ataGGTTAAAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa               <  2:1116396/100‑1 (MQ=255)
                                                                                   aGGTCAAAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGTCTGGCCTTTTTTGCCGCAAAATAGTCGCCCGTGTTTCCTTGCCCAtt             >  1:3532360/1‑100 (MQ=255)
                                                                                    ggTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt            >  1:2366800/1‑100 (MQ=255)
                                                                                    ggTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt            >  2:4121708/1‑100 (MQ=255)
                                                                                     gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  1:1739503/100‑1 (MQ=255)
                                                                                     gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           >  2:129609/1‑100 (MQ=255)
                                                                                     gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCCCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  2:2428042/100‑1 (MQ=255)
                                                                                       tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTcgg         <  2:2902196/100‑3 (MQ=255)
                                                                                         cAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt       >  2:4266193/1‑100 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  1:2057496/100‑1 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:2980360/100‑1 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:1634479/100‑1 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:2147642/100‑1 (MQ=255)
                                                                                              tCACGTTTTTTTTATCTCTTAAAGCACGAACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  1:122857/99‑1 (MQ=255)
                                                                                                   |                                                                                                
TCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665090‑1665282

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