Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,771 C→G 17.2% intergenic (+63/+36) 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,7710CG17.2% 390.2 / 68.1 154intergenic (+63/+36)acrD/ypfMaminoglycoside/multidrug efflux system/stress‑induced small enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base C (62/63);  new base G (7/19);  total (70/84)
Fisher's exact test for biased strand distribution p-value = 5.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 157 total aligned reads displayed.

GCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACA  >  NC_000913/2590672‑2590870
                                                                                                   |                                                                                                   
 cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCa                                                                                                    >  2:988859/1‑100 (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)
         cTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTc                                                                                            >  2:125519/1‑100 (MQ=255)
          ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                           <  1:2100073/100‑1 (MQ=255)
          ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCt                                                                                           <  2:4337182/100‑1 (MQ=255)
          ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCt                                                                                           <  2:4186193/100‑1 (MQ=255)
          ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGGGTACCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCtt                                                                                          <  1:2435310/100‑1 (MQ=255)
           tCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAACGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCtt                                                                                          <  1:1940745/100‑1 (MQ=255)
             ccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                        <  2:894299/100‑1 (MQ=255)
             ccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                        >  2:1310327/1‑100 (MQ=255)
              cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaat                                                                                       <  1:3990941/100‑4 (MQ=39)
               gCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                      <  1:3369016/100‑5 (MQ=38)
                cTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                     <  1:2661657/100‑6 (MQ=38)
                cTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                     <  2:1423833/100‑6 (MQ=38)
                cTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                     <  2:1377934/100‑6 (MQ=38)
                 tGAAGCTGCGCCCGGAATGAGCAATAAAAAATGCGACATGCCAATGCGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                      <  1:3573438/98‑5 (MQ=37)
                 tGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTTTCGCCTTTTTCAATTTTCCGATTAAGAACCTGCTTAGTAGGTTCttaatc                                                                                      <  1:756524/98‑5 (MQ=255)
                 tGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                    <  2:3007460/100‑1 (MQ=255)
                 tGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcgg                                                                                    <  2:3521415/100‑7 (MQ=38)
                 tGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTGCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                     <  1:2406990/99‑6 (MQ=38)
                  gAAGCCGCGCGCGGAATAAGCAATACAAAAGGCGACGTGCCAATGTGTCGCCTTGTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                     <  1:2339780/98‑6 (MQ=37)
                   aaGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTg                                                                                  <  1:2523864/100‑1 (MQ=255)
                   aaGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTg                                                                                  <  2:1888122/100‑1 (MQ=255)
                     gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtttccg                                                                                >  1:3488409/1‑97 (MQ=255)
                     gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTa                                                                                >  2:1605517/1‑100 (MQ=255)
                     gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTa                                                                                >  2:2166665/1‑100 (MQ=255)
                     gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTa                                                                                <  2:3385060/100‑1 (MQ=255)
                       cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTAc                                                                               <  1:35417/99‑1 (MQ=255)
                       cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACt                                                                              >  2:2887936/1‑100 (MQ=255)
                         cgcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGTTTCTTGCTCTTTGTACttt                                                                            >  1:614676/1‑100 (MQ=255)
                          gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTg                                                                           <  1:1175504/100‑1 (MQ=255)
                          gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTATACTTTg                                                                           >  2:146222/1‑100 (MQ=255)
                            ccGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtc                                                                         <  1:621924/100‑1 (MQ=255)
                               gAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCa                                                                      <  2:3991327/100‑1 (MQ=255)
                                 aTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAgg                                                                    <  1:3276556/100‑1 (MQ=255)
                                    aGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAt                                                                 >  1:2728171/1‑100 (MQ=255)
                                     gCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtt                                                                >  1:2008518/1‑100 (MQ=255)
                                       aaTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttac                                                              <  2:3147292/100‑1 (MQ=255)
                                        aTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttact                                                             >  2:108546/1‑100 (MQ=255)
                                            aaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTTCTGTTTGTTCTTTGTCTCAGGAATTACTTACg                                                         <  1:1171832/100‑1 (MQ=255)
                                             aaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTATTTTGTCTCAGGAATTACTTACgn                                                        <  1:2151214/100‑2 (MQ=255)
                                               aGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAg                                                      <  2:2008526/100‑1 (MQ=255)
                                                ggCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGc                                                     >  2:1147450/1‑100 (MQ=255)
                                                       tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTCCGAACCTGCCCAGCGCGTTCTTGCTGTTTGTCCCTTGGCCCAGCAATTACTTACCAAGCAtca                                                 >  1:1249581/1‑95 (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:963937/100‑1 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGCACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTATCAGGAATTACTTACGAAGCATAACTTc                                             >  1:2294808/1‑99 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                            >  2:1417886/1‑100 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                            >  2:4274848/1‑100 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                            <  2:3616522/100‑1 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                            <  1:4329798/100‑1 (MQ=255)
                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                            <  2:265794/100‑1 (MQ=255)
                                                           aaTGTGTCGCTTTTTGCAACTTTCCGATTAATCTCCTGCTTAGCGGGTTCTTGCTGTTTGTATTGCGTCTCAGGAATTACTTACGAAAAATAACTTCGa                                           <  1:3917699/99‑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                                     <  1:1883008/100‑1 (MQ=255)
                                                                  cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCt                                   >  2:389568/1‑100 (MQ=255)
                                                                   gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCtt                                  >  1:2062925/1‑100 (MQ=255)
                                                                    ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCttt                                 <  2:2728186/100‑1 (MQ=255)
                                                                        tttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGGAGCATAACTTCGATAAAGTCTTTCCAg                             >  2:2852301/1‑100 (MQ=255)
                                                                         ttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                            >  2:2791384/1‑100 (MQ=255)
                                                                           cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTc                          <  2:4314315/100‑1 (MQ=255)
                                                                           cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTc                          >  2:4012475/1‑100 (MQ=255)
                                                                           cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTATTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTc                          >  1:2874063/1‑100 (MQ=255)
                                                                            aaCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                         <  2:2493122/100‑1 (MQ=255)
                                                                            aaCGTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                         <  2:3962498/100‑1 (MQ=255)
                                                                                  ccGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGtt                   <  2:522136/100‑1 (MQ=255)
                                                                                   cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                  >  2:3239880/1‑100 (MQ=255)
                                                                                     aTTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCAc                >  2:1954795/1‑100 (MQ=255)
                                                                                      gcaaGAACCCGCTGAGCAGGTTCTGGCTGTTTGTACTTTGTCTCGGGAATTACTTGAGAAGCATAACTTCGATAAAGTCTTTACAGTTCCCCAGTTCACg               >  1:1891036/3‑100 (MQ=255)
                                                                                         aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc            >  2:231271/1‑100 (MQ=255)
                                                                                         aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc            <  2:3550053/100‑1 (MQ=255)
                                                                                         aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCACAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc            >  1:2600984/1‑100 (MQ=255)
                                                                                         aGAACCCGCGGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGGTAAAGTCTTTCCAGTTCCCCAGTTCACGTTc            >  1:3920861/1‑100 (MQ=255)
                                                                                          gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa           >  2:3059822/1‑100 (MQ=255)
                                                                                          gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa           >  1:1417053/1‑100 (MQ=255)
                                                                                          gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa           >  2:3739678/1‑100 (MQ=255)
                                                                                           aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCaa          >  1:1745020/1‑100 (MQ=255)
                                                                                           aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCaa          >  1:937631/1‑100 (MQ=255)
                                                                                            aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt         <  2:2702881/100‑1 (MQ=255)
                                                                                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt         <  1:692626/99‑1 (MQ=255)
                                                                                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc        >  2:3640036/1‑100 (MQ=255)
                                                                                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc        >  2:253680/1‑100 (MQ=255)
                                                                                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc        >  2:759566/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:3691027/100‑1 (MQ=255)
                                                                                                  tCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac   >  1:1594192/1‑100 (MQ=255)
                                                                                                   cAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaaca  <  1:2968817/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACA  >  NC_000913/2590672‑2590870

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: