Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,853,909 A→G 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,9090AG93.5% 297.9 / 4.9 93intergenic (+267/‑13)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (46/41);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853811‑3854005
                                                                                                  |                                                                                                
aTTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAg                                                                                                 <  1:1581430/100‑1 (MQ=255)
  tCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCa                                                                                               <  2:1375216/100‑1 (MQ=255)
   ccTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCAt                                                                                              >  2:267595/1‑100 (MQ=255)
    cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             <  1:1719248/100‑1 (MQ=255)
    cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  2:1463332/1‑100 (MQ=255)
     tGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATcc                                                                                            <  2:384555/100‑1 (MQ=255)
      ggCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCg                                                                                           >  2:1929017/1‑100 (MQ=255)
       gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGt                                                                                          <  2:1380200/100‑1 (MQ=255)
       gCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGt                                                                                          <  2:1515740/100‑1 (MQ=255)
         gTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGa                                                                                        >  1:415388/1‑100 (MQ=255)
         gTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGa                                                                                        <  1:1090248/100‑1 (MQ=255)
          tataTCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGAt                                                                                       <  1:1147973/100‑1 (MQ=255)
             atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:1876255/100‑1 (MQ=255)
             atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:1396640/100‑1 (MQ=255)
             atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:1053558/100‑1 (MQ=255)
              tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAtn                                                                                   <  2:1607147/100‑2 (MQ=255)
              tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   >  1:1117265/1‑100 (MQ=255)
              tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  2:54296/100‑1 (MQ=255)
               cTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGa                                                                                  >  1:1698830/1‑100 (MQ=255)
                tGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaa                                                                                 >  1:589601/1‑100 (MQ=255)
                   cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGCCTTATTCTACTTATTTTTATAAAGGCATACGTGGTAATGaaaa                                                                               >  1:717817/1‑99 (MQ=255)
                   cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  1:363268/1‑100 (MQ=255)
                   cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAg                                                                              >  1:1906639/1‑100 (MQ=255)
                    tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             >  1:661959/1‑100 (MQ=255)
                      aCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCa                                                                           >  2:1890542/1‑100 (MQ=255)
                       cATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaa                                                                          >  1:818105/1‑100 (MQ=255)
                        aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaa                                                                         <  1:823733/100‑1 (MQ=255)
                            atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     >  2:1081894/1‑100 (MQ=255)
                             tccaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCTGTGATAATGAAAAGGCAAAGa                                                                       >  1:444445/4‑97 (MQ=255)
                               aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGt                                                                  >  1:70865/1‑100 (MQ=255)
                                  gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTACTCTAATTATTTTTATAAAGGCATTCGTGATAATGAAATGGCAAAGaaa                                                                     >  1:573828/1‑94 (MQ=255)
                                  gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGGAGTGACTTATTCTCAGTATTTTTATAAAGGCAAGCGTGATAATGAAAAGCCAAAGACACGACa                                                                >  1:930177/1‑99 (MQ=255)
                                   tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAACGGCAAAGAAACGTCaa                                                               >  1:1616063/1‑99 (MQ=255)
                                    gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAt                                                              >  1:376126/1‑99 (MQ=255)
                                      tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           <  2:70865/100‑1 (MQ=255)
                                      tAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTg                                                           >  2:235485/1‑100 (MQ=255)
                                       aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                          <  1:1344339/100‑1 (MQ=255)
                                       aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                          <  2:1400867/100‑1 (MQ=255)
                                       aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                          <  2:1619201/100‑1 (MQ=255)
                                       aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                          >  1:377116/1‑100 (MQ=255)
                                        gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  2:314017/100‑1 (MQ=255)
                                        gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  1:1767815/100‑1 (MQ=255)
                                            ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     <  1:1865356/100‑1 (MQ=255)
                                            ccTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTg                                                     <  1:1528512/100‑1 (MQ=255)
                                             ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                    <  2:1787592/100‑1 (MQ=255)
                                              tctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGAt                                                   >  2:364602/1‑100 (MQ=255)
                                               ctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATg                                                  <  1:1678258/100‑1 (MQ=255)
                                                 cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                <  1:267593/100‑1 (MQ=255)
                                                 cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:1462135/1‑100 (MQ=255)
                                                   cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  1:1458314/100‑1 (MQ=255)
                                                   cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:1486188/100‑1 (MQ=255)
                                                   cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:1427332/100‑1 (MQ=255)
                                                    gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                             <  2:831953/100‑1 (MQ=255)
                                                    gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                             <  2:201248/100‑1 (MQ=255)
                                                         tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        <  1:298524/100‑1 (MQ=255)
                                                         tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                        >  2:1817564/1‑100 (MQ=255)
                                                          aTGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACt                                       >  2:918089/1‑100 (MQ=255)
                                                            gggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAACGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                      >  2:638369/1‑99 (MQ=255)
                                                             ggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt                                    <  2:1436814/100‑1 (MQ=255)
                                                              gAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTg                                   <  2:811333/100‑1 (MQ=255)
                                                               aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:128347/100‑1 (MQ=255)
                                                                aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  1:9836/100‑1 (MQ=255)
                                                                    gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTCGTTATTGATGTTGGTATTACTCGTGGCCGTc                             >  2:1624044/1‑100 (MQ=255)
                                                                      aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  2:1542850/100‑1 (MQ=255)
                                                                      aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           >  2:702853/1‑100 (MQ=255)
                                                                       atgaattaaacTAATTATTTTTATAAAAGCAGCGGAAATACTGAAAAGGTACAGAAACGTCAATAAGTTATTGACGTTGGTATGACGCGAGGCAGCCGGt                          >  1:184839/11‑100 (MQ=255)
                                                                         gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                        >  2:1387546/1‑100 (MQ=255)
                                                                          aCTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAg                       >  1:1750564/1‑100 (MQ=255)
                                                                           cTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGa                      >  1:979093/1‑100 (MQ=255)
                                                                           cTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGa                      >  2:349440/1‑100 (MQ=255)
                                                                             tATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATg                    <  1:702852/100‑1 (MQ=255)
                                                                                tCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcg                 >  1:1240892/1‑100 (MQ=255)
                                                                                 cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  1:435794/1‑100 (MQ=255)
                                                                                 cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                <  1:235483/100‑1 (MQ=255)
                                                                                 cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                <  1:649865/100‑1 (MQ=255)
                                                                                 cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGCTATCGATGTTGGTATTACTCGTGGCCTCCGGTCAGATCgcgc                >  2:635357/1‑100 (MQ=255)
                                                                                  tAATTTTTTTTATAAAGGCATCCGTGATAAGTAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcg                 <  2:184840/98‑1 (MQ=255)
                                                                                  tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               >  1:542225/1‑100 (MQ=255)
                                                                                  tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGATGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  1:1393799/1‑99 (MQ=255)
                                                                                  tAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               <  2:1906652/100‑1 (MQ=255)
                                                                                   aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGATGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  1:412138/1‑98 (MQ=255)
                                                                                     ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCTGTCAGATGGTGCAAAc            >  2:1734490/1‑100 (MQ=255)
                                                                                     ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            <  2:559868/100‑1 (MQ=255)
                                                                                      tattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGAGGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc            >  2:638272/1‑99 (MQ=255)
                                                                                       attTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCa          >  2:1200367/1‑100 (MQ=255)
                                                                                       attTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCa          >  1:34665/1‑100 (MQ=255)
                                                                                        tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         >  2:1665867/1‑100 (MQ=255)
                                                                                         ttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAtt        <  2:377118/100‑1 (MQ=255)
                                                                                          tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  1:1462127/100‑1 (MQ=255)
                                                                                          tttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAttt       <  2:950789/100‑1 (MQ=255)
                                                                                             ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata    >  1:132143/1‑100 (MQ=255)
                                                                                              taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  2:1460634/1‑100 (MQ=255)
                                                                                              taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGGTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  2:1914611/1‑100 (MQ=255)
                                                                                               aaagGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAtt  >  2:653834/1‑100 (MQ=255)
                                                                                                  |                                                                                                
ATTCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853811‑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: