Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,350 A→T 5.5% intergenic (+317/+31) 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,3500AT 5.5% 675.7 / 23.4 198intergenic (+317/+31)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (93/94);  new base T (3/8);  total (96/102)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 200 total aligned reads displayed.

TAAAAGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTC  >  USA300TCH1516_ALE/168211‑168487
                                                                                                                                           |                                                                                                                                         
tAAAAGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAg                                                                                                                                          >  1:573425/1‑141 (MQ=255)
 aaaaGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                                                                         <  2:909786/141‑1 (MQ=255)
 aaaaGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGa                                                                                                                                         <  1:448898/141‑1 (MQ=255)
  aaaGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                                                                        >  4:369860/1‑141 (MQ=255)
  aaaGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATAGAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGat                                                                                                                                        <  1:291182/141‑1 (MQ=255)
   aaGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatc                                                                                                                                       <  1:602312/141‑1 (MQ=255)
     gTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgat                                                                                                                                     <  2:931732/141‑4 (MQ=255)
      ttGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGatcatc                                                                                                                                    <  2:1307955/141‑1 (MQ=255)
       tgatgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                                                                   >  2:148432/1‑141 (MQ=255)
         atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                 <  2:26121/141‑1 (MQ=255)
         atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                 <  4:251778/141‑1 (MQ=255)
         atgaAACATTTTTGAATCATCTAGCAGGTATTAAGTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGatgatcttt                                                                                                                                 <  2:1261576/141‑4 (MQ=255)
             aaCATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                                                                             <  2:886885/141‑1 (MQ=255)
             aaCATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAg                                                                                                                             >  2:84568/1‑141 (MQ=255)
                 tttttGAATCATCTAGCAAGTATTAAGTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                                                                         <  1:1306565/141‑1 (MQ=255)
                   tttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAg                                                                                                                       >  1:77131/1‑141 (MQ=255)
                    ttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                                                                      >  3:343645/1‑141 (MQ=255)
                    ttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                                                                      >  2:159652/1‑141 (MQ=255)
                      gAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtt                                                                                                                    <  3:146885/141‑1 (MQ=255)
                       aaTCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                                                                   <  3:275572/141‑1 (MQ=255)
                       aaTCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                                                                   >  1:872623/1‑141 (MQ=255)
                       aaTCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                                                                   <  1:539712/141‑1 (MQ=255)
                            tcTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                                                                              >  3:461297/1‑141 (MQ=255)
                             cTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatt                                                                                                             >  4:538888/1‑141 (MQ=255)
                               aGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                                                           >  2:1464123/1‑141 (MQ=255)
                                gCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                                                          >  2:95466/1‑141 (MQ=255)
                                  aGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATc                                                                                                        >  1:154694/1‑141 (MQ=255)
                                       ttAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCtt                                                                                                   <  1:736960/141‑1 (MQ=255)
                                        tAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCttt                                                                                                  <  2:254805/141‑1 (MQ=255)
                                         aaTTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                                                           >  1:921353/1‑131 (MQ=255)
                                         aaTTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTg                                                                                                 >  1:999756/1‑141 (MQ=255)
                                             tAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATg                                                                                             <  1:386135/141‑1 (MQ=255)
                                               aaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                                                                           <  3:408138/141‑1 (MQ=255)
                                               aaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCCCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                                                                           <  3:159270/141‑1 (MQ=255)
                                                   tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                                                                       >  1:166424/1‑141 (MQ=255)
                                                   tctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAtt                                                                                       >  4:648747/1‑141 (MQ=255)
                                                    ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTc                                                                                      >  2:778523/1‑141 (MQ=255)
                                                       tAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCttt                                                                                   >  3:188410/1‑141 (MQ=255)
                                                         aGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTa                                                                                 >  2:640945/1‑141 (MQ=255)
                                                           gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAACAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                                                               >  1:333515/1‑141 (MQ=255)
                                                              tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                                                                            <  2:1640761/141‑1 (MQ=255)
                                                                tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                                                          <  4:668315/141‑1 (MQ=255)
                                                                 gTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAt                                                                         >  1:956172/1‑141 (MQ=255)
                                                                    aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                                                                      <  1:110701/141‑1 (MQ=255)
                                                                      tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtt                                                                    <  1:1283238/141‑1 (MQ=255)
                                                                        aaaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                                                                  >  4:16966/1‑141 (MQ=255)
                                                                          aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                                                                <  1:190621/141‑1 (MQ=255)
                                                                           aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                               <  1:1327598/141‑1 (MQ=255)
                                                                           aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                               <  4:107435/141‑1 (MQ=255)
                                                                            tAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                                                              >  3:308067/1‑141 (MQ=255)
                                                                            tAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                                                              <  2:329030/141‑1 (MQ=255)
                                                                             agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                                                             <  1:562386/141‑1 (MQ=255)
                                                                             agagAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGt                                                             >  1:963141/1‑141 (MQ=255)
                                                                                      aatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAtt                                                    >  4:387075/1‑141 (MQ=255)
                                                                                       atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg                                                   <  2:517568/141‑1 (MQ=255)
                                                                                       atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg                                                   <  1:1023483/141‑1 (MQ=255)
                                                                                       atataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTg                                                   <  3:700609/141‑1 (MQ=255)
                                                                                        tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt                                                  <  4:451627/141‑1 (MQ=255)
                                                                                        tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt                                                  >  2:1568029/1‑141 (MQ=255)
                                                                                         ataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTa                                                 <  1:658187/141‑1 (MQ=255)
                                                                                          taaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                 <  2:1131262/60‑1 (MQ=255)
                                                                                          taaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                 >  1:1131262/1‑60 (MQ=255)
                                                                                              gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt                                            <  3:256732/141‑1 (MQ=255)
                                                                                              gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt                                            <  2:161669/141‑1 (MQ=255)
                                                                                               aaTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt                                           >  2:609094/1‑141 (MQ=255)
                                                                                               aaTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAtttt                                           >  2:1492613/1‑141 (MQ=255)
                                                                                                aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttattt                                                                                                            <  3:375953/75‑1 (MQ=255)
                                                                                                aTCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttttt                                          <  1:551419/141‑1 (MQ=255)
                                                                                                  cTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTttat                                        <  1:1055556/141‑1 (MQ=255)
                                                                                                   tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTttatt                                       >  2:1146849/1‑141 (MQ=255)
                                                                                                    atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTttatta                                      >  1:819866/1‑141 (MQ=255)
                                                                                                     tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt                                                  <  2:215518/128‑1 (MQ=255)
                                                                                                     tataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTAc                                     <  2:1085099/141‑1 (MQ=255)
                                                                                                       taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACAt                                   <  4:558703/141‑1 (MQ=255)
                                                                                                       taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACAt                                   >  4:96277/1‑141 (MQ=255)
                                                                                                          ataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGt                                <  2:1484620/141‑1 (MQ=255)
                                                                                                             aGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATc                             >  1:1389737/1‑141 (MQ=255)
                                                                                                              ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATcc                            >  3:497619/1‑141 (MQ=255)
                                                                                                                 aCCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTaa                         >  1:1228959/1‑141 (MQ=255)
                                                                                                                     atatTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTta                     >  1:498222/1‑141 (MQ=255)
                                                                                                                       atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCATCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATGCTAAGTtcta                   >  4:133954/1‑141 (MQ=255)
                                                                                                                        tttAAACATCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTatat                  >  3:561522/1‑141 (MQ=255)
                                                                                                                        aaaaaaCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTatat                  >  1:1078276/4‑141 (MQ=38)
                                                                                                                         ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATc                 <  1:778523/141‑1 (MQ=255)
                                                                                                                         ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATc                 <  2:572793/141‑1 (MQ=255)
                                                                                                                         ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATc                 >  1:1479117/1‑141 (MQ=255)
                                                                                                                         ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATc                 <  2:1417951/141‑1 (MQ=255)
                                                                                                                          tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATcc                >  4:343167/1‑141 (MQ=255)
                                                                                                                          tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATcc                >  2:1419928/1‑141 (MQ=255)
                                                                                                                           aaaCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCt               >  2:852087/1‑141 (MQ=255)
                                                                                                                             aCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGa             <  1:1257331/141‑1 (MQ=38)
                                                                                                                               cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACa           <  1:1038236/141‑1 (MQ=255)
                                                                                                                               cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACa           <  3:648747/141‑1 (MQ=255)
                                                                                                                                tCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAt          <  3:446180/141‑1 (MQ=255)
                                                                                                                                 cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAtt         <  2:1187357/141‑1 (MQ=255)
                                                                                                                                 cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAtt         >  1:1314887/1‑141 (MQ=255)
                                                                                                                                    gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGt      >  2:714954/1‑141 (MQ=255)
                                                                                                                                     cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGtt     <  2:164987/141‑1 (MQ=255)
                                                                                                                                       gAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTtct   <  2:499596/141‑1 (MQ=255)
                                                                                                                                        aaaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTtctc  <  2:1554713/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TAAAAGTTGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTC  >  USA300TCH1516_ALE/168211‑168487

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: