Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 401,399 A→C 9.0% intergenic (‑184/‑57) 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,3990AC 9.0% 671.1 / 52.4 210intergenic (‑184/‑57)USA300HOU_RS01845/sutRhypothetical protein/HTH‑type transcriptional regulator SutR
Reads supporting (aligned to +/- strand):  ref base A (101/90);  new base C (5/14);  total (106/104)
Fisher's exact test for biased strand distribution p-value = 3.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 227 total aligned reads displayed.

AAACAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAG  >  USA300TCH1516_ALE/401310‑401489
                                                                                         |                                                                                          
 aaCAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAata                                                                                          >  1:1186358/1‑91 (MQ=255)
  aCAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAatat                                                                                         <  2:888860/91‑1 (MQ=255)
  aCAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAatat                                                                                         >  2:1899627/1‑91 (MQ=255)
  aCAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAatat                                                                                         >  2:362034/1‑91 (MQ=255)
   cAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACtata                                                                                        <  1:2227033/91‑1 (MQ=255)
    aaTTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACtatat                                                                                       <  1:2209706/91‑2 (MQ=255)
    aaTTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGACATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACtatat                                                                                       <  2:3266418/91‑2 (MQ=255)
     aTTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACtatatt                                                                                      <  1:3455461/91‑3 (MQ=255)
        aTGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGt                                                                                   <  1:53985/91‑1 (MQ=255)
        aTGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGt                                                                                   <  1:1513668/91‑1 (MQ=255)
         tGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGTc                                                                                  <  2:2699757/91‑1 (MQ=255)
         tGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGTc                                                                                  <  2:2529656/91‑1 (MQ=255)
          gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTTTTTTTTTTTTGACAAATATAGTTGTCa                                                                                 >  2:882175/1‑91 (MQ=255)
          gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAACTATATTTGTCa                                                                                 <  2:3304077/91‑1 (MQ=255)
          gACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCa                                                                                 >  1:1750152/1‑91 (MQ=255)
           aCTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCaa                                                                                >  1:3569842/1‑91 (MQ=255)
               gTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAg                                                                            <  2:3468461/91‑1 (MQ=255)
                 aGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGAc                                                                          <  1:1484830/91‑1 (MQ=255)
                  ggCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTATAGATTTGACAAATATAGTTGTCAAAAAGACg                                                                         >  2:2090382/1‑91 (MQ=255)
                  ggCGTTGAGGTGGTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTAAAAAAAAAGACg                                                                         <  2:3228376/91‑1 (MQ=255)
                     gTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCat                                                                      >  1:2367976/1‑91 (MQ=255)
                        gAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATc                                                                   >  1:699209/1‑91 (MQ=255)
                         aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                  >  1:931071/1‑91 (MQ=255)
                         aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                  <  1:1044827/91‑1 (MQ=255)
                         aGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATCTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCa                                                                  >  2:2995579/1‑91 (MQ=255)
                          ggTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCAt                                                                 <  2:1965358/91‑1 (MQ=255)
                              tttgtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTAc                                                             <  1:2601758/91‑1 (MQ=255)
                                tgtttCAAAATGTCATACAAAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAg                                                           <  2:1928205/91‑1 (MQ=255)
                                 gtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgt                                                          <  1:830835/91‑1 (MQ=255)
                                 gtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgt                                                          <  2:495223/91‑1 (MQ=255)
                                 gtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgt                                                          <  2:1226858/91‑1 (MQ=255)
                                 gtttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgt                                                          <  1:1113461/91‑1 (MQ=255)
                                  tttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtg                                                         >  2:1289347/1‑91 (MQ=255)
                                   ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtgt                                                        <  1:437622/91‑1 (MQ=255)
                                   ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtgt                                                        >  2:1442138/1‑91 (MQ=255)
                                   ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtgt                                                        >  2:2132191/1‑91 (MQ=255)
                                   ttCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAgtga                                                        <  2:5254/91‑2 (MQ=255)
                                    tCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTa                                                       >  1:1397833/1‑91 (MQ=255)
                                    tCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTa                                                       >  1:1166900/1‑91 (MQ=255)
                                     cAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTaa                                                      <  2:1267673/91‑1 (MQ=255)
                                      aaaaTGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAAt                                                     <  2:595425/91‑1 (MQ=255)
                                       aaaTGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAAtt                                                    >  1:179530/1‑91 (MQ=255)
                                          tGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAg                                                 <  1:2169027/91‑1 (MQ=255)
                                           gTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAgg                                                >  2:356138/1‑91 (MQ=255)
                                            tCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGa                                               >  2:3024219/1‑91 (MQ=255)
                                             catacataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAg                                              >  2:1710864/1‑91 (MQ=255)
                                              atacataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGt                                             <  2:2312079/91‑1 (MQ=255)
                                               tacataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAACGAGtt                                            >  2:3139407/1‑91 (MQ=255)
                                                acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGCCGCATATCATTTACAGTGTAATCAAGGAGtcg                                           >  2:1749972/1‑89 (MQ=255)
                                                acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                           <  2:3408749/91‑1 (MQ=255)
                                                acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                           <  2:2814411/91‑1 (MQ=255)
                                                acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                           <  1:3576898/91‑1 (MQ=255)
                                                acataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTg                                           <  1:209306/91‑1 (MQ=255)
                                                 cataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGa                                          >  1:2952792/1‑91 (MQ=255)
                                                  ataGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGAt                                         <  2:728630/91‑1 (MQ=255)
                                                   taGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATg                                        <  1:1969625/91‑1 (MQ=255)
                                                     ggTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCg                                      <  1:2398142/91‑1 (MQ=255)
                                                     ggTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCg                                      >  1:1126356/1‑91 (MQ=255)
                                                     ggTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCg                                      >  2:2163921/1‑91 (MQ=255)
                                                       tGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAt                                    >  1:2591437/1‑91 (MQ=255)
                                                       tGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAt                                    <  1:3257578/91‑1 (MQ=255)
                                                        ggATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtg                                   >  1:1539254/1‑91 (MQ=255)
                                                         gATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgt                                  <  1:1001214/91‑1 (MQ=255)
                                                         gATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgt                                  >  1:631053/1‑91 (MQ=255)
                                                          aTAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGAtgtg                                 <  1:3477255/91‑1 (MQ=255)
                                                           tAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGc                                <  2:3656801/91‑1 (MQ=255)
                                                           tAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGc                                >  1:1314281/1‑91 (MQ=255)
                                                           ggaTCATATTTTTTGTTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGc                                <  1:1883104/89‑1 (MQ=255)
                                                             aTCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGt                              <  2:3587940/91‑1 (MQ=255)
                                                               cATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTaa                            >  1:2319689/1‑91 (MQ=255)
                                                               cATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTaa                            >  2:2762148/1‑91 (MQ=255)
                                                               cATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTaa                            >  2:321688/1‑91 (MQ=255)
                                                                 tatTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATc                          >  2:1575711/1‑91 (MQ=255)
                                                                  atTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCg                         >  2:610145/1‑91 (MQ=255)
                                                                     ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                      >  1:3271654/1‑91 (MQ=255)
                                                                     ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                      <  2:723728/91‑1 (MQ=255)
                                                                     ttttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGAtt                      <  1:3214866/91‑1 (MQ=255)
                                                                      tttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTg                     <  1:337450/91‑1 (MQ=255)
                                                                      tttAGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTg                     <  1:3133629/91‑1 (MQ=255)
                                                                       ttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGa                    >  2:1056248/1‑91 (MQ=255)
                                                                       ttGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGa                    >  1:353061/1‑91 (MQ=255)
                                                                         ggTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGaaa                  >  2:1424717/1‑91 (MQ=255)
                                                                         ggTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGaaa                  >  2:1428156/1‑91 (MQ=255)
                                                                          gtgtAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAg                 <  1:1463836/91‑1 (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)
                                                                                    gACAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGc       >  1:3619546/1‑91 (MQ=255)
                                                                                     aCAACTATATTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCa      >  1:367779/1‑91 (MQ=255)
                                                                                     aCAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCa      <  2:931071/91‑1 (MQ=255)
                                                                                     aCAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCa      >  2:2005864/1‑91 (MQ=255)
                                                                                      cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc     <  1:1256683/91‑1 (MQ=255)
                                                                                      cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc     >  2:804559/1‑91 (MQ=255)
                                                                                      cAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCAc     <  2:1792845/91‑1 (MQ=255)
                                                                                       aaaTATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACg    <  1:1716370/91‑1 (MQ=255)
                                                                                        aaTATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACga   <  2:782433/91‑1 (MQ=255)
                                                                                        aaTATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACga   >  1:3724009/1‑91 (MQ=255)
                                                                                         atatAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACgag  >  2:3427889/1‑91 (MQ=255)
                                                                                         atatAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACgag  >  1:87351/1‑91 (MQ=255)
                                                                                         |                                                                                          
AAACAATTATGACTTGTAGGCGTTGAGGTGTTTGTTTCAAAATGTCATACATAGGTGGATAATCATATTTTTTGGTGTAGATTTGACAAATATAGTTGTCAAAAAGACGCATATCATTTACAGTGTAATTAAGGAGTTGATGCGATGTGCGTAATCGATTGAAAGAATTACGAGCACGAG  >  USA300TCH1516_ALE/401310‑401489

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