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,9090AG98.1% 381.7 / ‑3.0 120intergenic (+267/‑13)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (52/63);  total (55/65)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853813‑3854005
                                                                                                |                                                                                                
tCCTGTCGTATATCTGGCTAACATTCATCAATGTGATAGATTTCTCTCTCGCATTTATTGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCa                                                                                               <  1:1770279/100‑1 (MQ=255)
  cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             <  1:1126535/100‑1 (MQ=255)
  cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  1:1152596/1‑100 (MQ=255)
   tGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  1:675158/1‑98 (MQ=255)
      cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCg                                                                                           >  1:2363988/1‑98 (MQ=255)
        tataTCTGGCTAACATTCATCAATGTGATAGGTTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAATGCATCCGTGAt                                                                                       >  2:2663478/1‑100 (MQ=255)
        tataTCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGAt                                                                                       <  2:1628078/100‑1 (MQ=255)
        tataTCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGAt                                                                                       <  1:2594080/100‑1 (MQ=255)
         atatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATa                                                                                      >  1:1486503/1‑100 (MQ=255)
         atatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATa                                                                                      >  1:705753/1‑100 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     >  2:1895420/1‑100 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     >  2:336972/1‑100 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     <  2:386780/100‑1 (MQ=255)
          tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     <  1:2506717/100‑1 (MQ=255)
            tCTGGCTAACATTCATCAATGTGATAGATTCTTTTTCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  1:241393/100‑1 (MQ=255)
             cTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGTCATCCGTGATAATGa                                                                                  >  1:260723/1‑100 (MQ=255)
             cTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   >  2:1461130/1‑99 (MQ=255)
                gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               <  2:1549367/100‑1 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  1:2321525/1‑100 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              <  2:1583423/100‑1 (MQ=255)
                 cTAACATTCATCAATGTGATAGATTCATCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  2:537733/1‑100 (MQ=255)
                  tctcattcattAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATCATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             <  2:1995750/97‑1 (MQ=255)
                  tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             <  2:1883458/100‑1 (MQ=255)
                  tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             <  2:500848/100‑1 (MQ=255)
                   aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGc                                                                            <  2:1158537/100‑1 (MQ=255)
                     cATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGAGATAATGAAAAGGCaa                                                                          >  1:659626/1‑100 (MQ=255)
                     cATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTAAAAAGACATCCGTGATAATGAAAAGACa                                                                           >  1:854987/1‑99 (MQ=255)
                       ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        <  2:1232677/100‑1 (MQ=255)
                       ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        <  1:871801/100‑1 (MQ=255)
                         catcaATGTGATAGATTCCTCTACCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:2035963/1‑100 (MQ=255)
                           tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                    <  1:1913683/100‑1 (MQ=255)
                           tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                    >  2:2390780/1‑100 (MQ=255)
                            caATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  2:333839/100‑1 (MQ=255)
                            caATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  2:602508/100‑1 (MQ=255)
                              aTGTGATAGATTCCTCTCCCTCATTTATGGGAATGCGTAGTGTCTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                                                 <  1:1618433/100‑1 (MQ=255)
                              aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                                                 >  1:2233587/1‑100 (MQ=255)
                                  gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAtt                                                             >  1:5548/1‑100 (MQ=255)
                                  gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAt                                                              >  1:685031/1‑99 (MQ=255)
                                  gATAGATTCCTCTCCCGCATTTATGGGAAAGCGTAGTGACTTAGTCCAATTATTTTTATAAAGGCATACGTGATAATGAAAAGGCAAAGAAAcaccaa                                                               >  2:819845/1‑93 (MQ=255)
                                   aTAGATTCCTCTCCCTCATTTATGGGAATGCTTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGATAAGGCAAAGAAACGTCAAttt                                                            <  2:2703030/100‑1 (MQ=255)
                                    tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAcataaa                                                               >  1:1526412/1‑91 (MQ=255)
                                    tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           >  1:667984/1‑100 (MQ=255)
                                    tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           <  1:2763009/100‑1 (MQ=255)
                                      gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  2:2785490/100‑1 (MQ=255)
                                        ttCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       >  1:2098257/1‑100 (MQ=255)
                                        ttCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       <  2:2396663/100‑1 (MQ=255)
                                        ttCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       >  2:774197/1‑100 (MQ=255)
                                         tCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                                                      <  1:1277394/100‑1 (MQ=255)
                                           ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTGATTATTTTTATAACAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                    >  2:1870671/1‑100 (MQ=255)
                                             ctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGAAAAGAAACGTTAATTTGTTATTGa                                                    >  1:2687152/1‑98 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTCTTCTAATTATTTTTATAAAGGCCTCCGTGATAATGAAAAGGCAAAGCAACGTCAATTTGTTATTGATg                                                  >  1:2843816/1‑98 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:1972168/1‑100 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                <  2:581250/100‑1 (MQ=255)
                                               cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  1:452985/1‑100 (MQ=255)
                                                ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               >  1:2828172/1‑100 (MQ=255)
                                                 cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:775458/100‑1 (MQ=255)
                                                 cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:2928679/100‑1 (MQ=255)
                                                 cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:2683447/100‑1 (MQ=255)
                                                  gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  1:14303/99‑1 (MQ=255)
                                                  gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                             <  2:1467639/100‑1 (MQ=255)
                                                   cATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTa                                            >  2:368564/1‑100 (MQ=255)
                                                      ttATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTa                                         >  1:127343/1‑100 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAATGGCAAAGAAACGAAAATTTGTTATTGATGTTGGTATTAc                                        >  1:1823324/1‑100 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        <  2:2177340/100‑1 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        >  2:2967812/1‑100 (MQ=255)
                                                       tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGACAATTTGTTATTAAAGTTATTATTAc                                        >  1:2377046/1‑100 (MQ=255)
                                                         tGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGTAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                      <  2:1470818/100‑1 (MQ=255)
                                                          gggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCg                                     >  2:1128190/1‑100 (MQ=255)
                                                             aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:2785800/100‑1 (MQ=255)
                                                             aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:1851941/100‑1 (MQ=255)
                                                             aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  2:1486487/100‑1 (MQ=255)
                                                             aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCATGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  >  2:2242348/1‑100 (MQ=255)
                                                              aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  2:667980/100‑1 (MQ=255)
                                                              aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  1:3022287/100‑1 (MQ=255)
                                                               tGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGcc                                >  2:2211379/1‑100 (MQ=255)
                                                               tGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGcc                                <  1:1348642/100‑1 (MQ=255)
                                                                  gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  2:922131/100‑1 (MQ=255)
                                                                  gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  1:2361612/100‑1 (MQ=255)
                                                                  gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGATATTACTCGTGGCCg                               >  1:384039/1‑98 (MQ=255)
                                                                   tAGTTACTTATTCTAATTCTTTTTATAAAGGCATCCGTGATAATTCAAAGGCATAGAAACGTCTATTTGTTATTGATGTTGGTATTACTCg                                     <  2:2647492/91‑1 (MQ=255)
                                                                      tGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTc                         <  2:1558073/100‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTc                         <  2:240390/99‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        <  2:2508462/100‑1 (MQ=255)
                                                                       gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGAAGTTAATATTACTCGTGGCCGTCGCTCa                        >  1:861759/1‑100 (MQ=255)
                                                                        aCTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAg                       >  2:1232498/1‑100 (MQ=255)
                                                                           tATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATg                    >  2:314731/1‑100 (MQ=255)
                                                                            aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   <  1:2211402/100‑1 (MQ=255)
                                                                            aTTCTAATTATTTTTATAAAAGAAACCGTGAAAATGAAAAGGCGAAGACAAGTAAATCGGATATTGATGTTAGTACGACTCGTGGCAGTCGGTCAGAt                     >  2:298260/1‑98 (MQ=255)
                                                                             ttCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   >  2:652430/1‑99 (MQ=255)
                                                                               cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  2:716133/1‑100 (MQ=255)
                                                                               cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATAACTCATAGCCGTCAGACAGATGgcgc                >  1:1106187/1‑100 (MQ=255)
                                                                                  attattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaa             <  1:2121439/100‑1 (MQ=255)
                                                                                   ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  2:901730/100‑1 (MQ=255)
                                                                                      tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCCGATGGCGCAAACCAt         >  1:3000475/1‑100 (MQ=255)
                                                                                       ttttATAAATGCATCCGTGATAATGAAAACGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        <  2:1852141/100‑1 (MQ=255)
                                                                                         ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:815994/100‑1 (MQ=255)
                                                                                          tataAAGGCATCGGTGATAATGAAAAGTCAAAGAGACTGCAATTTGTTATTGAGGGTCGTATTTCTCGTGGTCATCGGTTAGATGGCGCAAACCATTta      <  2:1646188/99‑1 (MQ=255)
                                                                                           ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata    <  2:1709611/100‑1 (MQ=255)
                                                                                            taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  1:2149711/1‑100 (MQ=255)
                                                                                             aaagGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAtt  <  1:2812358/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: