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‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665089‑1665281
                                                                                                    |                                                                                                
ttCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                                                                   >  2:1640294/1‑100 (MQ=255)
  ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATAGTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tt                                                                                              <  2:2981599/100‑1 (MQ=255)
    tGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttt                                                                                            >  1:1436848/1‑100 (MQ=255)
     gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttt                                                                                           >  1:647786/1‑100 (MQ=255)
      tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttt                                                                                         <  1:1818905/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                        <  1:3780900/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                        <  2:2197828/100‑1 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttt                                                                                       <  2:3163881/100‑2 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttt                                                                                       <  1:2156622/100‑2 (MQ=255)
         cgGTTAAACCTCTGTTTCGATAGTATCTCCTTAATATTGCTGGATGGTGACTATTTCGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttt                                                                                         <  2:2847550/97‑1 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCACTCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttta                                                                                      <  2:2232766/100‑3 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttta                                                                                      >  1:4072289/1‑98 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tttttttt                                                                                       <  2:111324/99‑2 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttat                                                                                     >  1:2418946/1‑97 (MQ=255)
          ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttat                                                                                     >  2:631598/1‑97 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                    <  1:908597/100‑5 (MQ=255)
            ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatct                                                                                   <  1:3296467/100‑6 (MQ=255)
             tgaaCCGCGGTTGCGATAGGATCACCGTAATATTGCAGGATGGTGGCGATTACGCGCAATCCGACAATAATAGGTTACAGTGTCACG‑‑tttttttta                                                                                      <  2:2416198/94‑3 (MQ=255)
              aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATctct                                                                                  >  1:2192332/1‑100 (MQ=255)
              aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATctct                                                                                  >  1:576153/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCtt                                                                                 >  1:843338/1‑100 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGACGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTa                                                                                >  2:3979293/1‑100 (MQ=255)
                  cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTAGGCGCATTCCGGAAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaa                                                                              <  2:1137422/100‑1 (MQ=255)
                    cTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgc                                                                            <  1:3389051/100‑1 (MQ=255)
                     tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGAGTCACG‑‑TATTTTTTATCTCTTAAAgta                                                                           >  2:2405381/1‑98 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgc                                                                        >  2:3216258/1‑100 (MQ=255)
                         tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTCTCTTAAAgc                                                                            >  2:3372023/1‑95 (MQ=255)
                          cGATAGGATCACCGTAATTTTGCCGGATGTTGACGATTACTCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  1:3597128/100‑1 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  2:1976929/100‑1 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACt                                                                     >  2:3800176/1‑100 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  2:3476900/99‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACt                                                                     >  1:421438/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  2:1731558/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  2:2693030/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:4129188/1‑100 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtttt                                                               <  2:3784114/100‑1 (MQ=255)
                                  tCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                              <  1:2912863/100‑1 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtgtggc                                                             <  2:3287143/100‑6 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            >  2:2815507/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           <  1:2302489/100‑1 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  1:2474607/1‑100 (MQ=255)
                                     ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  2:388009/99‑1 (MQ=255)
                                      cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           <  1:3492732/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                         >  1:377278/1‑100 (MQ=255)
                                       gTAATATTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  2:4041797/99‑1 (MQ=255)
                                       gTAAGATTGCCGGATGGTGACGATTACGTGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                          <  2:3697821/99‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        <  1:3982713/100‑1 (MQ=255)
                                           tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  1:984488/100‑1 (MQ=255)
                                            atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                    >  1:3543182/1‑100 (MQ=255)
                                            atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                    >  1:979731/1‑100 (MQ=255)
                                             ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                   <  1:1283124/100‑1 (MQ=255)
                                              tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                  >  1:2697948/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTATGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  1:2274065/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATTTTTTAAAGTACACACTGCTTTTGCGGCTg                                                        >  2:1340214/1‑93 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 >  1:1399954/1‑100 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTGTTATCTCTTAAAGCAAGCACTGCTTTTGCGGATGGCCTCtt                                                 >  1:102028/1‑100 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                               >  2:3585667/1‑100 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                               >  1:3319313/1‑100 (MQ=255)
                                                  ggATGGTGACGATTATGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                              <  1:4110684/100‑1 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             >  2:4167818/1‑100 (MQ=255)
                                                   gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  2:2091885/100‑1 (MQ=255)
                                                    aTGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGcc                                            <  1:3352200/100‑1 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  1:646229/1‑100 (MQ=255)
                                                        tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTTTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                        >  2:1623775/1‑100 (MQ=255)
                                                          aCGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa                                      <  1:3449054/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  1:928523/1‑100 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  1:2613952/100‑1 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  1:1306382/1‑100 (MQ=255)
                                                            gATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  2:2958680/99‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATTTCTTAATTCACGCACTGCTTTTGCGGCTGGCCCCTTTTGCCGCAACATCGTc                                 >  1:623479/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               <  2:866492/100‑1 (MQ=255)
                                                                  gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                              >  2:374407/1‑100 (MQ=255)
                                                                   cgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCcccc                             >  2:3346254/1‑100 (MQ=255)
                                                                    gcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                            <  2:3467885/100‑1 (MQ=255)
                                                                     cAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGCCTGGCTTCTTTTGCCGCAAACTAGTCGccc                             >  2:4039305/1‑98 (MQ=255)
                                                                        tCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        >  1:3585095/1‑100 (MQ=255)
                                                                         ccGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTTTTTTTTAAAGCACGCACTGCTTTTGCGGCTTTCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                        >  2:2630655/1‑100 (MQ=255)
                                                                           ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                     <  2:844097/100‑1 (MQ=255)
                                                                             cAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   <  2:699683/100‑1 (MQ=255)
                                                                             cAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                   >  1:2713628/1‑100 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 <  2:2896487/100‑1 (MQ=255)
                                                                                taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  1:548183/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc               <  1:3150305/100‑1 (MQ=255)
                                                                                  ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                <  2:2754179/98‑1 (MQ=255)
                                                                                     ggTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCCTTGCCCAttt           >  2:2306923/1‑100 (MQ=255)
                                                                                     ggTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt           >  2:4137417/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  1:3129367/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  1:3633257/1‑100 (MQ=255)
                                                                                      gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          >  1:601611/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGTTTTTTTTTTTTTTTTAAGCACGCACTGCTTTTACGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        >  2:1853264/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        >  2:1992338/1‑100 (MQ=255)
                                                                                        tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        <  1:1731568/100‑1 (MQ=255)
                                                                                        tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        <  1:218660/100‑1 (MQ=255)
                                                                                         aCAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc       <  1:2451544/100‑1 (MQ=255)
                                                                                            gtgtCGCGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc     <  1:1823928/99‑1 (MQ=255)
                                                                                            gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt      <  1:3446905/98‑1 (MQ=255)
                                                                                            gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa    >  1:2388273/1‑100 (MQ=255)
                                                                                              gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  1:925374/100‑1 (MQ=255)
                                                                                               tCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  <  2:1620435/99‑1 (MQ=255)
                                                                                               tCACGTTTTTTTTATCTCTTAAAGCACCCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCCAACTAGTCGCCCGTCTTTCCTTGCCCCTTTCtccccctg  >  2:3650322/1‑92 (MQ=255)
                                                                                                    |                                                                                                
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  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: