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‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665090‑1665282
                                                                                                   |                                                                                                
tCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                   >  1:1951244/1‑99 (MQ=255)
 gaaacagcgGTTTAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑t                                                                                              >  1:1703128/7‑100 (MQ=255)
  cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt                                                                                             >  2:2620853/1‑100 (MQ=255)
   tGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑ttt                                                                                            >  2:3574100/1‑100 (MQ=255)
    gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt                                                                                           >  2:1428853/1‑100 (MQ=255)
      cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttt                                                                                         <  1:3653776/100‑1 (MQ=255)
      cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttt                                                                                         <  2:1512681/100‑1 (MQ=255)
       gctgTTGAATCGCGGTTTCGATAGTGTCACCGTAATGTTGCCGGATCGTGTCGATTGCGCGCTATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttt                                                                                         <  1:1008507/96‑1 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑ttttttt                                                                                        >  2:3415831/1‑99 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑ttttttt                                                                                        >  2:3502197/1‑99 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttt                                                                                         <  1:3733107/99‑1 (MQ=255)
        cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑ttttttta                                                                                       <  2:345740/100‑3 (MQ=255)
          gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttttatc                                                                                     >  2:2443043/1‑96 (MQ=255)
            tAAACGGCTGTTTCGATAGGATCAGCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttttttatc                                                                                     <  1:2714855/98‑5 (MQ=255)
            tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTAtctc                                                                                    >  1:429986/1‑100 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCtt                                                                                  <  2:3722783/100‑1 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCtt                                                                                  <  1:1802639/100‑1 (MQ=255)
              aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCtt                                                                                  >  1:8328/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTaa                                                                                >  2:3599404/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTaa                                                                                >  2:3138245/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTaa                                                                                >  1:2313700/1‑100 (MQ=255)
                 cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTaaa                                                                               <  2:283451/100‑1 (MQ=255)
                 cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTaaa                                                                               >  1:2287931/1‑100 (MQ=255)
                  gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAg                                                                              <  2:2434246/100‑1 (MQ=255)
                   cTGGTTCGATAGGATCACCGTAACATTTTCGGATTTTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAg                                                                              <  1:1208057/99‑1 (MQ=255)
                      tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAgcacg                                                                          <  2:1133078/100‑1 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAgcacgcac                                                                       >  2:3938313/1‑100 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAgcacgcac                                                                       <  1:2708607/100‑1 (MQ=255)
                          ggtAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGTAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:2889926/3‑100 (MQ=255)
                          gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:1050954/1‑100 (MQ=255)
                          gATAGGATCACCGTAATATTGCCGGATGGTGACGAGTACGCGCAATCCGGCAATAAGAGGTTACAGTGTCAGGT‑T‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  1:844805/1‑100 (MQ=255)
                           aTAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTg                                                                     <  1:3471403/100‑1 (MQ=255)
                            tAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGc                                                                    >  2:2259184/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCCATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCAAGCACTGCtt                                                                  >  1:1773914/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATTTCTTAAAGAACGCACtac                                                                    >  1:431593/1‑96 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  2:2489232/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  1:2797808/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  <  2:3806278/100‑1 (MQ=255)
                               gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAATCACGCACTGCttt                                                                 <  2:2904236/100‑1 (MQ=255)
                                 tCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                               <  2:3154422/100‑1 (MQ=255)
                                   aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:3990957/100‑1 (MQ=255)
                                   aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:3595723/100‑1 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            >  2:2173936/1‑100 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            <  2:3514041/100‑1 (MQ=255)
                                      gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                          >  1:2752718/1‑100 (MQ=255)
                                       tAATATTGTAGGATGGTGACGATTACGCGCAAGCCGTCAAGAATAGGTTACAGTGTCTCGT‑T‑TTTTTTATCTCTGAAAGCACGCACTGCTTTTGCGGCTg                                                         <  1:1477972/100‑1 (MQ=255)
                                        aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  1:2166925/100‑1 (MQ=255)
                                            ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                     <  1:4207781/100‑1 (MQ=255)
                                            ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TGTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                    <  1:1290116/100‑1 (MQ=255)
                                             tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                   <  2:3271109/100‑1 (MQ=255)
                                              gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                  >  1:2823809/1‑100 (MQ=255)
                                               ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 >  1:4089258/1‑100 (MQ=255)
                                               ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 >  2:3062686/1‑100 (MQ=255)
                                                cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                 <  1:931537/99‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                               <  2:303231/100‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCAATCCGGCAAGAATAGGTTACAGTGTCACGT‑T‑TTTTTTATGTCTTAAAGCACGCACGGCTTTTggtgctagc                                                       >  1:2447972/1‑84 (MQ=255)
                                                  gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                              <  1:3000537/100‑1 (MQ=255)
                                                    tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                              <  2:3364622/98‑1 (MQ=255)
                                                       tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                         >  2:543901/1‑100 (MQ=255)
                                                        gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                        <  1:3574087/100‑1 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                      <  1:3256350/100‑1 (MQ=255)
                                                            aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTTTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                             >  1:2408759/1‑91 (MQ=255)
                                                            aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTACGTTTTAAAGCACCCCCTGCTTTTGCGGTTGGCCTCTTTTGCGGCAAAATAg                                    >  2:2987836/1‑100 (MQ=255)
                                                             ttACGCGGAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCTCTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                   <  1:2289146/100‑1 (MQ=255)
                                                              tACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc                                  <  1:4177109/100‑1 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 <  1:1645582/100‑1 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 >  2:1407568/1‑100 (MQ=255)
                                                               aCGCGCAATCCGGCAACAATAGGTTACAGTGACACGT‑T‑TTTTTTATCTCTTAAAGCAAGCAATGCTTGTGTGGATAGCCTCTTTTGCCCCAAAAGAGt                                   >  1:108840/1‑98 (MQ=255)
                                                                cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  2:1903707/1‑100 (MQ=255)
                                                                   gcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                             <  1:1714644/100‑1 (MQ=255)
                                                                    cAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                            <  1:1575408/100‑1 (MQ=255)
                                                                       tCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                         <  2:2910498/100‑1 (MQ=255)
                                                                        ccGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                         <  2:1887767/99‑1 (MQ=255)
                                                                         cGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                       >  2:3251750/1‑100 (MQ=255)
                                                                          ggCAATAATAGGTTACGGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                      <  1:128651/100‑1 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTTTCTCTTAAAGCACGCACTGCTTGTGCGGCTGGCCTCGTTTTCCGCGAAATAATCGCCCgtgt                          >  1:1603025/1‑96 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCTCCCGTGTTTCa                      >  1:1535798/1‑100 (MQ=255)
                                                                          ggCAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTGTTGCCGCAAAATAGTCGCCCGTGTTTc                       <  1:2436796/99‑1 (MQ=255)
                                                                            cAATAATAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                    >  2:280761/1‑100 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGT‑T‑TTTTTTTTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                    <  1:1280359/99‑1 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                   <  2:2938713/100‑1 (MQ=255)
                                                                              ataataGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                  <  2:3376095/100‑1 (MQ=255)
                                                                               taataGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                 <  2:2735439/100‑1 (MQ=255)
                                                                                aataGGTTACAGTGTCACGT‑T‑TTTTTTTTCTTTTAAAGCACGTACTGCTTTTGCGGCTGGCCTCTTTTTCCGCAAAATAGTCGCCCTTGTTTCAtt                    >  1:3841817/1‑96 (MQ=255)
                                                                                aataGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                <  1:3958157/100‑1 (MQ=255)
                                                                                 ataGGTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                >  1:4214411/1‑99 (MQ=255)
                                                                                 ataGGTTAAAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa               <  2:1116396/100‑1 (MQ=255)
                                                                                   aGGTCAAAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGTCTGGCCTTTTTTGCCGCAAAATAGTCGCCCGTGTTTCCTTGCCCAtt             >  1:3532360/1‑100 (MQ=255)
                                                                                    ggTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt            >  1:2366800/1‑100 (MQ=255)
                                                                                    ggTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt            >  2:4121708/1‑100 (MQ=255)
                                                                                     gTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  1:1739503/100‑1 (MQ=255)
                                                                                     gTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           >  2:129609/1‑100 (MQ=255)
                                                                                     gTTACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCCCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  2:2428042/100‑1 (MQ=255)
                                                                                       tACAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTcgg         <  2:2902196/100‑3 (MQ=255)
                                                                                         cAGTGTCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt       >  2:4266193/1‑100 (MQ=255)
                                                                                             gtCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  1:2057496/100‑1 (MQ=255)
                                                                                             gtCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:2980360/100‑1 (MQ=255)
                                                                                             gtCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:1634479/100‑1 (MQ=255)
                                                                                             gtCACGT‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:2147642/100‑1 (MQ=255)
                                                                                              tCACGT‑T‑TTTTTTATCTCTTAAAGCACGAACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  1:122857/99‑1 (MQ=255)
                                                                                                   |                                                                                                
TCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  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: ‑