Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000731 5,784 A→G 8.6% 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,7840AG 8.6% 1899.7 / 127.0 648intergenic (+301/‑1187)cadX/USA300HOU_pUSA300HOUMR0010cadmium efflux regulator/MarR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (305/287);  new base G (45/11);  total (350/298)
Fisher's exact test for biased strand distribution p-value = 3.26e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

Only 100 of 681 total aligned reads displayed.

GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTA  >  CP000731/5694‑5874
                                                                                          |                                                                                          
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            >  1:2026778‑M1/1‑91 (MQ=255)
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            <  1:1951015‑M1/91‑1 (MQ=255)
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCaa                                                                                            <  1:1853029‑M1/91‑1 (MQ=255)
 CGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAAt                                                                                           <  1:2874114‑M1/91‑1 (MQ=255)
 CGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTACAAt                                                                                           <  1:2768480‑M1/91‑1 (MQ=255)
 CGTATTATACTTATTGGCTTAGGAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAAt                                                                                           >  1:2242722‑M1/1‑91 (MQ=255)
    attatACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAgt                                                                                        <  2:3393846/91‑1 (MQ=255)
     ttatACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTATATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAgtg                                                                                       <  2:2993819/91‑1 (MQ=255)
       ataCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGtgtg                                                                                     <  1:3240471/91‑1 (MQ=255)
       ataCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGtgtg                                                                                     >  2:2824868/1‑91 (MQ=255)
       ataCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATAGAAAAAATCCCCACACTACTGCAATAGtgtg                                                                                     >  2:632377/1‑91 (MQ=255)
        taCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTATTGCAGTAGTGTgg                                                                                    <  2:35298/91‑1 (MQ=255)
        taCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTATTGCAGTAGTGTgg                                                                                    <  1:3301170/91‑1 (MQ=255)
         aCTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGCggg                                                                                   >  2:2232949/1‑91 (MQ=255)
          cTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTgggg                                                                                  <  1:1704630/91‑1 (MQ=255)
          cTTATTGGCGTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAAGCCCCACACTATTGCAGTAGTGTgggg                                                                                  <  1:1766072/91‑1 (MQ=38)
                ggCTTAGTAAACGCAACAAATAACTCATAGGCGTCCTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGt                                                                            <  1:1283351/91‑1 (MQ=255)
                 gCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTg                                                                           <  2:1865789/91‑1 (MQ=255)
                  cTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGc                                                                          <  2:2720142/91‑1 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         <  1:2824868/91‑1 (MQ=255)
                   ttAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCa                                                                         <  2:80999/91‑1 (MQ=255)
                    tAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAt                                                                        <  1:38902/91‑1 (MQ=255)
                    tAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAt                                                                        <  1:1367758/91‑1 (MQ=255)
                    tAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAt                                                                        >  1:1817451/1‑91 (MQ=255)
                     aGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAtt                                                                       >  2:21315/1‑91 (MQ=255)
                     aGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAtt                                                                       >  2:3310495/1‑91 (MQ=255)
                      gTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTACAATAGTGTGGGGATTGGTGCATTa                                                                      >  2:3616160/1‑91 (MQ=255)
                       tAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTaa                                                                     >  2:1388934/1‑91 (MQ=255)
                         aaCGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAAtt                                                                   >  1:111857/1‑91 (MQ=255)
                            gcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGtt                                                                <  1:693219‑M2/85‑1 (MQ=255)
                            gcaacaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGtt                                                                <  1:3743597‑M2/85‑1 (MQ=255)
                               acaAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAAt                                                             >  1:810636‑M2/4‑91 (MQ=255)
                                caAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATc                                                            <  2:2759252‑M2/89‑1 (MQ=255)
                                  aataatTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCAt                                                          >  2:1585401/1‑91 (MQ=255)
                                   ataatTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCAtt                                                         <  2:860788/91‑1 (MQ=255)
                                     aattaatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCAtttt                                                       >  2:1670475/1‑91 (MQ=255)
                                       ttaatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCg                                                     <  1:294589/91‑1 (MQ=255)
                                        taatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATGGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGt                                                    >  1:862925/1‑91 (MQ=255)
                                         aatAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGtt                                                   >  2:1331455/1‑91 (MQ=255)
                                            aGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtata                                                <  2:3827164/91‑1 (MQ=255)
                                            aGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCAGTAATTGTTAATCATTTTCGTtata                                                >  2:3595766/1‑91 (MQ=255)
                                             ggCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtatat                                               <  2:3761976/91‑1 (MQ=255)
                                             ggCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtatat                                               <  2:522661/91‑1 (MQ=255)
                                              gCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtatata                                              >  2:92397/1‑91 (MQ=255)
                                               cGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaa                                             <  1:2349474/91‑1 (MQ=255)
                                               cGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaa                                             <  1:3715796/91‑1 (MQ=255)
                                                gACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaaa                                            <  1:2549496/91‑1 (MQ=255)
                                                gACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATaaa                                            <  1:761480/91‑1 (MQ=255)
                                                 aCTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAAt                                           >  1:264832/1‑91 (MQ=255)
                                                  cTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAAtt                                          >  1:2213578/1‑91 (MQ=255)
                                                  cTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAAtt                                          <  2:2627701/91‑1 (MQ=255)
                                                  cTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAAtt                                          <  1:3194594/91‑1 (MQ=255)
                                                    ttATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtat                                        <  1:2740271/91‑1 (MQ=255)
                                                    ttATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGAGTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtat                                        <  2:3187826/91‑1 (MQ=255)
                                                     tataAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATtata                                       >  2:2472948/1‑91 (MQ=255)
                                                       taAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtat                                                 >  2:1166251/1‑79 (MQ=255)
                                                       taAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTtat                                                 <  1:1166251/79‑1 (MQ=255)
                                                       taAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAac                                     <  1:1529881/91‑1 (MQ=255)
                                                        aaGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAaca                                    <  1:3046473/91‑1 (MQ=255)
                                                         aGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAacac                                   <  1:3666377/91‑1 (MQ=255)
                                                         aGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAacac                                   >  2:1800326/1‑91 (MQ=255)
                                                         aGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAacac                                   <  2:2072401/91‑1 (MQ=255)
                                                            cAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTgg                                <  1:956589/91‑1 (MQ=255)
                                                            aaGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTgg                                <  2:2942810/90‑1 (MQ=255)
                                                              gCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGt                                                                            >  2:2324484/1‑45 (MQ=255)
                                                               cATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTta                             >  2:2881450/1‑91 (MQ=255)
                                                                aTCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtag                            <  2:2005564/91‑1 (MQ=255)
                                                                   gAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtagtag                         >  2:68527/1‑91 (MQ=255)
                                                                   gAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTtagtag                         >  2:1232738/1‑91 (MQ=255)
                                                                     aaaaaTCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATGAATTATAACACTGGTTAGTagag                       >  1:1000571/1‑91 (MQ=255)
                                                                        aaTCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGaaa                    >  1:2394184/1‑91 (MQ=255)
                                                                           ccccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGc                 <  2:3249704/91‑1 (MQ=255)
                                                                           ccccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGc                 <  1:2904453/91‑1 (MQ=255)
                                                                             ccACACTATTGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               >  2:343113/1‑91 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCCTTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               >  1:3088490/1‑91 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               <  1:2194906/91‑1 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               <  1:3592127/91‑1 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               >  2:1003842/1‑91 (MQ=255)
                                                                             ccACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaa               >  1:2449809/1‑91 (MQ=255)
                                                                              cacaCTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaaa              >  1:625734/1‑91 (MQ=255)
                                                                               acacTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCaaaa             >  1:2851308/1‑91 (MQ=255)
                                                                                cacTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAAt            <  1:572981/91‑1 (MQ=255)
                                                                                  ctactGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAACTTATAACACTGGTTAGTAGAGAAATGCAAAATGa          >  2:2533542/1‑91 (MQ=255)
                                                                                  ctactGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGa          <  2:1746648/91‑1 (MQ=255)
                                                                                   tattGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGAt         >  1:3392099/4‑91 (MQ=255)
                                                                                   tactGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGAt         <  1:3145728/91‑1 (MQ=255)
                                                                                   tactGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTACAGAAATGCAAAATGAt         >  2:3709483/1‑91 (MQ=255)
                                                                                    attGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGAtt        >  1:1619941/3‑91 (MQ=255)
                                                                                     ttGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTCGTAGAGAAATGCAAAATGATTg       >  2:3119186/2‑91 (MQ=255)
                                                                                     ttGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTg       >  1:3470297/2‑91 (MQ=255)
                                                                                     ttGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTg       >  2:430470/2‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:2124030/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:3374153/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  1:487864/1‑91 (MQ=255)
                                                                                      tGCAGTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  2:1505768/1‑91 (MQ=255)
                                                                                      tGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGa      >  2:763577/1‑91 (MQ=255)
                                                                                       gCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGaa     <  2:1447314/91‑1 (MQ=255)
                                                                                         aaTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAAtt   <  1:147621/91‑1 (MQ=255)
                                                                                         aaTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAAtt   >  1:3846387/1‑91 (MQ=255)
                                                                                          aTAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTa  <  2:1924609/91‑1 (MQ=255)
                                                                                          |                                                                                          
GCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGCATTAATTGTTAATCATTTTCGTTATATAAATTATAACACTGGTTAGTAGAGAAATGCAAAATGATTGAATTA  >  CP000731/5694‑5874

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