Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,767 T→C 16.7% intergenic (+59/+40) acrD → / ← ypfM aminoglycoside/multidrug efflux system/stress‑induced small enterobacterial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,590,7670TC16.7% 345.9 / 57.3 155intergenic (+59/+40)acrD/ypfMaminoglycoside/multidrug efflux system/stress‑induced small enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base T (65/64);  new base C (7/19);  total (72/83)
Fisher's exact test for biased strand distribution p-value = 3.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 158 total aligned reads displayed.

GGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAT  >  NC_000913/2590669‑2590866
                                                                                                  |                                                                                                   
ggTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTg                                                                                                    >  2:2407198/1‑100 (MQ=255)
 gTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGc                                                                                                   <  1:3347816/100‑1 (MQ=255)
   gcgcCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTc                                                                                                 >  1:3383351/1‑100 (MQ=255)
   gcgcCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTc                                                                                                 >  2:2564048/1‑100 (MQ=255)
    cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCa                                                                                                >  2:988859/1‑100 (MQ=255)
     gccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCACTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCt                                                                                                  >  1:3613074/1‑97 (MQ=255)
     gccgccgcTTCCCGCTGAAACCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAg                                                                                               >  2:152140/1‑100 (MQ=255)
      ccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTACACTTTCCGATTAAGAACCTGCTCa                                                                                                >  1:399097/1‑98 (MQ=255)
       cgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCg                                                                                             >  1:3233629/1‑100 (MQ=255)
        gccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCgg                                                                                            <  1:3860752/100‑1 (MQ=255)
        gccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCgg                                                                                            <  1:3070093/100‑1 (MQ=255)
         ccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                           >  2:2396005/1‑100 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                       <  1:2100073/100‑1 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                       <  1:920444/100‑1 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCt                                                                                       <  2:4337182/100‑1 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCt                                                                                       <  2:4186193/100‑1 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTACCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCtt                                                                                      <  1:2280823/100‑1 (MQ=255)
             ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGGGTACCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCtt                                                                                      <  1:2435310/100‑1 (MQ=255)
               cccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                       <  1:2952151/98‑1 (MQ=255)
                ccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                    <  2:894299/100‑1 (MQ=255)
                ccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                    >  2:1310327/1‑100 (MQ=255)
                 cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCt                                                                                   >  1:2727450/1‑100 (MQ=255)
                 cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaat                                                                                   <  1:3990941/100‑4 (MQ=39)
                  gCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                  <  1:3884580/100‑5 (MQ=38)
                  gCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                  <  1:3331286/100‑5 (MQ=38)
                  gCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                  <  1:1907068/100‑5 (MQ=38)
                      aaGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTg                                                                              <  2:1888122/100‑1 (MQ=255)
                         ccGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTAc                                                                           >  1:1514702/1‑100 (MQ=255)
                          cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTAc                                                                           <  1:35417/99‑1 (MQ=255)
                          cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACt                                                                          >  2:2887936/1‑100 (MQ=255)
                          cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACt                                                                          <  2:1446938/100‑1 (MQ=255)
                          cgcgCCCGGAATAAGCAATAAAAAAGGCGAAATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAACAACCTGCTACGCGGGTTCTTGCTATTTGTAc                                                                           >  1:351973/1‑99 (MQ=255)
                            cgcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGTTTCTTGCTCTTTGTACttt                                                                        >  1:614676/1‑100 (MQ=255)
                             gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTg                                                                       <  1:1175504/100‑1 (MQ=255)
                                  gAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCa                                                                  <  2:3991327/100‑1 (MQ=255)
                                    aTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAgg                                                                <  1:4340994/100‑1 (MQ=255)
                                        gCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtt                                                            >  1:2008518/1‑100 (MQ=255)
                                        gCAATAAAAAAGGCGACATGCCAATGTGTCGCCTATTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtt                                                            >  1:1787469/1‑100 (MQ=255)
                                         cAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACATTGTCTCAGGAAtta                                                           >  2:2998234/1‑100 (MQ=255)
                                          aaTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttac                                                          <  2:3147292/100‑1 (MQ=255)
                                             aaaaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttactta                                                       <  2:3291158/100‑1 (MQ=255)
                                              aaaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGAGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttacttac                                                      >  1:2388296/1‑100 (MQ=255)
                                               aaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTTCTGTTTGTTCTTTGTCTCAGGAATTACTTACg                                                     <  1:1171832/100‑1 (MQ=255)
                                                  aGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAg                                                  <  2:2008526/100‑1 (MQ=255)
                                                   ggCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGc                                                 >  1:2091225/1‑100 (MQ=255)
                                                   ggCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGc                                                 >  2:1147450/1‑100 (MQ=255)
                                                          tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCCGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACtt                                          >  1:3405119/1‑100 (MQ=255)
                                                          tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACtt                                          >  1:313996/1‑100 (MQ=255)
                                                          cattggcatgtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACtt                                          >  1:3626973/9‑100 (MQ=255)
                                                           gCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTc                                         >  1:645263/1‑100 (MQ=255)
                                                            ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGCACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTATCAGGAATTACTTACGAAGCATAACTTc                                         >  1:2294808/1‑99 (MQ=255)
                                                            ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                        <  2:3616522/100‑1 (MQ=255)
                                                            ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                        >  2:1417886/1‑100 (MQ=255)
                                                            ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                        <  1:4329798/100‑1 (MQ=255)
                                                               aTGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGAt                                      <  1:2963263/99‑1 (MQ=255)
                                                               aTGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATa                                     <  1:4221962/100‑1 (MQ=255)
                                                                  tgtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAg                                  <  2:2731917/100‑1 (MQ=255)
                                                                   gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                 >  1:3692676/1‑100 (MQ=255)
                                                                   gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                 >  2:1030432/1‑100 (MQ=255)
                                                                   gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                 <  1:1883008/100‑1 (MQ=255)
                                                                     cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCt                               >  2:389568/1‑100 (MQ=255)
                                                                       ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCttt                             <  2:2728186/100‑1 (MQ=255)
                                                                         tttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                           <  2:3502708/100‑1 (MQ=255)
                                                                            ttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                        >  2:2791384/1‑100 (MQ=255)
                                                                             tCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                        <  1:2581240/99‑1 (MQ=255)
                                                                              cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTc                      <  2:4314315/100‑1 (MQ=255)
                                                                              cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTc                      >  2:4012475/1‑100 (MQ=255)
                                                                               aaCGTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                     <  2:3962498/100‑1 (MQ=255)
                                                                                     ccGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGtt               <  2:522136/100‑1 (MQ=255)
                                                                                      cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc              >  2:3239880/1‑100 (MQ=255)
                                                                                       cagcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCa             >  2:370700/5‑100 (MQ=38)
                                                                                       cagcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGCTCa             >  2:1617830/5‑100 (MQ=38)
                                                                                        agcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCAc            >  2:1289822/4‑100 (MQ=38)
                                                                                        agcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCAc            >  1:242465/4‑100 (MQ=38)
                                                                                        aTTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCAc            <  2:1112333/100‑1 (MQ=255)
                                                                                         gcaagaagCCGCTGAGCAGGTTATTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAGGCATAACTTCGATAAAGTGTTTCAAGATACCCAGGTc              >  1:1998820/3‑97 (MQ=255)
                                                                                            aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc        <  2:3550053/100‑1 (MQ=255)
                                                                                            aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc        <  2:1390516/100‑1 (MQ=255)
                                                                                            aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc        <  2:1514707/100‑1 (MQ=255)
                                                                                            aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCACAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc        >  1:2600984/1‑100 (MQ=255)
                                                                                            aGAACCCGCGGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGGTAAAGTCTTTCCAGTTCCCCAGTTCACGTTc        >  1:3920861/1‑100 (MQ=255)
                                                                                             gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa       >  1:1417053/1‑100 (MQ=255)
                                                                                             gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa       <  1:2852285/100‑1 (MQ=255)
                                                                                             gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa       <  2:2802200/100‑1 (MQ=255)
                                                                                             gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa       >  2:3059822/1‑100 (MQ=255)
                                                                                              aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCaa      <  2:1973256/100‑1 (MQ=255)
                                                                                              aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCaa      >  1:937631/1‑100 (MQ=255)
                                                                                               aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt     <  2:2702881/100‑1 (MQ=255)
                                                                                               aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt     >  2:2179904/1‑100 (MQ=255)
                                                                                                ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt     <  1:692626/99‑1 (MQ=255)
                                                                                                ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc    >  2:253680/1‑100 (MQ=255)
                                                                                                ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc    >  2:759566/1‑100 (MQ=255)
                                                                                                ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc    >  2:3640036/1‑100 (MQ=255)
                                                                                                 cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc    <  1:4092838/99‑1 (MQ=255)
                                                                                                 cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa   <  2:3753351/100‑1 (MQ=255)
                                                                                                 cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa   <  2:3128479/100‑1 (MQ=255)
                                                                                                 cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa   <  2:3408547/100‑1 (MQ=255)
                                                                                                 cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa   <  1:2798112/100‑1 (MQ=255)
                                                                                                  tGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAt  >  2:4235376/1‑100 (MQ=255)
                                                                                                  tGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAt  <  2:3691027/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAT  >  NC_000913/2590669‑2590866

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: