Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 401,404 G→T 8.6% intergenic (‑189/‑52) USA300HOU_RS01845 ← / → sutR hypothetical protein/HTH‑type transcriptional regulator SutR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE401,4040GT 8.6% 736.4 / 43.4 219intergenic (‑189/‑52)USA300HOU_RS01845/sutRhypothetical protein/HTH‑type transcriptional regulator SutR
Reads supporting (aligned to +/- strand):  ref base G (102/98);  new base T (11/8);  total (113/106)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 224 total aligned reads displayed.

AATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAGATGGC  >  USA300TCH1516_ALE/401314‑401494
                                                                                          |                                                                                          
aaTTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAg                                                                                            >  1:1402352/1‑91 (MQ=255)
   tATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTg                                                                                         <  2:1736033/91‑1 (MQ=255)
    aTGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGt                                                                                        <  1:53985/91‑1 (MQ=255)
    aTGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGt                                                                                        >  1:1381300/1‑91 (MQ=255)
    aTGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGt                                                                                        <  1:1513668/91‑1 (MQ=255)
      gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGTCa                                                                                      <  1:2399610/91‑1 (MQ=255)
      gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCa                                                                                      >  2:882206/1‑91 (MQ=255)
      gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCAGATTTTTTGGTGTAGATTTGACAACTATATTTGTCa                                                                                      <  2:3304172/91‑1 (MQ=255)
              ggCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACg                                                                              >  2:209306/1‑91 (MQ=255)
                cGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCa                                                                            <  1:1524121/91‑1 (MQ=255)
                  ttGAGGTGTTTGTTTCAAAATGTCATGCATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCata                                                                          <  1:23636/91‑1 (MQ=255)
                  ttGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCata                                                                          >  2:373578/1‑91 (MQ=255)
                   tGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCatat                                                                         <  1:1464048/91‑1 (MQ=255)
                     aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                       <  1:1044827/91‑1 (MQ=255)
                     aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                       >  1:931071/1‑91 (MQ=255)
                     aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATCTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                       >  2:2995579/1‑91 (MQ=255)
                        tgtttgttTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCAttt                                                                    >  2:2357227/1‑91 (MQ=255)
                          tttgtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTAc                                                                  <  1:2601733/91‑1 (MQ=255)
                            tgtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAg                                                                <  2:1925130/91‑1 (MQ=255)
                             gtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgt                                                               <  2:1226858/91‑1 (MQ=255)
                              tttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtg                                                              >  2:1289347/1‑91 (MQ=255)
                               ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtgt                                                             >  2:1442138/1‑91 (MQ=255)
                               ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtga                                                             <  2:5254/91‑2 (MQ=255)
                                tCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTa                                                            >  1:1166900/1‑91 (MQ=255)
                                 cAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTaa                                                           <  2:1267673/91‑1 (MQ=255)
                                  aaaaTGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAAt                                                          <  2:2534192/91‑1 (MQ=255)
                                   aaaTGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAGGACGCATATCATTTACAGTGTAAtt                                                         >  1:179295/1‑91 (MQ=255)
                                   aaaTGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAAtt                                                         >  1:179530/1‑91 (MQ=255)
                                      tGTCATACGTAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAg                                                      <  1:1739870/91‑1 (MQ=255)
                                      tGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAg                                                      >  1:432973/1‑91 (MQ=255)
                                      tGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAg                                                      <  1:2212985/91‑1 (MQ=255)
                                      tGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAg                                                      <  1:2169027/91‑1 (MQ=255)
                                       gTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAgg                                                     >  2:356139/1‑91 (MQ=255)
                                       gTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAgg                                                     >  2:502494/1‑91 (MQ=255)
                                        tCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGa                                                    >  2:3024219/1‑91 (MQ=255)
                                          atacataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGt                                                  <  1:3508260/91‑1 (MQ=255)
                                           tacataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAACGAGtt                                                 >  2:3139407/1‑91 (MQ=255)
                                            acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                                <  1:3576689/91‑1 (MQ=255)
                                            acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                                <  1:209306/91‑1 (MQ=255)
                                            acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                                <  1:3578217/91‑1 (MQ=255)
                                            acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                                <  2:3408749/91‑1 (MQ=255)
                                            acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGACTTg                                                >  1:3009575/1‑91 (MQ=255)
                                             cataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGa                                               >  1:2952792/1‑91 (MQ=255)
                                               taGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCCTTTACCGTGTAATTAAGGAGTTGATg                                             >  1:1287332/1‑91 (MQ=255)
                                               taGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATg                                             <  1:1982603/91‑1 (MQ=255)
                                                 ggTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCg                                           >  1:1125330/1‑91 (MQ=255)
                                                   tGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAt                                         <  1:3257578/91‑1 (MQ=255)
                                                   tGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAt                                         >  1:847500/1‑91 (MQ=255)
                                                    ggATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtg                                        >  2:2065123/1‑91 (MQ=255)
                                                    ggATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtg                                        <  2:351802/91‑1 (MQ=255)
                                                    ggATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtg                                        <  2:2152101/91‑1 (MQ=255)
                                                    ggATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtg                                        <  1:2369140/91‑1 (MQ=255)
                                                     gATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgt                                       >  1:631053/1‑91 (MQ=255)
                                                     gATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgt                                       <  1:1001214/91‑1 (MQ=255)
                                                     gATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgt                                       >  2:1963517/1‑91 (MQ=255)
                                                      aTAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgtg                                      <  1:3477255/91‑1 (MQ=255)
                                                       tAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGc                                     >  2:2787898/1‑91 (MQ=255)
                                                       ggaTCATATTTTTTGTTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGc                                     <  1:1883104/89‑1 (MQ=255)
                                                         aTCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGt                                   <  2:3587940/91‑1 (MQ=255)
                                                           cATATTTTTTGGTGTAGATTTGACAAATGTAGTTGTCAAAAAGAAGCATGTCATTTACAGTGTAATTAAGGAGTGGGTGCGATGTGCGTaa                                 <  2:1573771/91‑1 (MQ=255)
                                                           cATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTaa                                 <  2:817344/91‑1 (MQ=255)
                                                           cATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTaa                                 <  2:2738946/91‑1 (MQ=255)
                                                              atTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCg                              >  2:610145/1‑91 (MQ=255)
                                                                 ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                           <  1:3214866/91‑1 (MQ=255)
                                                                 ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                           >  1:3271654/1‑91 (MQ=255)
                                                                 ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                           <  1:74354/91‑1 (MQ=255)
                                                                 ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                           <  2:724019/91‑1 (MQ=255)
                                                                   ttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGTGTAATCGATTGa                         >  2:761810/1‑91 (MQ=255)
                                                                    tgttgtAGATTTGACAAATATAGTTGTCAAAAAGGCGCATATCAGTTACAGTGTAATTAAGGAGGTGATGCGATGTGCGTAATCGATTGaa                        <  2:2257302/88‑1 (MQ=255)
                                                                     ggTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGaaa                       >  2:1424717/1‑91 (MQ=255)
                                                                     ggTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGaaa                       >  2:285518/1‑91 (MQ=255)
                                                                     ggTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGaaa                       >  2:1428156/1‑91 (MQ=255)
                                                                      gtgtAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAg                      <  2:1972649/91‑1 (MQ=255)
                                                                        gtAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGaa                    >  1:1137068/1‑91 (MQ=255)
                                                                        gtAGATTTGACAAATATAGTTGTCAAAAAGACGCATACCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGaa                    >  1:3281563/1‑91 (MQ=255)
                                                                          aGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGCTGCGATGTGCGTAATCGATTGAAAGAAtt                  >  1:3472433/1‑91 (MQ=255)
                                                                          aGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAAtt                  >  1:2048103/1‑91 (MQ=255)
                                                                            ttttGACAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTAc                >  1:1909880/2‑91 (MQ=255)
                                                                               tGACAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAg             >  2:1311147/1‑91 (MQ=255)
                                                                               tGACAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAg             >  2:572883/1‑91 (MQ=255)
                                                                               tGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAg             >  2:3656677/1‑91 (MQ=255)
                                                                               tGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAg             >  2:165969/1‑91 (MQ=255)
                                                                                gACAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGc            >  1:3619546/1‑91 (MQ=255)
                                                                                 aCAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCa           >  1:367779/1‑91 (MQ=255)
                                                                                 aCAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCa           <  2:931071/91‑1 (MQ=255)
                                                                                  caactataTTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  1:3534044/5‑91 (MQ=255)
                                                                                  caactataTTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  1:2749989/5‑91 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  2:3267183/1‑91 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  2:3228321/1‑91 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          <  2:1792845/91‑1 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  2:668573/1‑91 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  1:460993/1‑91 (MQ=255)
                                                                                  cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc          >  2:804559/1‑91 (MQ=255)
                                                                                   aactataTTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACg         >  2:1202506/4‑91 (MQ=255)
                                                                                     ctataTTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACgag       >  1:2295448/2‑91 (MQ=255)
                                                                                     atatAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACgag       >  2:3427889/1‑91 (MQ=255)
                                                                                     atatAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACgag       >  1:87351/1‑91 (MQ=255)
                                                                                       ataGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAGAt     >  1:1680250/1‑91 (MQ=255)
                                                                                          gTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGATATGGc  <  2:175283/91‑1 (MQ=255)
                                                                                          gTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAGATGGc  >  2:3173461/1‑91 (MQ=255)
                                                                                          |                                                                                          
AATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAGATGGC  >  USA300TCH1516_ALE/401314‑401494

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: