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,9090AG97.1% 374.9 / ‑0.6 119intergenic (+267/‑13)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (54/58);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 120 total aligned reads displayed.

TTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853812‑3854005
                                                                                                 |                                                                                                
ttCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGc                                                                                                <  2:291198/100‑1 (MQ=255)
   cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  2:1076929/1‑100 (MQ=255)
   cTGGCGTATATCTGGCTAACACTAATCAATGTGATAGATTCCTCTCCCGCATTTATGTGAATGCGTCGTGACTTATTCTAATCATTTTTATAAAGGc                                                                                                >  2:1893447/1‑97 (MQ=255)
     ggCGTATATCTGGCTAACATTCATTAATGTGATAGATTCCTCTCACGCATTTATGGGAATGCGTAGTGACTTATTATAATTATTTTTATAAAGGAATCCg                                                                                           >  2:493984/1‑100 (MQ=255)
       cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATACGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTg                                                                                         <  2:80665/100‑1 (MQ=255)
         tataTCTGGCTAACATTCATCAATGTGATAGATTCTTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGa                                                                                        <  2:949826/99‑1 (MQ=255)
          atatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATa                                                                                      <  2:1756007/100‑1 (MQ=255)
           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                     <  2:1842845/100‑1 (MQ=255)
            atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTTATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:1298092/100‑1 (MQ=255)
            atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:498787/100‑1 (MQ=255)
             tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  1:786214/100‑1 (MQ=255)
               tGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaa                                                                                 >  1:582778/1‑100 (MQ=255)
                 gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                               >  2:2191067/1‑100 (MQ=255)
                  cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAACGCATCCGTGATAATGAAAAg                                                                              >  2:2273082/1‑100 (MQ=255)
                  cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAACGCATCCGTGATAATGAAAAg                                                                              <  2:1468007/100‑1 (MQ=255)
                   tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaagatcggaa                                                                             >  1:244424/1‑92 (MQ=255)
                    aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGc                                                                            >  2:168453/1‑100 (MQ=255)
                      cATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaa                                                                          >  2:1870739/1‑100 (MQ=255)
                      cATTCATCAATGTGATAGATTCATCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGTCATCCGTGATAATGAAAAGGCaa                                                                          >  2:1422939/1‑100 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:1767167/1‑100 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:1461101/1‑100 (MQ=255)
                          catcaATGTGATAGATTCCTCTCCAGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaata                                                                       >  2:730860/1‑97 (MQ=255)
                              aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGCGATAATGAAAAGGCAAAGAAACGt                                                                  >  1:2181978/1‑100 (MQ=255)
                               aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                   <  2:2196878/98‑1 (MQ=255)
                                 gtgATATCTTCCCTTCCCGCATTTATGTGAATGCTTAGTGACTTATTCTAATTATTTTTATAAAGACATCCGTCATAATGATAACGCAAAGaaatatca                                                                >  2:1183632/1‑94 (MQ=255)
                                     tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           <  1:1020012/100‑1 (MQ=255)
                                       gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  1:2078055/100‑1 (MQ=255)
                                       gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         >  1:660980/1‑100 (MQ=255)
                                         ttCCTCTCCCGCATTTATGGGAATGCGTACTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       <  2:1345529/100‑1 (MQ=255)
                                         ctCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       <  1:483466/99‑1 (MQ=255)
                                           ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtag                                                     >  2:483463/1‑98 (MQ=255)
                                           ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     <  1:799609/100‑1 (MQ=255)
                                           ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     >  2:156693/1‑100 (MQ=255)
                                            ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAATTGCATTCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                    >  2:967733/1‑100 (MQ=255)
                                             tctttcGCCTTTATGGGAATGCGTAGTGACTTATTCTAATGATTTTTATATAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                    <  2:522656/94‑1 (MQ=255)
                                             tctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGACAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGAt                                                   <  1:1687428/100‑1 (MQ=255)
                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:400510/1‑100 (MQ=255)
                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                <  1:190462/100‑1 (MQ=255)
                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:520477/1‑100 (MQ=255)
                                                 ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                               >  1:2243901/1‑100 (MQ=255)
                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:2006314/100‑1 (MQ=255)
                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              >  1:300954/1‑100 (MQ=255)
                                                   gaaTTTATGGGAATGAGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:1402084/97‑1 (MQ=255)
                                                   gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAACGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                             >  2:1151664/1‑100 (MQ=255)
                                                    cATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTa                                            <  1:1461117/100‑1 (MQ=255)
                                                     aTTTATTGGAATGCGTAGTTACTTATTCTAATTATATTTATAAAGGCATCAGTGATAGTGAAAAGGCAAAGAAACGTAAATTTGTTATTGATGttagaat                                           >  2:16226/1‑95 (MQ=255)
                                                     aTTTATGGGAATGCGTAGTGACTTATTCTAGTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAt                                           <  2:2388811/100‑1 (MQ=255)
                                                     aTTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGGGATAATGAAAAGGCAAAGAAACGCCAATTAGTTATTGATGtt                                                >  1:2116924/1‑95 (MQ=255)
                                                      tttATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAtt                                          <  2:175266/100‑1 (MQ=255)
                                                        tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGACAATTTGTTATTGATGTTGGTGtt                                          >  1:234193/1‑98 (MQ=255)
                                                             gAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg                                   <  2:583271/100‑1 (MQ=255)
                                                             gAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg                                   <  2:579699/100‑1 (MQ=255)
                                                             gAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg                                   <  1:2073314/100‑1 (MQ=255)
                                                              aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:1804513/100‑1 (MQ=255)
                                                              aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  2:2329006/100‑1 (MQ=255)
                                                               aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  1:2348261/100‑1 (MQ=255)
                                                               aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  1:2130773/100‑1 (MQ=255)
                                                               aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  2:2351340/100‑1 (MQ=255)
                                                               aTGCGTAGTGACTTATGCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  1:2246328/100‑1 (MQ=255)
                                                                tGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAGGGCAAAGAAACGTCAATTTGTTATTGATGTTGCTATTACTCGTGGc                                 >  1:1532982/1‑99 (MQ=255)
                                                                tGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  2:1322968/99‑1 (MQ=255)
                                                                 gCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCg                               <  2:402035/100‑1 (MQ=255)
                                                                  cGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              >  1:1969146/1‑100 (MQ=255)
                                                                   gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  1:2401110/100‑1 (MQ=255)
                                                                   gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             >  1:1423393/1‑100 (MQ=255)
                                                                    tAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              >  2:723345/1‑98 (MQ=255)
                                                                    tAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             >  1:1351581/1‑99 (MQ=255)
                                                                     aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  1:1795639/100‑1 (MQ=255)
                                                                      gttaCTTATTCTAATTATTTTTATAAAGGCAGCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  2:1296519/96‑1 (MQ=255)
                                                                       tGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  1:309665/98‑1 (MQ=255)
                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        >  2:35672/1‑100 (MQ=255)
                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        >  2:931820/1‑100 (MQ=255)
                                                                            tATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATg                    >  2:1936138/1‑100 (MQ=255)
                                                                            tATTCTAATTATTTTTATAAAAGCATCCGTAATAATTAAATGGCAAATTAACGTCAATTTGTTATATATGTTTGTATTACTCGTGGCCGTATGTCAg                       >  2:1078961/1‑97 (MQ=255)
                                                                             aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGAtcg                   >  1:1211716/1‑98 (MQ=255)
                                                                              ttCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGAATTACTCTTGTTCGTCGGTCAGATg                    >  1:946162/1‑98 (MQ=255)
                                                                               tCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcg                 <  1:2297098/100‑1 (MQ=255)
                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                <  2:1241193/100‑1 (MQ=255)
                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  2:2354397/1‑100 (MQ=255)
                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcg                 >  2:202978/1‑99 (MQ=255)
                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcg                 >  1:1968226/1‑99 (MQ=255)
                                                                                 tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               <  1:1029238/100‑1 (MQ=255)
                                                                                 tAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               >  2:1859872/1‑100 (MQ=255)
                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaa              >  1:120644/1‑100 (MQ=255)
                                                                                   attattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGAGCaaa             <  1:1853209/100‑1 (MQ=255)
                                                                                    ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  1:1051856/100‑1 (MQ=255)
                                                                                      attTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGATATTGATGTTGGTATTACTCGCGGCCGTCGGTCAGATGGCGCAAACCa          >  2:131154/1‑100 (MQ=255)
                                                                                       tttttATAAAGGCATCCGTGATAATGTAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         >  1:1459172/1‑100 (MQ=255)
                                                                                       tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTCGCCGTCGGTCAGATGGCGCAAACCAt         <  2:689193/100‑1 (MQ=255)
                                                                                        ttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        <  1:392195/100‑1 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  1:400513/100‑1 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  2:830488/100‑1 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       >  2:1425046/1‑100 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       >  2:1755280/1‑100 (MQ=255)
                                                                                         tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  2:2116900/100‑1 (MQ=255)
                                                                                          ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  2:323691/100‑1 (MQ=255)
                                                                                            ataAAGGCATCCGTGTTAATGAACAGCGAATGACACGTCAATTTGGTATTGATGTTGGTATTAATCGTAGCCGTCGGTCAGATGGCGCATACCATTtat     <  2:1648097/99‑1 (MQ=255)
                                                                                             taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  2:1160360/1‑100 (MQ=255)
                                                                                             taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  1:1429195/1‑100 (MQ=255)
                                                                                              aaagGCATCCGTGTTAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAtt  >  2:1411397/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: