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

Only 100 of 103 total aligned reads displayed.

CCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665091‑1665282
                                                                                                  |                                                                                               
ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑t                                                                                              >  2:2982656/1‑100 (MQ=255)
 cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tt                                                                                             >  1:1508503/1‑100 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttt                                                                                         <  2:1616080/100‑1 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttt                                                                                         >  1:3483654/1‑100 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttt                                                                                         >  1:3540439/1‑100 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttt                                                                                         <  1:832596/100‑1 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttt                                                                                         <  2:4104808/100‑1 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑ttttt                                                                                          <  1:3344864/99‑1 (MQ=255)
      gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑ttttttt                                                                                        >  2:3220723/1‑99 (MQ=255)
      gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑ttttttt                                                                                        >  1:182699/1‑99 (MQ=255)
      gcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑ttttt                                                                                          <  2:3942264/98‑1 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCGCCGTTATTTTGCCGGATGGTGACGATTACGTGCAATGCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttttat                                                                                      <  2:589624/100‑4 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttttat                                                                                      >  2:636096/1‑97 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttttatc                                                                                     <  1:1825523/100‑5 (MQ=255)
          ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttttatc                                                                                     <  1:471948/99‑5 (MQ=255)
           tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTAtctc                                                                                    <  1:561775/100‑1 (MQ=255)
            aaaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATctct                                                                                   <  1:3939234/100‑1 (MQ=255)
             aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCtt                                                                                  >  1:2398669/1‑100 (MQ=255)
              aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTa                                                                                 >  1:3658879/1‑100 (MQ=255)
              aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTa                                                                                 >  1:1824683/1‑100 (MQ=255)
              aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t‑tttttttatctttt                                                                                 >  2:2429411/1‑91 (MQ=255)
               ccGCTGTTTCGATAGGATCAACGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTaa                                                                                <  1:3875931/99‑1 (MQ=255)
                cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTaaa                                                                               <  1:112126/100‑1 (MQ=255)
                cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAg                                                                              <  2:2995965/100‑1 (MQ=255)
                   tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgca                                                                            >  1:414690/1‑100 (MQ=255)
                   tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAgcac                                                                           >  1:2750411/1‑100 (MQ=255)
                    ntttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcac                                                                           >  1:1110091/2‑100 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcac                                                                           >  2:2430635/1‑100 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACTATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgc                                                                         <  1:2897409/100‑1 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgc                                                                         >  1:2628973/1‑100 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgc                                                                         <  1:1673668/100‑1 (MQ=255)
                        cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgcac                                                                       <  1:835987/100‑1 (MQ=255)
                         gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  1:2430154/1‑100 (MQ=255)
                         gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:2478665/1‑100 (MQ=255)
                         gATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACt                                                                      >  2:67943/1‑100 (MQ=255)
                         gATAGGATCACCGTAATATTGCCGGATGGTGACGATGACGCGCAATCCGGCAATACTAGGTTACAGTGTGACGTT‑TTTTTTATATCTTAAAGAACGCACt                                                                      >  2:890486/1‑100 (MQ=255)
                            aGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGc                                                                    <  1:642308/99‑1 (MQ=255)
                             ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCtt                                                                  >  1:3295231/1‑100 (MQ=255)
                             ggATCACAGTACTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTGACGTT‑TTTTTTATCTCTTAAATCACGCACt                                                                      >  2:2170601/1‑96 (MQ=255)
                              gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                 >  1:902280/1‑100 (MQ=255)
                              gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                 <  1:884812/100‑1 (MQ=255)
                              gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                 >  1:2542550/1‑100 (MQ=255)
                              gATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                 >  1:1779341/1‑100 (MQ=255)
                               aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCtttt                                                                >  2:1531749/1‑100 (MQ=255)
                                  aCCGTAATATTGCCGGATGGTGACGATTACGTGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             >  1:312427/1‑100 (MQ=255)
                                  aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  1:1614245/100‑1 (MQ=255)
                                  aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                             <  2:209929/100‑1 (MQ=255)
                                   ncGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTCCTTTTGCgg                                                            >  2:1482014/2‑100 (MQ=255)
                                   ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            >  2:4050219/1‑100 (MQ=255)
                                   ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                            <  1:1564805/100‑1 (MQ=255)
                                    cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                           >  2:2655127/1‑100 (MQ=255)
                                     gTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCt                                                          >  2:1956454/1‑100 (MQ=255)
                                       aaTATTGCCGGATGGTGACGATTACTCGCACTCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  2:3325601/100‑1 (MQ=255)
                                       aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  2:1971611/100‑1 (MQ=255)
                                       aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                        <  2:2310500/100‑1 (MQ=255)
                                        atatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                       <  2:182698/100‑1 (MQ=255)
                                         tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                      <  1:764381/100‑1 (MQ=255)
                                            tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCTCGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                   <  2:1106729/100‑1 (MQ=255)
                                            tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                   <  1:2478679/100‑1 (MQ=255)
                                             gCCGGATGGTGACGATTACGCGCAATCCGGCAATCATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAACGCCAGAACTGCTTTTGCGGCTGGcgtatt                                                  >  1:2828364/1‑95 (MQ=255)
                                             gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCTGCTGGc                                                       >  2:1117949/1‑95 (MQ=255)
                                                  aTGGTGACGATTACGCGCAATCCGGCAATAATAGGTGACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCGTTTGCGGCTGGCCTCTGTGGcc                                             <  2:1957241/100‑1 (MQ=255)
                                                   tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                            <  1:3219890/100‑1 (MQ=255)
                                                   tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                            <  1:3111001/100‑1 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                           >  1:759047/1‑100 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                           <  2:420122/100‑1 (MQ=255)
                                                     gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                          >  2:118970/1‑100 (MQ=255)
                                                      tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                         >  1:994784/1‑100 (MQ=255)
                                                      tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                         >  1:1463589/1‑100 (MQ=255)
                                                      tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                          <  1:2648160/99‑1 (MQ=255)
                                                       gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                        >  1:1866416/1‑100 (MQ=255)
                                                           aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                    >  1:1560149/1‑100 (MQ=255)
                                                               cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  1:1266859/1‑100 (MQ=255)
                                                               cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  1:369295/1‑100 (MQ=255)
                                                               cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                                >  1:1439437/1‑100 (MQ=255)
                                                                gcgcACTCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAACGCACGCACTGCTTTCGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                               <  2:2731761/100‑1 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTTCCGCAAAATaat                                   >  2:1437017/1‑94 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                               <  1:422746/100‑1 (MQ=255)
                                                                 cgcAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGccc                              >  2:3464750/1‑100 (MQ=255)
                                                                 cgcAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGccc                              >  1:491814/1‑100 (MQ=255)
                                                                   cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                            <  1:987716/100‑1 (MQ=255)
                                                                   cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg                             >  1:3266484/1‑99 (MQ=255)
                                                                       ccGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                        <  1:2291916/100‑1 (MQ=255)
                                                                        cGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                       >  2:1800880/1‑100 (MQ=255)
                                                                        cGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTc                       >  1:3726918/1‑100 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                      <  1:2692132/100‑1 (MQ=255)
                                                                              taataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                 <  1:1979691/100‑1 (MQ=255)
                                                                               aataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                >  1:2837518/1‑100 (MQ=255)
                                                                               aataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                >  1:2583964/1‑100 (MQ=255)
                                                                               aataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACGGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                 <  1:132955/99‑1 (MQ=255)
                                                                                   ggTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCCAAATAGTCGCCCGTGTTTCATTGCCCAttt            >  1:12434/1‑100 (MQ=255)
                                                                                      tACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg         <  2:3981785/100‑1 (MQ=255)
                                                                                      tACACTGTCTCG‑TGTTGTTTATCTCTTCAAGGACGCTCTGCTTTTGCGGCTCGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAtt             <  2:1071793/96‑1 (MQ=255)
                                                                                        cAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCt       >  1:2704290/1‑100 (MQ=255)
                                                                                        cAGTGCCACGTT‑TTTTTTATCTCTTAATGCACGCATTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg         <  2:1655156/98‑1 (MQ=255)
                                                                                          gtgtCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa     >  1:3897660/1‑100 (MQ=255)
                                                                                            gtCACGTT‑TTTTTTATTTTTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCTCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt    >  2:2412487/1‑99 (MQ=255)
                                                                                            gtCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   >  1:2483441/1‑100 (MQ=255)
                                                                                            gtCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg   >  1:2482881/1‑100 (MQ=255)
                                                                                             tCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGc  <  1:4050246/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGC  >  NC_000913/1665091‑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: ‑