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.3% 564.2 / ‑5.5 165intergenic (+267/‑13)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (78/83);  minor base T (2/1);  total (81/84)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

TTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853812‑3854005
                                                                                                 |                                                                                                
ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGc                                                                                                <  2:81030/100‑1 (MQ=255)
ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGc                                                                                                <  1:1774815/100‑1 (MQ=255)
  cgtGGCGTATATCAGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  1:2804706/3‑100 (MQ=255)
  ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  1:1747871/1‑100 (MQ=255)
  ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  2:633784/1‑100 (MQ=255)
  ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  1:3460005/1‑100 (MQ=255)
      gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGt                                                                                          <  1:1295986/100‑1 (MQ=255)
      gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGt                                                                                          >  2:689070/1‑100 (MQ=255)
       cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTg                                                                                         <  2:1912905/100‑1 (MQ=255)
          atatCTGGCTAACATTCATTAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATa                                                                                      <  2:1080280/100‑1 (MQ=255)
           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     >  1:1803568/1‑100 (MQ=255)
           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     >  1:3652186/1‑100 (MQ=255)
              cTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGa                                                                                  >  1:1340535/1‑100 (MQ=255)
                ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaa                                                                                >  2:2688981/1‑100 (MQ=255)
                ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaa                                                                                >  1:2429548/1‑100 (MQ=255)
                ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaa                                                                                >  1:3456907/1‑100 (MQ=255)
                 gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               <  1:3862005/100‑1 (MQ=255)
                  cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               >  2:1397900/1‑99 (MQ=255)
                   tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             <  1:726401/100‑1 (MQ=255)
                   tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAGCAGTGATAATGAAAAgg                                                                             >  2:3340256/1‑100 (MQ=255)
                    aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGc                                                                            >  2:4123164/1‑100 (MQ=255)
                     aCATTCATCAATGTGATAGATGCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCa                                                                           <  1:3340276/100‑1 (MQ=255)
                        ttCATCAATTTGATAGATTCCTTTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        <  2:3728559/100‑1 (MQ=255)
                        ttCATCAATTTGATAGATGCCTCTCCCGCATTTATGTGAATGCGTAGTTACTTATTCTAATTATTTTTATAAATGCATCCGTGATAATGAAAAGGCaaat                                                                        <  2:1012205/100‑2 (MQ=255)
                        ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        >  2:3358099/1‑100 (MQ=255)
                        ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        <  2:3078643/100‑1 (MQ=255)
                         tcatcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGa                                                                       <  2:3819849/100‑1 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGTCAAAGaa                                                                      >  2:3291929/1‑100 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:1156095/1‑100 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:717831/1‑100 (MQ=255)
                          cataaATGTGATAGATTCCTCTCCCGACTTTATGGGAATGTGTAGTGATTTATCATAAATATTCTTATAAAGAACTCCGTGGTAATGCAGAGGCCAAg                                                                        >  2:3631306/1‑98 (MQ=255)
                           atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     >  1:2074214/1‑100 (MQ=255)
                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                    <  2:2074201/100‑1 (MQ=255)
                             caATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  2:2427027/100‑1 (MQ=255)
                             caATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  2:553154/100‑1 (MQ=255)
                             caAGGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  1:519199/100‑1 (MQ=255)
                              aaTGTGATAGATTCCTCTCCCGCATTTATGTGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATTAAAAGGCAAAGACACGt                                                                  >  2:827022/1‑100 (MQ=255)
                              aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAATAAACGt                                                                  >  1:872104/1‑100 (MQ=255)
                               aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                    <  2:4016041/97‑1 (MQ=255)
                                 gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGTCAAAGAAACGTCaa                                                               >  1:3403932/1‑100 (MQ=255)
                                 gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCaa                                                               <  1:3169946/100‑1 (MQ=255)
                                   gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAtt                                                             >  2:3977195/1‑100 (MQ=255)
                                     tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           <  1:1121174/100‑1 (MQ=255)
                                       gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         >  2:287015/1‑100 (MQ=255)
                                           ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     >  2:1431740/1‑100 (MQ=255)
                                           ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     >  2:3456096/1‑100 (MQ=255)
                                            ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTCTTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCACTGTGTTAt                                                       >  2:1916011/1‑97 (MQ=255)
                                            ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTTTTTTTTTTTCTGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     >  2:3326396/1‑99 (MQ=255)
                                             tctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGAt                                                   >  2:1305122/1‑100 (MQ=255)
                                              ctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGGGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATg                                                  >  1:2210620/1‑100 (MQ=255)
                                              ctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGGGATAATGAAAAGGCAAAGAAACGTAAATTTGTTAttgttg                                                  >  1:2784528/1‑100 (MQ=255)
                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  1:485032/1‑100 (MQ=255)
                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGt                                                 <  2:3042728/99‑1 (MQ=255)
                                                 ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGttt                                               >  1:2343382/1‑99 (MQ=255)
                                                 ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               >  1:4098575/1‑100 (MQ=255)
                                                 ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               <  1:2663324/100‑1 (MQ=255)
                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  1:2987911/100‑1 (MQ=255)
                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              >  2:1875320/1‑100 (MQ=255)
                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTAGTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              >  2:2207094/1‑100 (MQ=255)
                                                      tttATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAtt                                          <  1:3816162/100‑1 (MQ=255)
                                                         aTGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACt                                       <  1:3677556/100‑1 (MQ=255)
                                                          tGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                      >  2:13856/1‑100 (MQ=255)
                                                            ggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt                                    >  2:3590138/1‑100 (MQ=255)
                                                            ggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt                                    >  1:2484830/1‑100 (MQ=255)
                                                              aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:1063037/100‑1 (MQ=255)
                                                              aaTGCGTAATTAATTATTCTAATTATTTTTATAAAGGCATCCATGATAATGCAACGGCAAAGAAACGTCAATTTGTTATTGATGTTAGTATTACt                                       <  2:2312253/95‑1 (MQ=255)
                                                                  cGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              <  1:3901162/100‑1 (MQ=255)
                                                                  cGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              >  2:3846397/1‑100 (MQ=255)
                                                                   gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  1:352996/100‑1 (MQ=255)
                                                                     aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgn                           <  1:2520423/100‑2 (MQ=255)
                                                                     aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  2:1550491/100‑1 (MQ=255)
                                                                       tGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTc                         <  2:1140109/100‑1 (MQ=255)
                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        >  2:859782/1‑100 (MQ=255)
                                                                         aCTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAATGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAg                       <  2:1773941/100‑1 (MQ=255)
                                                                         aCTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAg                       <  1:1710251/100‑1 (MQ=255)
                                                                          cTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATGGATGTTGGTATTACTCGTTGCCGTCGTTCata                      >  1:2678317/1‑98 (MQ=255)
                                                                           ttATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGAt                     <  1:3962710/100‑1 (MQ=255)
                                                                             aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   <  2:2922408/100‑1 (MQ=255)
                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAATAGGCATATTAACGTTAATTTGTTATTGATGTTGGTATTACTCGTGTATGTCGGTCAGTTGgtgc                >  1:428175/1‑100 (MQ=255)
                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGCTGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaa              >  1:2564090/1‑100 (MQ=255)
                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAATCGTGGCCGTCGGTCAGATGGCGCaa              >  2:2310311/1‑100 (MQ=255)
                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATAACTCGTGGCCGTCGGTGAGATGGCGCaa              >  2:2476319/1‑100 (MQ=255)
                                                                                   attattTTTATAAAGGCATCCGTGATAATGAAAATGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaa             >  2:1784864/1‑100 (MQ=255)
                                                                                   attattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaa             <  1:1122091/100‑1 (MQ=255)
                                                                                    ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  2:2020545/100‑1 (MQ=255)
                                                                                    ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  1:1949359/100‑1 (MQ=255)
                                                                                    ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            >  1:2499914/1‑100 (MQ=255)
                                                                                      attTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCa          >  1:3864591/1‑100 (MQ=255)
                                                                                      attTTGATGAAGGCATCCGTCATGATGAATAGGCAAAGAAACGTCAATTTGTTATTGATCTTCGTATTACTCGTGGCCGTCGGTCAGATGgcg                 <  2:2562329/93‑1 (MQ=255)
                                                                                       tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         >  1:3321511/1‑100 (MQ=255)
                                                                                        ttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        <  1:2441287/100‑1 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  1:1068917/100‑1 (MQ=255)
                                                                                          ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCCATTTGTTATTGCTGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      >  2:1040832/1‑100 (MQ=255)
                                                                                          ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtttn      <  1:721788/100‑2 (MQ=255)
                                                                                          ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  1:2207107/100‑1 (MQ=255)
                                                                                          ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:207776/100‑1 (MQ=255)
                                                                                           tataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtat     >  2:1314901/1‑100 (MQ=255)
                                                                                           tataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtat     >  1:1178862/1‑100 (MQ=255)
                                                                                             taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   <  1:1470317/100‑1 (MQ=255)
                                                                                              aaagGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAtt  >  1:3966062/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853812‑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: