Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 5,784 A→G 10.7% intergenic (+301/‑1187) cadX → / → USA300HOU_pUSA300HOUMR0010 cadmium efflux regulator/MarR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007315,7840AG10.7% 2474.5 / 244.4 899intergenic (+301/‑1187)cadX/USA300HOU_pUSA300HOUMR0010cadmium efflux regulator/MarR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (457/340);  new base G (90/6);  total (552/347)
Fisher's exact test for biased strand distribution p-value = 2.35e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-02

Only 100 of 935 total aligned reads displayed.

GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTA  >  CP000731/5694‑5874
                                                                                          |                                                                                          
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            >  1:2081288‑M1/1‑91 (MQ=255)
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            >  1:162397‑M1/1‑91 (MQ=255)
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAAGAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            >  2:3111117‑M1/1‑91 (MQ=255)
 CGTATTATACTTATTGGCTTAGTAAACGCAACAATTAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATACACACACTACTGCAAt                                                                                           <  2:821236‑M1/91‑1 (MQ=255)
 CGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAAt                                                                                           <  1:2948775‑M1/91‑1 (MQ=255)
 CGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATAGAAAAAATCCCAACACTACTGCAAt                                                                                           <  2:3839784‑M1/91‑1 (MQ=255)
    attatACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAgt                                                                                        >  1:1874015/1‑91 (MQ=255)
    attatACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAgt                                                                                        >  1:1713776/1‑91 (MQ=255)
     ttatACTTATTGGCTTAGTAAACGCAACAAATAGTTACTAGGCGACATTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAgtg                                                                                       >  1:3159226/1‑91 (MQ=255)
       ataCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGtgtg                                                                                     >  1:76737/1‑91 (MQ=255)
        taCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTATTGCAGTAGTGTgg                                                                                    <  1:227337/91‑1 (MQ=255)
         aCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCTACTTTATAAGCCAGCATCGAAACAATCCCCACACTACTGCAATAGTGTggg                                                                                   >  2:2373922/1‑91 (MQ=255)
         aCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTggg                                                                                   <  2:2277783/91‑1 (MQ=255)
         aCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTggg                                                                                   >  2:162176/1‑91 (MQ=255)
          cTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTgggg                                                                                  >  1:512033/1‑91 (MQ=255)
            tATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGAt                                                                                <  1:3592169/91‑1 (MQ=255)
             aTTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACCCTACTGCAATAGTGTGGGGAtt                                                                               >  1:442072/1‑91 (MQ=255)
             aTTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGAtt                                                                               >  2:2955045/1‑91 (MQ=255)
               tGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTgg                                                                             <  1:632955/91‑1 (MQ=255)
                ggCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGt                                                                            >  1:2785768/1‑91 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAATCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           <  2:3477051/91‑1 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           >  1:2490435/1‑91 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           <  1:3286240/91‑1 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           >  2:3600800/1‑91 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           >  1:3936997/1‑91 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGAATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           >  2:3546323/1‑91 (MQ=255)
                  cTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGc                                                                          <  2:3687323/91‑1 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         <  1:2581728/91‑1 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         >  1:848581/1‑91 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         <  1:930870/91‑1 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         <  2:1124102/91‑1 (MQ=255)
                    tAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAt                                                                        >  2:60883/1‑91 (MQ=255)
                        aaaCGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAAt                                                                    <  1:379480/91‑1 (MQ=255)
                          acgcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTg                                                                  >  1:311741‑M2/9‑91 (MQ=255)
                           cgcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGt                                                                 <  1:1420092‑M2/84‑1 (MQ=255)
                            gcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGtt                                                                >  1:3048428‑M2/7‑91 (MQ=255)
                            gcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGtt                                                                <  2:76480‑M2/85‑1 (MQ=255)
                             caacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTa                                                               <  1:2112688‑M2/86‑1 (MQ=255)
                             caacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTa                                                               <  1:1325360‑M2/86‑1 (MQ=255)
                              aacaAATAAGTAAGAGGCGACTTTATAAGCCAGCAGCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTaa                                                              <  2:171633‑M2/87‑1 (MQ=255)
                               acaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAAt                                                             >  1:195826‑M2/4‑91 (MQ=255)
                               acaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGATGCATTAATTGTTAAt                                                             >  2:3032309‑M2/4‑91 (MQ=255)
                               acaAATAACTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAACAGAGCGGGGATTGGTGCATCAATCGTtaac                                                             >  2:181276‑M2/4‑90 (MQ=255)
                                caAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATc                                                            >  1:2712524‑M2/3‑91 (MQ=255)
                                caAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATc                                                            >  1:2376203‑M2/3‑91 (MQ=255)
                                  aataatTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATCGTGTGGGGATTGGTGCATTAATTGTTAATCAt                                                          >  2:2449085/1‑91 (MQ=255)
                                  aataatTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCAt                                                          >  1:3338539/1‑91 (MQ=255)
                                       ttaatAGGCGACTTTATAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGt                                                    >  2:3732638/1‑91 (MQ=255)
                                         aatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGtt                                                   >  1:118911/1‑91 (MQ=255)
                                         aatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGtt                                                   >  1:1368742/1‑91 (MQ=255)
                                          atAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTta                                                  >  1:895750/1‑91 (MQ=255)
                                          atAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGCGGGGATTGGTGCATTAATTGTTAATCATTTTCGTta                                                  >  2:698945/1‑91 (MQ=255)
                                             ggCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtatat                                               <  1:712960/91‑1 (MQ=255)
                                             ggCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtatat                                               <  2:2614693/91‑1 (MQ=255)
                                               cGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaa                                             <  1:60292/91‑1 (MQ=255)
                                                gACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaaa                                            <  2:1023466/91‑1 (MQ=255)
                                                gACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaaa                                            >  1:2659980/1‑91 (MQ=255)
                                                    ttATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtat                                        <  2:1960348/91‑1 (MQ=255)
                                                    ttATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtat                                        <  2:1876149/91‑1 (MQ=255)
                                                    ttATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtat                                        <  2:153113/91‑1 (MQ=255)
                                                     tataAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtata                                       >  1:1128192/1‑91 (MQ=255)
                                                     tataAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtata                                       >  1:2877537/1‑91 (MQ=255)
                                                      ataAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATaa                                      >  1:1924453/1‑91 (MQ=255)
                                                      ataAGCCAGCATCGAAAAAAGCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATaa                                      <  1:1253552/91‑1 (MQ=255)
                                                       taAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAac                                     >  2:1441397/1‑91 (MQ=255)
                                                        aaGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAaca                                    >  1:3231570/1‑91 (MQ=255)
                                                        aaGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAaca                                    >  1:2919158/1‑91 (MQ=255)
                                                        aaGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAAGTGTTAATCATTTTCGTTATATAAATTATAaca                                    <  2:2852510/91‑1 (MQ=255)
                                                         aGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAacac                                   <  1:200726/91‑1 (MQ=255)
                                                          gCCAGCATCGAAAAAATCCCCACACTACTGCACTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACt                                  <  2:89054/91‑1 (MQ=255)
                                                            cAGCATCGAAAAAATCGCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTgg                                >  1:3894407/1‑91 (MQ=255)
                                                            aaGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTgg                                <  2:490060/90‑1 (MQ=255)
                                                             aGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGt                               <  1:1664637/91‑1 (MQ=255)
                                                                aTCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATTATTTTCGTTATATCAATTATAACACTGGttcg                            >  2:389923/1‑89 (MQ=255)
                                                                 tCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtagt                           >  2:20984/1‑91 (MQ=255)
                                                                  cGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtagta                          <  1:2544754/91‑1 (MQ=255)
                                                                  cGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtagta                          <  2:319690/91‑1 (MQ=255)
                                                                    aaaaaaTCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTaga                        >  2:2678090/1‑91 (MQ=255)
                                                                          tCCCCACACTACTGCAATAGTGTGGGGATTGGTACATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATg                  <  1:3505056/91‑1 (MQ=255)
                                                                            cccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCa                >  1:2193828/1‑91 (MQ=255)
                                                                            cccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCa                >  1:3444115/1‑91 (MQ=255)
                                                                            cccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCa                >  2:3271391/1‑91 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               >  2:2549494/1‑91 (MQ=255)
                                                                              cacaCTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaaa              >  1:1772156/1‑91 (MQ=255)
                                                                               acacTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaaaa             <  1:1489611/91‑1 (MQ=255)
                                                                                cacTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAAt            >  2:2059965/1‑91 (MQ=255)
                                                                                  ctattGCAGTAGTGTGGGGATTGGTGCATTAATTGTCAATTATTTTCGTTATATAAATTATAACACTGGTTAGTCGAGAAATGCAAAATGa          >  2:1284171/1‑91 (MQ=255)
                                                                                   tactGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGAt         <  2:3352733/91‑1 (MQ=255)
                                                                                    actGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGAtt        >  1:1093958/1‑91 (MQ=255)
                                                                                     ttGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTg       >  2:2757801/2‑91 (MQ=255)
                                                                                     ctGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTg       >  1:3724991/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCCTTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGCAGAGAAATGCAAAATGATTGa      >  2:547065/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:1550838/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  2:2175457/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:3656082/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:384378/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:409262/1‑91 (MQ=255)
                                                                                       gCAATCGTGTGGGGATTGGAGCATTAATTGTTAACCATTTTCGCCATATAACTTATAACACTGGTTAGTAGAGAAATGCACAACGATTGaa     >  2:1387719/1‑91 (MQ=255)
                                                                                       gCAATAGTGTGGGGATTGGTGCATAAATTGTTAATCATGTTCGATATTTAGATTATAACACCGGCTAGTAGAGAAATGCAACATGACTGaa     >  2:1212922/1‑91 (MQ=255)
                                                                                          aTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTa  <  1:101288/91‑1 (MQ=255)
                                                                                          |                                                                                          
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTA  >  CP000731/5694‑5874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: