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.T96.5% 282.4 / 2.8 86intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base T (45/38);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                    |                                                                                                
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                   >  2:1640294/1‑100 (MQ=255)
  ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATAGTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑t                                                                                              <  2:2981599/100‑1 (MQ=255)
    tGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑ttt                                                                                            >  1:1436848/1‑100 (MQ=255)
     gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt                                                                                           >  1:647786/1‑100 (MQ=255)
      tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑t                                                                                         <  1:1818905/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tt                                                                                        <  1:3780900/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tt                                                                                        <  2:2197828/100‑1 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttt                                                                                       <  2:3163881/100‑2 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttt                                                                                       <  1:2156622/100‑2 (MQ=255)
         cgGTTAAACCTCTGTTTCGATAGTATCTCCTTAATATTGCTGGATGGTGACTATTTCGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑t                                                                                         <  2:2847550/97‑1 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCACTCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttta                                                                                      <  2:2232766/100‑3 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttta                                                                                      >  1:4072289/1‑98 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttt                                                                                       <  2:111324/99‑2 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tttat                                                                                     >  1:2418946/1‑97 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tttat                                                                                     >  2:631598/1‑97 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tttatc                                                                                    <  1:908597/100‑5 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑tttatct                                                                                   <  1:3296467/100‑6 (MQ=255)
             tgaaCCGCGGTTGCGATAGGATCACCGTAATATTGCAGGATGGTGGCGATTACGCGCAATCCGACAATAATAGGTTACAGTGTCACG‑‑t‑tttt‑ttta                                                                                      <  2:2416198/94‑3 (MQ=255)
              aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATctct                                                                                  >  1:2192332/1‑100 (MQ=255)
              aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATctct                                                                                  >  1:576153/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCtt                                                                                 >  1:843338/1‑100 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGACGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTa                                                                                >  2:3979293/1‑100 (MQ=255)
                  cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTAGGCGCATTCCGGAAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTaaa                                                                              <  2:1137422/100‑1 (MQ=255)
                    cTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAgc                                                                            <  1:3389051/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGAGTCACG‑‑TATTTT‑TTATCTCTTAAAgta                                                                           >  2:2405381/1‑98 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAgcacgc                                                                        >  2:3216258/1‑100 (MQ=255)
                         tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTTTCTCTTAAAgc                                                                            >  2:3372023/1‑95 (MQ=255)
                          cGATAGGATCACCGTAATTTTGCCGGATGTTGACGATTACTCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAgcacgcac                                                                      <  1:3597128/100‑1 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAgcacgcac                                                                      <  2:1976929/100‑1 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TTTT‑TTATCTCTTAAAGCACGCACt                                                                     >  2:3800176/1‑100 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAgcacgcac                                                                      <  2:3476900/99‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACt                                                                     >  1:421438/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCtt                                                                 >  2:1731558/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCtt                                                                 >  2:2693030/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:4129188/1‑100 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCtttt                                                               <  2:3784114/100‑1 (MQ=255)
                                  tCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTg                                                              <  1:2912863/100‑1 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCtgtggc                                                             <  2:3287143/100‑6 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            >  2:2815507/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           <  1:2302489/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  1:2474607/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  2:388009/99‑1 (MQ=255)
                                      cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           <  1:3492732/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                         >  1:377278/1‑100 (MQ=255)
                                       gTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  2:4041797/99‑1 (MQ=255)
                                       gTAAGATTGCCGGATGGTGACGATTACGTGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  2:3697821/99‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  1:3982713/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  1:984488/100‑1 (MQ=255)
                                            atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                    >  1:3543182/1‑100 (MQ=255)
                                            atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                    >  1:979731/1‑100 (MQ=255)
                                             ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                   <  1:1283124/100‑1 (MQ=255)
                                              tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                  >  1:2697948/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTATGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  1:2274065/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATTTTTTAAAGTACACACTGCTTTTGCGGCTg                                                        >  2:1340214/1‑93 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  1:1399954/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TTTTGTTATCTCTTAAAGCAAGCACTGCTTTTGCGGATGGCCTCtt                                                 >  1:102028/1‑100 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                               >  2:3585667/1‑100 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                               >  1:3319313/1‑100 (MQ=255)
                                                  ggATGGTGACGATTATGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                              <  1:4110684/100‑1 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             >  2:4167818/1‑100 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  2:2091885/100‑1 (MQ=255)
                                                    aTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            <  1:3352200/100‑1 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  1:646229/1‑100 (MQ=255)
                                                        tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTTTTTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                        >  2:1623775/1‑100 (MQ=255)
                                                          aCGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa                                      <  1:3449054/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  1:928523/1‑100 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  1:2613952/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  1:1306382/1‑100 (MQ=255)
                                                            gATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  2:2958680/99‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATTTCTTAATTCACGCACTGCTTTTGCGGCTGGCCCCTTTTGCCGCAACATCGTc                                 >  1:623479/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               <  2:866492/100‑1 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              >  2:374407/1‑100 (MQ=255)
                                                                   cgcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCcccc                             >  2:3346254/1‑100 (MQ=255)
                                                                    gcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                            <  2:3467885/100‑1 (MQ=255)
                                                                     cAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGCCTGGCTTCTTTTGCCGCAAACTAGTCGccc                             >  2:4039305/1‑98 (MQ=255)
                                                                        tCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        >  1:3585095/1‑100 (MQ=255)
                                                                         ccGGCAATAATAGGTTACAGTGTCACGTTT‑TTTT‑TTTTTTTTTAAAGCACGCACTGCTTTTGCGGCTTTCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        >  2:2630655/1‑100 (MQ=255)
                                                                           ggCAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                     <  2:844097/100‑1 (MQ=255)
                                                                             cAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   <  2:699683/100‑1 (MQ=255)
                                                                             cAATAATAGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   >  1:2713628/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 <  2:2896487/100‑1 (MQ=255)
                                                                                taataGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  1:548183/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc               <  1:3150305/100‑1 (MQ=255)
                                                                                  ataGGTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  2:2754179/98‑1 (MQ=255)
                                                                                     ggTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCCTTGCCCAttt           >  2:2306923/1‑100 (MQ=255)
                                                                                     ggTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt           >  2:4137417/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  1:3129367/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  1:3633257/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  1:601611/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGT‑T‑TTTT‑TTTTTTTTTTAAGCACGCACTGCTTTTACGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        >  2:1853264/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        >  2:1992338/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        <  1:1731568/100‑1 (MQ=255)
                                                                                        tACAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        <  1:218660/100‑1 (MQ=255)
                                                                                         aCAGTGTCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc       <  1:2451544/100‑1 (MQ=255)
                                                                                            gtgtCGCGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc     <  1:1823928/99‑1 (MQ=255)
                                                                                            gtgtCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt      <  1:3446905/98‑1 (MQ=255)
                                                                                            gtgtCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa    >  1:2388273/1‑100 (MQ=255)
                                                                                              gtCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:925374/100‑1 (MQ=255)
                                                                                               tCACGT‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  2:1620435/99‑1 (MQ=255)
                                                                                               tCACGT‑T‑TTTT‑TTATCTCTTAAAGCACCCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCCAACTAGTCGCCCGTCTTTCCTTGCCCCTTTCtccccctg  >  2:3650322/1‑92 (MQ=255)
                                                                                                    |                                                                                                
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TTTT‑TTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281

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: ‑