Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,358 A→T 5.6% intergenic (+325/+23) yfkN_1 → / ← USA300TCH1516_00148 Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE168,3580AT 5.6% 441.8 / 13.2 126intergenic (+325/+23)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (57/62);  new base T (4/3);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 147 total aligned reads displayed.

AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTT  >  USA300TCH1516_ALE/168268‑168447
                                                                                          |                                                                                         
aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                           >  1:2823368/1‑91 (MQ=255)
aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                           >  1:247129/1‑91 (MQ=255)
aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                           >  1:2284914/1‑91 (MQ=255)
aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                           >  1:2111773/1‑91 (MQ=255)
 ggATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATcatcat                                                                                          >  1:2950388/1‑91 (MQ=255)
  gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                         <  1:2649266/91‑1 (MQ=255)
   aTTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                        <  1:432457/91‑1 (MQ=38)
    ttATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTc                                                                                       <  1:3798275/91‑1 (MQ=37)
       tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                                    >  1:3169752/1‑91 (MQ=255)
       tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                                    >  1:3167026/1‑91 (MQ=255)
       tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGt                                                                                    >  1:3742347/1‑91 (MQ=255)
        gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGtt                                                                                   >  1:1923203/1‑91 (MQ=255)
         ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                                  <  1:1757696/91‑1 (MQ=255)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                <  1:2685490/91‑1 (MQ=255)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                <  1:1485845/91‑1 (MQ=255)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                >  1:1512637/1‑91 (MQ=255)
           aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCTTTTCAGCTCGa                                                                                <  1:666004/91‑1 (MQ=25)
                aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                           >  1:2719580/1‑91 (MQ=255)
                aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                           >  1:3538921/1‑91 (MQ=255)
                aaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                           <  1:1391000/91‑1 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                          <  1:874976/91‑1 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                          <  1:155164/91‑1 (MQ=255)
                 aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtttt                                                                          <  1:2166680/91‑1 (MQ=255)
                  aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttttt                                                                         <  1:3953785/91‑1 (MQ=255)
                    agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTtttta                                                                       >  1:2404449/1‑91 (MQ=255)
                     gagaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                                      >  1:1496318/1‑91 (MQ=255)
                      agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                     >  1:2599597/1‑91 (MQ=255)
                      agaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                     >  1:467926/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAACAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    >  1:3532671/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    >  1:1592567/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    >  1:2526917/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    >  1:2612246/1‑91 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    <  1:2460874/91‑1 (MQ=255)
                       gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                    >  1:3615799/1‑91 (MQ=255)
                        aatataatatAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                   >  1:3005504/1‑91 (MQ=255)
                         atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                  >  1:1624878/1‑91 (MQ=255)
                         atataatataACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                  >  1:2504852/1‑91 (MQ=255)
                            taatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATcc                                                               <  1:1575291/91‑1 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                             >  1:446932/1‑91 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                             <  1:678389/91‑1 (MQ=255)
                              atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCctc                                                             >  1:2987127/1‑91 (MQ=255)
                                ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                           <  1:2824575/91‑1 (MQ=255)
                                 taaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCttt                                                          <  1:2095825/91‑1 (MQ=255)
                                   aCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGt                                                        <  1:1474307/91‑1 (MQ=255)
                                       aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                    <  1:1172964/91‑1 (MQ=255)
                                       aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                    >  1:96507/1‑91 (MQ=255)
                                       aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATgg                                                    >  1:624871/1‑91 (MQ=255)
                                         cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGt                                                  <  1:2871825/91‑1 (MQ=255)
                                         cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGt                                                  <  1:4127482/91‑1 (MQ=255)
                                         cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATAATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGt                                                  <  1:3133079/91‑1 (MQ=255)
                                          tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTa                                                 <  1:3400376/91‑1 (MQ=255)
                                           atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                                >  1:2257514/1‑91 (MQ=255)
                                            tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                               >  1:826304/1‑91 (MQ=255)
                                            tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                               >  1:222210/1‑91 (MQ=255)
                                             ataaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTc                                              >  1:2275094/1‑91 (MQ=255)
                                              caaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCt                                             <  1:4097435/90‑1 (MQ=255)
                                                  taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                         >  1:2979845/1‑91 (MQ=255)
                                                   aaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                        <  1:1050916/91‑1 (MQ=255)
                                                    aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                       <  1:1021796/91‑1 (MQ=255)
                                                        aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                                   >  1:3609822/1‑91 (MQ=255)
                                                         ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                  >  1:2899776/1‑91 (MQ=255)
                                                         ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                  <  1:1847281/91‑1 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                >  1:152503/1‑91 (MQ=255)
                                                           aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                <  1:1718931/91‑1 (MQ=255)
                                                             tatTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                              >  1:3787203/1‑91 (MQ=255)
                                                                ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGttt                           >  1:1985483/1‑91 (MQ=255)
                                                                 tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                          <  1:621451/91‑1 (MQ=255)
                                                                 tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                          >  1:4141772/1‑91 (MQ=255)
                                                                 tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                          >  1:3715657/1‑91 (MQ=255)
                                                                  aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc                         <  1:2872601/91‑1 (MQ=255)
                                                                  aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc                         <  1:2203359/91‑1 (MQ=255)
                                                                  aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTc                         <  1:3986316/91‑1 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                        >  1:2764186/1‑91 (MQ=255)
                                                                   aaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                        <  1:2398677/91‑1 (MQ=255)
                                                                    aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                       <  1:1600144/91‑1 (MQ=255)
                                                                    aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                       <  1:1375657/91‑1 (MQ=255)
                                                                     ccTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                      <  1:4098634/91‑1 (MQ=255)
                                                                      cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                     <  1:1009095/91‑1 (MQ=255)
                                                                          aGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTCTCCTGttttt                 >  1:4117191/1‑91 (MQ=255)
                                                                           gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttttt                >  1:4046648/1‑91 (MQ=255)
                                                                           gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttttt                >  1:4047049/1‑91 (MQ=255)
                                                                           gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttaaa                <  1:1228253/91‑4 (MQ=255)
                                                                              gAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt             >  1:2760187/1‑91 (MQ=37)
                                                                              gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt             <  1:2914442/91‑1 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt            <  1:2498388/91‑1 (MQ=255)
                                                                               aaaaGATCATCATTTCAGTTCGAGGTTCTTTATTTTATCCTCTTTGTATGGGTATTCTTTTCCTTGAAATTTGTTTTCCTGTTTTTTTAtc            >  1:1500283/1‑90 (MQ=255)
                                                                                aatgatgatCTTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg           >  1:2557324/4‑91 (MQ=255)
                                                                                aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg           >  1:776378/1‑91 (MQ=255)
                                                                                aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg           <  1:397438/91‑1 (MQ=255)
                                                                                   gATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaa        <  1:1806481/91‑1 (MQ=255)
                                                                                   gATCATCATTTCAGTTCGAGGTTTTTTATTTTAGCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaa        <  1:3979488/91‑1 (MQ=255)
                                                                                     tcatcaTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAt      <  1:2494152/91‑1 (MQ=255)
                                                                                     tcatcaTTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAt      <  1:2493799/91‑1 (MQ=255)
                                                                                      catcatTTCAGTTCGAGGTTTTTTATTTTATTCTCTTTGTATGGCTATTCTTTTACTAGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtt     <  1:3317168/91‑1 (MQ=255)
                                                                                      catcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtt     <  1:462958/91‑1 (MQ=255)
                                                                                      catcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtt     <  1:2059654/91‑1 (MQ=255)
                                                                                        tcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt   <  1:3026972/91‑1 (MQ=255)
                                                                                        tcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt   <  1:2235611/91‑1 (MQ=255)
                                                                                        tcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt   <  1:300184/91‑1 (MQ=255)
                                                                                         catTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttttt  <  1:2201060/91‑1 (MQ=255)
                                                                                          |                                                                                         
AGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTT  >  USA300TCH1516_ALE/168268‑168447

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