Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 168,354 C→G 5.1% intergenic (+321/+27) 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,3540CG 5.1% 383.4 / 11.3 119intergenic (+321/+27)yfkN_1/USA300TCH1516_00148Trifunctional nucleotide phosphoesterase protein YfkN/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (53/59);  new base G (4/2);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 123 total aligned reads displayed.

TGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTCCAGTATT  >  USA300TCH1516_ALE/168218‑168494
                                                                                                                                        |                                                                                                                                            
tgatgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTCAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                                                                          >  1:926710/1‑141 (MQ=255)
tgatgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAtcatca                                                                                                                                          <  2:120864/141‑1 (MQ=255)
  atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                        <  2:925671/141‑1 (MQ=255)
  atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                        <  3:119138/141‑1 (MQ=255)
  atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                        <  1:567964/141‑1 (MQ=255)
  atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                        <  3:77881/141‑1 (MQ=255)
  atgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAtt                                                                                                                                        >  1:288793/1‑141 (MQ=255)
   tgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCAttt                                                                                                                                       >  1:30334/1‑141 (MQ=255)
   tgaAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAACTGAAATGATGATCtttt                                                                                                                                       <  4:32037/141‑1 (MQ=255)
         aTTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTc                                                                                                                                 <  2:400679/141‑1 (MQ=255)
          tttttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCg                                                                                                                                <  1:107647/141‑1 (MQ=255)
           ttttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAATGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                                                               >  2:63602/1‑141 (MQ=255)
           ttttGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGa                                                                                                                               >  4:436045/1‑141 (MQ=255)
             tggAATCATCTAGCAGGTATTAATTTAAAACTCTTGAAGGATTATGTTAAGGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATGTTTAAATCTCGAGCTGAAAAGATCATCATTTCAGTTCGAgg                                                                                                                             <  2:99876/139‑1 (MQ=255)
              tGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGt                                                                                                                            >  1:165343/1‑141 (MQ=255)
              tGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGAAAATAATTTCAGTTCGAgga                                                                                                                            >  3:195695/1‑140 (MQ=255)
               gAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtt                                                                                                                           >  2:1044873/1‑141 (MQ=255)
               gAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGtt                                                                                                                           >  2:177582/1‑141 (MQ=255)
                atttatcTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                                                                          >  4:124223/5‑141 (MQ=255)
                aaTCATCTAGCAGGTATTAATTTAAAACTCTTTAATGATTATGTTAATGAAAATAGAAAATATAATATAACGAATCTATATAAATAAGCCACCAATATTTACACCTCGATCTGAAAAGGTCATCATTTCAGTTCGAGGttt                                                                                                                          <  4:301232/141‑1 (MQ=255)
                aaTCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGttt                                                                                                                          >  4:38615/1‑141 (MQ=255)
                     tcTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttat                                                                                                                     >  2:719414/1‑141 (MQ=255)
                        aGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttt                                                                                                                  <  1:128176/141‑1 (MQ=255)
                         gaaGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                                                                 >  1:528573/3‑141 (MQ=255)
                         gCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTttttatttta                                                                                                                 >  3:171593/1‑141 (MQ=255)
                                    tttAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATTTTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTa                                                                                                      <  4:343630/141‑1 (MQ=255)
                                        aaaCTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATggg                                                                                                  >  2:656653/1‑141 (MQ=255)
                                           ctctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTAt                                                                                               >  2:724032/1‑141 (MQ=255)
                                             ctTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTc                                                                                             >  2:42973/1‑141 (MQ=255)
                                              tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGATAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCt                                                                                            <  1:717584/141‑1 (MQ=255)
                                              tttAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCt                                                                                            >  3:126494/1‑141 (MQ=255)
                                                 aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGGAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCtttt                                                                                         >  2:630983/1‑141 (MQ=255)
                                                 aaGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCtttt                                                                                         >  3:307846/1‑141 (MQ=255)
                                                   ggATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTAc                                                                                       <  2:602389/141‑1 (MQ=255)
                                                    gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                                                                      <  4:59279/141‑1 (MQ=255)
                                                    gATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACt                                                                                      >  1:568218/1‑141 (MQ=255)
                                                       tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                                                                                   <  4:293860/141‑1 (MQ=255)
                                                       tATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGa                                                                                   <  4:74275/141‑1 (MQ=255)
                                                        aTGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaa                                                                                  >  3:32448/1‑141 (MQ=255)
                                                         tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGtaa                                                                                 >  2:91965/1‑141 (MQ=255)
                                                         tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                                                                 <  2:351000/141‑1 (MQ=255)
                                                         tGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGaaa                                                                                 <  1:258192/141‑1 (MQ=255)
                                                           ttAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAAtt                                                                               <  1:839600/141‑1 (MQ=255)
                                                             aaTGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTg                                                                             >  1:426253/1‑141 (MQ=255)
                                                               tGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtt                                                                           <  2:232383/141‑1 (MQ=255)
                                                                 aaaaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGtttt                                                                         <  1:452003/141‑1 (MQ=255)
                                                                   aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                                                                       >  2:961438/1‑141 (MQ=255)
                                                                   aaTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTcc                                                                       <  2:575620/141‑1 (MQ=255)
                                                                    aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGGGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                                      <  1:212820/141‑1 (MQ=255)
                                                                    aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                                      >  4:288486/1‑141 (MQ=255)
                                                                    aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                                      >  2:298552/1‑141 (MQ=255)
                                                                    aTAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCt                                                                      <  2:194176/141‑1 (MQ=255)
                                                                     tAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                                                                     <  1:823247/141‑1 (MQ=255)
                                                                     tAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTg                                                                     >  4:259032/1‑141 (MQ=255)
                                                                         gaATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGtttt                                                                 <  2:866924/141‑1 (MQ=255)
                                                                              taatataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTAt                                                            <  3:47858/141‑1 (MQ=255)
                                                                                 tataaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGt                                                         <  1:1014301/141‑1 (MQ=255)
                                                                                    aaCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTaaa                                                      <  1:1085311/141‑1 (MQ=255)
                                                                                     aCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAt                                                     <  4:32448/141‑1 (MQ=255)
                                                                                     aCGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAt                                                     >  1:389315/1‑141 (MQ=255)
                                                                                       gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt                                                   <  1:653480/141‑1 (MQ=255)
                                                                                       gAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAAttt                                                   <  1:320803/141‑1 (MQ=255)
                                                                                          tCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTtta                                                >  2:170939/1‑141 (MQ=255)
                                                                                            tatataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTATATTTttatt                                              >  2:867108/1‑141 (MQ=255)
                                                                                             atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTttatta                                             >  3:68923/1‑141 (MQ=255)
                                                                                             atataAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTAAATGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTttatta                                             >  3:7980/1‑141 (MQ=255)
                                                                                                taaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACAt                                          >  2:389898/1‑141 (MQ=255)
                                                                                                 aaataaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATa                                         <  2:529560/141‑1 (MQ=255)
                                                                                                    taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTa                                      <  4:129783/141‑1 (MQ=255)
                                                                                                    taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTa                                      <  4:15401/141‑1 (MQ=255)
                                                                                                    taaGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTa                                      <  1:153741/141‑1 (MQ=255)
                                                                                                       ggCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATct                                   >  2:938645/1‑140 (MQ=255)
                                                                                                         cacaAATATTTAAACCTCGAGCTGAAAAGATCATTATTTCAGTTCGAAGTTTTTTATTTTATACTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTa                                 <  4:79/141‑1 (MQ=255)
                                                                                                           ccAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGATTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTTCATAGTATCCTAAg                               >  2:937895/1‑141 (MQ=255)
                                                                                                            cAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGt                              <  1:406710/141‑1 (MQ=255)
                                                                                                             aaTATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTTTTTTTTTTTTTTATTACATATTATCCTAAttt                             >  4:398129/1‑141 (MQ=255)
                                                                                                                atTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTTTTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTtata                          >  2:38010/1‑141 (MQ=255)
                                                                                                                  ttAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATc                        <  4:171593/141‑1 (MQ=255)
                                                                                                                   tAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATcc                       <  1:554169/141‑1 (MQ=255)
                                                                                                                      aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGa                    <  4:126494/141‑1 (MQ=255)
                                                                                                                      aCCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGa                    <  1:724032/141‑1 (MQ=255)
                                                                                                                      aCCTCGAACTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTTTTGTATTTTTTTATTACATCGTATCCTAAGTTATATCCTgt                    >  2:287175/1‑140 (MQ=37)
                                                                                                                        cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACa                  <  1:709958/141‑1 (MQ=255)
                                                                                                                        cTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACa                  >  1:806905/1‑141 (MQ=255)
                                                                                                                         tCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAt                 <  2:410864/141‑1 (MQ=255)
                                                                                                                         tCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAt                 <  3:291894/141‑1 (MQ=255)
                                                                                                                          cgaacTGAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAACTTTTTATTACATAGTATCCTAAGTTATATCCTGACAtt                >  3:373848/5‑141 (MQ=38)
                                                                                                                          cGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACAtt                <  4:175238/141‑1 (MQ=255)
                                                                                                                             gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGt             >  2:15613/1‑141 (MQ=255)
                                                                                                                             gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGt             >  1:480142/1‑141 (MQ=255)
                                                                                                                             gCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATCTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGt             >  1:783122/1‑141 (MQ=255)
                                                                                                                              cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATGCTAAGTTATATCCTGACATTTGtt            >  4:245083/1‑141 (MQ=255)
                                                                                                                              cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGtt            >  2:456928/1‑141 (MQ=255)
                                                                                                                              cTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGtt            >  2:642502/1‑141 (MQ=255)
                                                                                                                                gAAATGATGATCTTTTCAGCTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTACTGTTTTTTATTGTACATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTtctc         >  1:972687/1‑141 (MQ=255)
                                                                                                                                  aaaGATCATCATTTCAGTTCGTGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTcc        <  1:652884/141‑1 (MQ=255)
                                                                                                                                  aaaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTATcc        >  1:1071164/1‑141 (MQ=255)
                                                                                                                                   aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTCCa       <  4:413095/141‑1 (MQ=255)
                                                                                                                                   aaGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTCCa       <  1:1080744/141‑1 (MQ=255)
                                                                                                                                        catcatTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTCCAGTAtt  <  1:900850/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
TGATGAAACATTTTTGAATCATCTAGCAGGTATTAATTTAAAACTCTTTAAGGATTATGTTAATGAAAATAGAGAATATAATATAACGAATCTATATAAATAAGGCACCAATATTTAAACCTCGAGCTGAAAAGATCATCATTTCAGTTCGAGGTTTTTTATTTTATCCTCTTTGTATGGGTATTCTTTTACTTGAAATTTGTTTTCCTGTTTTTTTATTGTAAATTTTTATTACATAGTATCCTAAGTTATATCCTGACATTTGTTCTCCAGTATT  >  USA300TCH1516_ALE/168218‑168494

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