Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,195 (T)7→8 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.4% 358.2 / 6.0 111intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base T (43/63);  total (47/64)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665089‑1665282
                                                                                                    |                                                                                                
 tCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t                                                                                                 >  1:2687918/1‑100 (MQ=255)
  ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑t                                                                                               >  2:1309745/1‑100 (MQ=255)
      tCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑ttt                                                                                          <  2:1548334/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttt                                                                                         >  2:1493739/1‑100 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttt                                                                                         <  1:310345/100‑1 (MQ=255)
       cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttt                                                                                         >  2:19041/1‑100 (MQ=255)
        ggtgTTAGACCGCTGTGTCGAGAGGATTAGCGTAATAGTGACGAATGGTGATGATTATGCTCGTTCCGGCAATAATAGGCAACAGGGTCACG‑t‑tt‑tttt                                                                                         <  1:1419025/96‑1 (MQ=255)
        gcgcTTAAACCGCTGTTTCGATAGTATCACCGGAATATTGCCGGATGGTGACTATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttt                                                                                         <  2:1156820/99‑1 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑ttttt                                                                                        <  2:2509978/100‑2 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑ttttt                                                                                        >  1:290609/1‑99 (MQ=255)
        gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttt                                                                                         <  1:2003483/99‑1 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑ttttta                                                                                       >  2:3072093/1‑98 (MQ=255)
         cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑ttttta                                                                                       >  2:1171734/1‑98 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttttatc                                                                                     <  1:1045461/100‑5 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttttatc                                                                                     <  1:1831011/100‑5 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttttatc                                                                                     >  1:733371/1‑96 (MQ=255)
           gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt‑tttttatc                                                                                     >  2:2865325/1‑96 (MQ=255)
             tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTAtctc                                                                                    <  1:566815/100‑1 (MQ=255)
             tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTAtctc                                                                                    <  1:2960432/100‑1 (MQ=255)
             tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTAtctc                                                                                    >  1:2363554/1‑100 (MQ=255)
              aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATctct                                                                                   >  1:1256399/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCtt                                                                                  <  2:1602486/100‑1 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCtt                                                                                  >  1:3245321/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCtt                                                                                  >  1:3321408/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCtt                                                                                  >  1:3323925/1‑100 (MQ=255)
               aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCtt                                                                                  <  1:458657/100‑1 (MQ=255)
                aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTa                                                                                 >  1:737525/1‑100 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTaa                                                                                >  2:1558213/1‑100 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTaa                                                                                <  2:1297277/100‑1 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTaa                                                                                <  2:1063225/100‑1 (MQ=255)
                 ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑TATT‑TTTTATCTCTTaa                                                                                >  2:3418222/1‑100 (MQ=255)
                   gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAg                                                                              <  1:1433350/100‑1 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGTCGGAGGGTGCCGATTCCGCGCAATCCGGCAATAATAGGTTAGAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgc                                                                             <  1:1893202/98‑1 (MQ=255)
                      gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcac                                                                           >  2:3107140/1‑100 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgc                                                                         <  1:2018820/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgc                                                                         <  2:191457/100‑1 (MQ=255)
                        ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgc                                                                         <  1:1644776/100‑1 (MQ=255)
                         tgtgtAGGATCACCGTAATGTTGCCTGATGGTGACGATTACGCGCAATCGGGCAATAATAGGTTACAGGGTCACG‑T‑TTGTTTTATCTCTTAAAgcacgc                                                                         <  1:3434447/95‑1 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgcac                                                                       >  2:1466191/1‑100 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgcac                                                                       >  2:2133318/1‑100 (MQ=255)
                          cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAgcacgcac                                                                       <  1:1746905/100‑1 (MQ=255)
                           gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACt                                                                      <  1:3383944/100‑1 (MQ=255)
                             tAGGATCACCGTAATATTCCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGc                                                                    <  2:2769656/100‑1 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTTTTTCTTAAAGCACGCACTGCt                                                                   >  1:3355836/1‑99 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  2:3115900/1‑100 (MQ=255)
                               ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCttt                                                                 >  2:3282212/1‑100 (MQ=255)
                                gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCttt                                                                 >  2:3327398/1‑100 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCtttt                                                                <  1:1722013/100‑1 (MQ=255)
                                 aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCtaga                                                                >  1:18222/1‑97 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                              <  1:1589728/100‑1 (MQ=255)
                                   cACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                              <  2:1913064/100‑1 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:2030021/100‑1 (MQ=255)
                                    aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:205434/100‑1 (MQ=255)
                                      cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                           <  2:2140429/100‑1 (MQ=255)
                                       gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTg                                                               <  1:480517/95‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                         <  1:642395/100‑1 (MQ=255)
                                        tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                         <  2:2914554/100‑1 (MQ=255)
                                         aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  2:526703/100‑1 (MQ=255)
                                          atatTGCCGTATGGAGACGATTACTCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                       <  2:2236358/100‑1 (MQ=255)
                                          atatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                       <  2:2682584/100‑1 (MQ=255)
                                           tatTGCCGGATGGTAACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                      <  2:1819124/100‑1 (MQ=255)
                                               gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                  <  1:842602/100‑1 (MQ=255)
                                                 cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtttt                                                >  2:2914954/1‑100 (MQ=255)
                                                  ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                               <  1:2268807/100‑1 (MQ=255)
                                                     tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                            <  2:2182144/100‑1 (MQ=255)
                                                       gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                          >  1:282215/1‑100 (MQ=255)
                                                        tgatgaTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTTTTTTTTTTAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                         >  1:604428/1‑100 (MQ=255)
                                                         gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                        >  2:746191/1‑100 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                      >  1:2964483/1‑100 (MQ=255)
                                                           cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                      >  2:2152012/1‑100 (MQ=255)
                                                              ttACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaaa                                       <  1:2111851/96‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 <  2:821813/100‑1 (MQ=255)
                                                                aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 <  1:1226915/100‑1 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTTTTAAATTACGCAATGCTTCTGCAGCTACCCTCTTTTGATGCAAAATTGTCg                                 >  1:329085/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                 >  1:1419826/1‑99 (MQ=255)
                                                                 cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  2:2585560/1‑100 (MQ=255)
                                                                  tcaccaTCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                   >  1:3309702/6‑96 (MQ=255)
                                                                  gcgcAGTCCGGCAATAATAGGTTACAGAGTGACGCT‑TT‑TTTTACCTCTTAGAGCACGACCTGCTTTGGCGGCTGGCCTCTCGTGCCGGAACCTAgccgcc                               >  1:3229447/1‑100 (MQ=255)
                                                                    gcAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                             <  1:1169191/100‑1 (MQ=255)
                                                                     cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                            >  2:78111/1‑100 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAGAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                          <  1:1781822/100‑1 (MQ=255)
                                                                       aTCCGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                          <  2:1309686/100‑1 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAACACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                          >  1:2364045/1‑97 (MQ=255)
                                                                          cGGCAATAATAGGTTACAGTGTCACG‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                      >  1:2012538/1‑100 (MQ=255)
                                                                            gCAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                     <  2:1926245/100‑1 (MQ=255)
                                                                             cAATAATAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                    <  1:2177942/100‑1 (MQ=255)
                                                                              aataatAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                   <  2:1741576/100‑1 (MQ=255)
                                                                              aataatAGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAGAGTCGCCCGTGTTTCAtt                    <  1:578031/99‑1 (MQ=255)
                                                                               ataataGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                  <  2:3007016/100‑1 (MQ=255)
                                                                                 aataGGTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                <  1:2991051/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  1:2042574/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           <  2:1264695/100‑1 (MQ=255)
                                                                                      gTTACAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc           >  2:2213149/1‑100 (MQ=255)
                                                                                         aCAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc        >  1:2916340/1‑100 (MQ=255)
                                                                                         aCAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc        <  2:3462397/100‑1 (MQ=255)
                                                                                          cAGTGTCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt       <  2:2363586/100‑1 (MQ=255)
                                                                                            gtgtCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa     >  1:1791573/1‑100 (MQ=255)
                                                                                             tgtCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt    <  2:3333732/100‑1 (MQ=255)
                                                                                              gtCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   <  2:1325872/100‑1 (MQ=255)
                                                                                               tCACGTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGc  <  2:2580027/100‑1 (MQ=255)
                                                                                                    |                                                                                                
TTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665089‑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: ‑