Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,909 A→G 100% intergenic (+267/‑13) tisB → / → emrD toxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,9090AG99.4% 542.8 / ‑4.8 157intergenic (+267/‑13)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (82/73);  total (83/74)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 158 total aligned reads displayed.

TCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853813‑3854005
                                                                                                |                                                                                                
tCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCa                                                                                               >  2:751269/1‑100 (MQ=255)
  cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  2:1098924/1‑100 (MQ=255)
      cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTg                                                                                         >  1:3210874/1‑100 (MQ=255)
      cGTATATCTGGCTAACATTCATCAATGTGATAGAGTCCTCTCCCGCATTTATGGTAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTg                                                                                         >  1:1477825/1‑100 (MQ=255)
        tataTCTAGCTAACATTCATCAAAGTGATAGATTCAACTTGAGCATTTATGGGCATGAGTAGTGACTTATTCTGATTATTTTTATAAAGGCGTCCGTg                                                                                         <  1:354041/98‑1 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     <  1:3080797/100‑1 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     <  2:3770325/100‑1 (MQ=255)
            tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   >  2:1362084/1‑100 (MQ=255)
            tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  1:1090584/100‑1 (MQ=255)
             cTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGa                                                                                  <  2:2704542/100‑1 (MQ=255)
               ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  1:1090872/97‑1 (MQ=255)
               ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaa                                                                                <  2:183548/100‑1 (MQ=255)
                gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               >  2:3761457/1‑100 (MQ=255)
                gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGGAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               >  1:2924572/1‑100 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTTTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              <  2:3679400/100‑1 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  1:3153080/1‑100 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              <  2:3443834/100‑1 (MQ=255)
                 cTAACATGCATGAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  1:2244545/1‑100 (MQ=255)
                  tgacacTCATCAGTCTGATAGATCCCTCTACCGCATTTATGCGCATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              <  2:3842393/97‑1 (MQ=255)
                  tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               >  2:312874/1‑98 (MQ=255)
                  tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  2:2775353/1‑99 (MQ=255)
                    aCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGTCATCCGTGATAATGAAAAGGCa                                                                           <  1:2477770/100‑1 (MQ=255)
                     cATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaa                                                                          <  1:1386182/100‑1 (MQ=255)
                      aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaa                                                                         <  1:759825/100‑1 (MQ=255)
                      aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaa                                                                         >  2:1685371/1‑100 (MQ=255)
                        tcatcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGa                                                                       <  2:1862882/100‑1 (MQ=255)
                          atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     <  2:4277796/100‑1 (MQ=255)
                          atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     <  2:2537364/100‑1 (MQ=255)
                          atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     >  2:788241/1‑100 (MQ=255)
                           tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                    >  1:2459175/1‑100 (MQ=255)
                                gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATACGTGATAATGAAAGTGCAAAGAAACGTCaa                                                               >  1:1816655/1‑100 (MQ=255)
                                 tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAATGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAt                                                              >  1:96824/1‑100 (MQ=255)
                                 tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAt                                                              >  2:1654443/1‑100 (MQ=255)
                                  gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTGTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAtt                                                             >  2:2466472/1‑100 (MQ=255)
                                  gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAtt                                                             >  2:3600840/1‑100 (MQ=255)
                                      tccgatcttcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  1:3515651/92‑1 (MQ=255)
                                      gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  1:693935/100‑1 (MQ=255)
                                       aTTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTa                                                        <  2:2436715/100‑1 (MQ=255)
                                        ttCCTCTCCCGCATTCATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       >  1:3439633/1‑100 (MQ=255)
                                         tCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                                                      >  1:2461444/1‑100 (MQ=255)
                                          ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     <  2:2346879/100‑1 (MQ=255)
                                           ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTCATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                    >  1:3852232/1‑100 (MQ=255)
                                           ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGACACGGCAAATAAACTTCCATCTGTTATTg                                                     >  1:3852183/1‑99 (MQ=255)
                                           ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAGTTGTTATTGa                                                    >  1:2798368/1‑100 (MQ=255)
                                            tctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGACAAAGAAACGTCAATTTGTTATTGAt                                                   >  1:2439902/1‑100 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTGTAATGGTTTCATTGATAATGAAAAGACAAAGAAACTTCAATTTGTTATTGAt                                                   >  1:1030437/1‑97 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  1:4149116/1‑100 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:2590063/1‑100 (MQ=255)
                                                ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               >  2:3350406/1‑100 (MQ=255)
                                                ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               >  2:1019337/1‑100 (MQ=255)
                                                 cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTCATGTTgg                                              <  2:2244551/100‑1 (MQ=255)
                                                  gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                             >  2:2370434/1‑100 (MQ=255)
                                                   cATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTa                                            <  2:52762/100‑1 (MQ=255)
                                                     tttATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAtt                                          <  1:704054/100‑1 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        >  2:2982529/1‑100 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        <  2:1483156/100‑1 (MQ=255)
                                                        aTGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACt                                       <  1:1294796/100‑1 (MQ=255)
                                                         tGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                      >  1:3579008/1‑100 (MQ=255)
                                                         tGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                      >  1:239071/1‑100 (MQ=255)
                                                          gggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCg                                     <  1:229367/100‑1 (MQ=255)
                                                           ggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt                                    >  1:2462074/1‑100 (MQ=255)
                                                            gAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg                                   <  1:481041/100‑1 (MQ=255)
                                                              aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 >  1:2151237/1‑100 (MQ=255)
                                                               tGCGTAGTGCCTTATGCTAATTATGTTTATAAAGGCATCCGTGATAATGAAAAGACAAAGAAACGTCAATTTGTTATTGATGTTGGTATTCCTCGTGGcc                                <  1:2578013/100‑1 (MQ=255)
                                                                  gTAGTGACTTATTCTAATTATTTTTATAAAGGGATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              >  1:1587094/1‑99 (MQ=255)
                                                                  gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  1:1935420/100‑1 (MQ=255)
                                                                    aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  1:3514393/100‑1 (MQ=255)
                                                                    aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           >  1:2344310/1‑100 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        <  1:1395113/100‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        >  1:1565571/1‑100 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        <  2:2235203/100‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        <  1:3325277/100‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        <  1:2982517/100‑1 (MQ=255)
                                                                           tATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATg                    <  1:2381489/100‑1 (MQ=255)
                                                                            aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   <  2:3876679/100‑1 (MQ=255)
                                                                            aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   >  1:4239436/1‑100 (MQ=255)
                                                                             ttCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgc                  <  2:4234317/100‑1 (MQ=255)
                                                                             ttCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgc                  <  2:4070589/100‑1 (MQ=255)
                                                                              tCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgc                  >  1:979329/1‑99 (MQ=255)
                                                                               cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                <  2:1673091/100‑1 (MQ=255)
                                                                               cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  1:3306071/1‑100 (MQ=255)
                                                                               cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  1:1891946/1‑100 (MQ=255)
                                                                               cTAATTATTATTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                <  1:2671993/100‑1 (MQ=255)
                                                                                tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGTCGCa               >  2:2010532/1‑100 (MQ=255)
                                                                                 aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAATGAAATGTCCATTTGTTATTGATGTTGGTATTACACGTGGCCGTCGGTCAGATGgcg                 >  1:2093704/1‑97 (MQ=255)
                                                                                 aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaa              >  2:3630993/1‑100 (MQ=255)
                                                                                  attattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaa             <  2:4149138/100‑1 (MQ=255)
                                                                                   ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  1:2940554/100‑1 (MQ=255)
                                                                                   ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  2:3458558/100‑1 (MQ=255)
                                                                                    tattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAcc           <  1:821031/100‑1 (MQ=255)
                                                                                      tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         >  1:61153/1‑100 (MQ=255)
                                                                                       ttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        >  1:1081127/1‑100 (MQ=255)
                                                                                       ttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        >  2:1742135/1‑100 (MQ=255)
                                                                                         ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGATACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:3210886/100‑1 (MQ=255)
                                                                                         ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:1983885/100‑1 (MQ=255)
                                                                                         ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:878163/100‑1 (MQ=255)
                                                                                         atataAAGGAATCCGTGATAATGATAAGGCAAAGAAACATCAATTTGTTATTGAGGGTGGTATAACTCGTGGCCGTCGATCAGATGGCGAACACCATTta      >  1:1507655/2‑100 (MQ=255)
                                                                                          tataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtat     >  2:1005487/1‑100 (MQ=255)
                                                                                           ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata    >  2:3203424/1‑100 (MQ=255)
                                                                                           ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata    <  2:2902813/100‑1 (MQ=255)
                                                                                                |                                                                                                
TCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853813‑3854005

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: