Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 837,540:1 +T 7.7% intergenic (+34/‑235) USA300HOU_RS04100 → / → uvrB hypothetical protein/UvrABC system protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE837,5401.T 7.7% 963.3 / 71.0 246intergenic (+34/‑235)USA300HOU_RS04100/uvrBhypothetical protein/UvrABC system protein B
Reads supporting (aligned to +/- strand):  ref base . (110/117);  new base T (9/9);  total (120/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 274 total aligned reads displayed.

GGTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGCAC  >  USA300TCH1516_ALE/837450‑837630
                                                                                           |                                                                                          
ggTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttg                                                                                             <  2:417373/91‑1 (MQ=255)
ggTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttg                                                                                             >  1:664022/1‑91 (MQ=255)
 gTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttgt                                                                                           >  1:1981832/1‑91 (MQ=255)
   tATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttgttg                                                                                         <  2:2439750/91‑1 (MQ=255)
   tATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttgttg                                                                                         <  1:2083542/91‑1 (MQ=255)
   tATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATttgttg                                                                                         >  2:2180834/1‑91 (MQ=255)
    aTTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGa                                                                                        <  1:2785101/91‑1 (MQ=255)
    aTTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAACAAATttgtttga                                                                                       <  1:1342065/91‑4 (MQ=255)
      ttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGaaa                                                                                      <  2:3404032/91‑1 (MQ=255)
      ttAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGaaa                                                                                      >  2:3949594/1‑91 (MQ=255)
       tAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAAc                                                                                     <  1:2368950/91‑1 (MQ=255)
        aaTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACa                                                                                    >  2:3241390/1‑91 (MQ=255)
        aaTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAACAAATTTGTTTGAAACa                                                                                    <  1:1172317/91‑1 (MQ=39)
         aTTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAACAAATTTGTTTGAAACat                                                                                   <  1:479752/91‑2 (MQ=38)
          ttATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACaaa                                                                                  >  2:2699559/1‑91 (MQ=255)
            ataGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAat                                                                                >  2:2156461/1‑91 (MQ=255)
                 tCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAAtgatga                                                                           >  1:118973/1‑91 (MQ=255)
                 tCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAAtgatga                                                                           >  2:1029270/1‑91 (MQ=255)
                  cAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATgatgat                                                                          <  1:1025065/91‑1 (MQ=255)
                  cAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATgatgat                                                                          <  1:1811567/91‑1 (MQ=255)
                    aCTCGAATTTCGTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTa                                                                        <  1:3579889/91‑1 (MQ=255)
                     cTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTaa                                                                       >  1:221254/1‑91 (MQ=255)
                     cTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTaa                                                                       >  2:2524726/1‑91 (MQ=255)
                      tCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAat                                                                      <  1:1107263/91‑1 (MQ=255)
                       cGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAata                                                                     <  2:3263841/91‑1 (MQ=255)
                       cGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAata                                                                     <  2:3137876/91‑1 (MQ=255)
                          aTTTCATAAACATCACGAGGAAGATCGCTCAGATTAACTTTAAACGAATGTTTCAAACAAATTTGTTCAAACAAAACGTTGATTAGCATaa                                                                  >  2:3725976/1‑91 (MQ=255)
                          aTTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATaa                                                                  <  1:1760439/91‑1 (MQ=255)
                          aTTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATaa                                                                  >  1:3528030/1‑91 (MQ=255)
                           tttCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAAt                                                                 <  1:3885391/91‑1 (MQ=255)
                              cATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAAtgtg                                                              <  1:804629/91‑1 (MQ=255)
                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATgtgt                                                             >  2:2369279/1‑91 (MQ=255)
                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATgtgt                                                             <  1:3793455/91‑1 (MQ=255)
                               aTAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATgtgt                                                             >  2:3380441/1‑91 (MQ=255)
                                  aaCATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTAtt                                                          <  2:2789452/91‑1 (MQ=255)
                                   aCATCAAGATGAGCATCGCTAAGATTAACTTGAAAATAATGTTTCAAACAAATTTGTTGAAACAATAGTATGATTAATATAATGTGTAttt                                                         <  2:1892229/91‑1 (MQ=255)
                                   aCATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTAttt                                                         <  2:3696794/91‑1 (MQ=255)
                                     aTCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTAc                                                       <  1:1011730/91‑1 (MQ=255)
                                      tCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACa                                                      <  2:2492527/91‑1 (MQ=255)
                                       cAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACAt                                                     >  1:1476771/1‑91 (MQ=255)
                                         aGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATAc                                                   >  1:3615344/1‑91 (MQ=255)
                                          gATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACt                                                  <  1:3502182/91‑1 (MQ=255)
                                          gATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACt                                                  >  2:327766/1‑91 (MQ=255)
                                            tGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaa                                                >  2:3008073/1‑91 (MQ=255)
                                            tGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaa                                                <  1:1418048/91‑1 (MQ=255)
                                            tGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAAAGATGATTAATATAATGTGTATTTACATACTaa                                                <  2:211463/91‑1 (MQ=255)
                                             gAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaa                                               <  1:3340890/91‑1 (MQ=255)
                                             gAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaa                                               <  2:3607851/91‑1 (MQ=255)
                                             gAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaa                                               >  1:1519850/1‑91 (MQ=255)
                                              aaGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaaa                                              >  1:890361/1‑91 (MQ=255)
                                               aGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAATCAAAATGATGATTAATATAATGTGTATTTACATACTaaaaa                                             >  2:3472285/1‑91 (MQ=255)
                                               aGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTaaaaa                                             <  2:2013758/91‑1 (MQ=255)
                                                gATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAAAATAATGTGTATTTACATACTAAAAAt                                            <  2:2237066/91‑1 (MQ=255)
                                                  tCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATaa                                          <  1:2623491/91‑1 (MQ=255)
                                                  tCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATaa                                          <  1:66943/91‑1 (MQ=255)
                                                  tCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATaa                                          <  2:2077159/91‑1 (MQ=255)
                                                   cGCTAAGATTAACTTTCAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAAc                                         >  1:3923994/1‑91 (MQ=255)
                                                   cGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTGGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAAc                                         <  2:1301693/91‑1 (MQ=255)
                                                     cTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGCTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACaa                                       >  2:1280705/1‑91 (MQ=255)
                                                      tAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAg                                      <  2:1219944/91‑1 (MQ=255)
                                                       aaGATTATCTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGa                                     <  2:3506382/91‑1 (MQ=255)
                                                       aaGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGa                                     <  2:3662286/91‑1 (MQ=255)
                                                       aaGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGa                                     >  1:2145440/1‑91 (MQ=255)
                                                         gATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATa                                   <  1:2637523/91‑1 (MQ=255)
                                                            tAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAAc                                <  1:2681988/91‑1 (MQ=255)
                                                            tAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAAc                                <  2:2971253/91‑1 (MQ=255)
                                                            tAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAAc                                <  1:700108/91‑1 (MQ=255)
                                                             aaCTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACg                               <  2:1941566/91‑1 (MQ=255)
                                                               cTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgat                             >  2:2277583/1‑91 (MQ=255)
                                                                 ttAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgattt                           >  1:1549655/1‑91 (MQ=255)
                                                                   aaaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttga                         <  2:1463164/91‑1 (MQ=255)
                                                                   aaaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttga                         >  2:46753/1‑91 (MQ=255)
                                                                   aaaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttga                         >  1:3212658/1‑91 (MQ=255)
                                                                    aaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttgat                        >  2:2329257/1‑91 (MQ=255)
                                                                    aaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttgat                        <  2:316674/91‑1 (MQ=255)
                                                                    aaaTAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttgat                        >  2:1595338/1‑91 (MQ=255)
                                                                     aataatGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACgatttgatt                       <  1:1277618/91‑1 (MQ=255)
                                                                       taatGTTTCGTACAAATTTGTTGAAACAAGATGATGATTAATGTAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTa                     <  1:3190944/91‑1 (MQ=255)
                                                                         ttGTTTCAACAAATTTGCTTGAAACAAAATGATGATTAATACAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAg                   >  2:2667151/2‑91 (MQ=255)
                                                                         atGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAg                   >  1:1366201/1‑91 (MQ=255)
                                                                           gTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 <  1:74276/91‑1 (MQ=255)
                                                                           gTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACACACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 >  2:3900792/1‑91 (MQ=255)
                                                                           gTTTCAACAAATTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGc                 >  1:405456/1‑91 (MQ=39)
                                                                            tttCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCa                >  2:3105237/1‑91 (MQ=255)
                                                                             ttCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCaa               <  1:595369/91‑1 (MQ=255)
                                                                               tcaacaaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAg              >  2:2938384/3‑91 (MQ=255)
                                                                                aaacaaaTTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCa            <  2:1554255/91‑1 (MQ=255)
                                                                                 aacaaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCa            >  2:1777615/1‑91 (MQ=255)
                                                                                 aacaaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCa            >  1:3015019/1‑91 (MQ=255)
                                                                                 aacaaaTTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCAt           <  2:1702439/91‑1 (MQ=255)
                                                                                   caaaTTTGTTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATa          >  2:564746/1‑91 (MQ=255)
                                                                                     aaTTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGtt       >  1:2887364/1‑91 (MQ=255)
                                                                                     aaTTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGtt       >  2:905337/1‑91 (MQ=255)
                                                                                       tttGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAg     <  2:1638927/91‑1 (MQ=255)
                                                                                         tgttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGCTAAACGATTTGATTTAAGGCAAGCATAGTTAGca   >  2:1084602/1‑91 (MQ=255)
                                                                                         tgttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGca   >  1:2823461/1‑91 (MQ=255)
                                                                                         tgttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGca   >  2:2065011/1‑91 (MQ=255)
                                                                                          gttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGcac  <  1:2444175/91‑1 (MQ=255)
                                                                                          gttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGcac  >  1:3813718/1‑91 (MQ=255)
                                                                                          gttgAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGACTTAAGGCAAGCATAGTGAGcac  >  2:3197093/1‑91 (MQ=255)
                                                                                           |                                                                                          
GGTTATTTAATTATAGATCAACTCGAATTTCATAAACATCAAGATGAAGATCGCTAAGATTAACTTTAAAATAATGTTTCAAACAAATTTGTTGAAACAAAATGATGATTAATATAATGTGTATTTACATACTAAAAATAACAAGATAAACGATTTGATTTAAGGCAAGCATAGTTAGCAC  >  USA300TCH1516_ALE/837450‑837630

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