Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 455,079 A→T 6.3% intergenic (+79/‑301) tst_3 → / → speC_5 Toxic shock syndrome toxin‑1/Exotoxin type C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE455,0790AT 6.3% 703.0 / 28.1 207intergenic (+79/‑301)tst_3/speC_5Toxic shock syndrome toxin‑1/Exotoxin type C
Reads supporting (aligned to +/- strand):  ref base A (95/99);  new base T (6/7);  total (101/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 230 total aligned reads displayed.

AATTTGAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACTAATCGGAAA  >  USA300TCH1516_ALE/454989‑455169
                                                                                          |                                                                                          
 aTTTGAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTataa                                                                                           <  2:918259/91‑1 (MQ=255)
 aTTTGAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTataa                                                                                           <  1:3851901/91‑1 (MQ=255)
   ttGAAATAATCAATGATATATATAGAATGAAAGATTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTGTTataata                                                                                         >  2:3364431/1‑91 (MQ=255)
    tGAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTAtaataa                                                                                        >  1:830188/1‑91 (MQ=255)
     gAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATTATaaa                                                                                       <  1:2840573/91‑1 (MQ=255)
     gAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATaaa                                                                                       >  2:3757777/1‑91 (MQ=255)
      aaaTAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATaaaa                                                                                      <  1:155536/91‑1 (MQ=255)
       aataatCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATTATaaaaa                                                                                     <  2:1635154/91‑1 (MQ=255)
       aataatCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATaaaaa                                                                                     >  1:93322/1‑91 (MQ=255)
         taatCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATTATAAAAACa                                                                                   <  2:91469/91‑1 (MQ=255)
          aatCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATTATAAAAACac                                                                                  <  1:2036107/91‑2 (MQ=255)
          aatCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACat                                                                                  <  1:2186551/91‑1 (MQ=255)
           atCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACata                                                                                 <  1:1402271/91‑1 (MQ=255)
           atCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACata                                                                                 <  1:2912158/91‑1 (MQ=255)
           atCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACata                                                                                 <  2:2076897/91‑1 (MQ=255)
               aTGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGa                                                                             >  1:3206348/1‑91 (MQ=255)
                 gATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAAc                                                                           <  1:3634168/91‑1 (MQ=255)
                 gATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAAc                                                                           <  1:2046498/91‑1 (MQ=255)
                  atatatatAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACa                                                                          <  1:3229104/91‑1 (MQ=255)
                  atatatatAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACa                                                                          >  1:2724411/1‑91 (MQ=255)
                   tatatataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACAt                                                                         >  2:1521610/1‑91 (MQ=255)
                    atatataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACAtt                                                                        <  2:1905360/91‑1 (MQ=255)
                    atatataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACAtt                                                                        <  2:3047431/91‑1 (MQ=255)
                    atatataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCAATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACAtt                                                                        <  2:633496/91‑1 (MQ=255)
                     tatataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTg                                                                       <  1:3806196/91‑1 (MQ=255)
                      atataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGa                                                                      <  1:1274553/91‑1 (MQ=255)
                       tataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGAc                                                                     >  2:1323990/1‑91 (MQ=255)
                       tataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGAc                                                                     <  2:390705/91‑1 (MQ=255)
                        ataGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAAGATATCGAACATTGACt                                                                    >  2:2974406/1‑91 (MQ=255)
                          aGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTAc                                                                  <  1:1117980/91‑1 (MQ=255)
                             aTGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGtt                                                               <  1:1768148/91‑1 (MQ=255)
                              tGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGtta                                                              >  1:1793514/1‑91 (MQ=255)
                               gAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATaaaaa                                                                                     <  1:100280/67‑1 (MQ=255)
                               gAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATaaaaa                                                                                     >  2:100280/1‑67 (MQ=255)
                               gAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGttat                                                             >  2:3237380/1‑91 (MQ=255)
                                 aaGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGttatta                                                           >  2:356036/1‑91 (MQ=255)
                                 aaGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGttatta                                                           <  2:1675311/91‑1 (MQ=255)
                                 aaGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGttatta                                                           >  1:1747577/1‑91 (MQ=255)
                                 aaGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGttatta                                                           >  2:2480372/1‑91 (MQ=255)
                                  aGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTaa                                                          >  2:2834157/1‑91 (MQ=255)
                                  aGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTaa                                                          <  1:419169/91‑1 (MQ=255)
                                  aGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGAGTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTaa                                                          <  2:3810379/91‑1 (MQ=255)
                                    cTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgc                                                        >  2:2800920/1‑91 (MQ=255)
                                    cTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgc                                                        <  1:97423/91‑1 (MQ=255)
                                    cTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgc                                                        <  1:1344612/91‑1 (MQ=255)
                                    cTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgc                                                        <  1:681141/91‑1 (MQ=255)
                                     ttAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgct                                                       <  1:1703951/91‑1 (MQ=255)
                                     ttAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAgct                                                       <  2:3506072/91‑1 (MQ=255)
                                         gaagCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCtt                                                   <  1:3142516/91‑1 (MQ=255)
                                         gaagCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCtt                                                   >  2:1174480/1‑91 (MQ=255)
                                         gaagCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCtt                                                   >  2:1173704/1‑91 (MQ=255)
                                           agCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCtttt                                                 >  1:2251318/1‑91 (MQ=255)
                                              ggTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTg                                              >  1:788478/1‑91 (MQ=255)
                                               gTTTAATGATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGt                                             <  2:2901945/91‑1 (MQ=255)
                                               gTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGt                                             <  1:3495697/91‑1 (MQ=255)
                                               gTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGt                                             >  2:1021865/1‑91 (MQ=255)
                                                tttAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTa                                            >  1:1773150/1‑91 (MQ=255)
                                                tttAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTa                                            >  2:1770814/1‑91 (MQ=255)
                                                 ttAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTac                                           <  2:636359/91‑1 (MQ=255)
                                                 ttAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTac                                           <  2:158896/91‑1 (MQ=255)
                                                  taataaTCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTaca                                          >  2:1932133/1‑91 (MQ=255)
                                                  taataaTCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTaca                                          >  2:3662594/1‑91 (MQ=255)
                                                   aataatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTacac                                         >  1:2210699/1‑91 (MQ=255)
                                                   aataatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTacac                                         >  1:527210/1‑91 (MQ=255)
                                                     taatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACtt                                       >  1:2627311/1‑91 (MQ=255)
                                                     taatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACtt                                       >  1:398351/1‑91 (MQ=255)
                                                     taatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACtt                                       >  2:2517964/1‑91 (MQ=255)
                                                      aatCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACttt                                      <  1:1201219/91‑1 (MQ=255)
                                                        tCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTtat                                    <  1:3583907/91‑1 (MQ=255)
                                                         cccATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACAGTTtata                                   >  1:2368046/1‑91 (MQ=255)
                                                         cccATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTtata                                   <  2:1264914/91‑1 (MQ=255)
                                                            aTGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAAcc                                >  1:3353181/1‑91 (MQ=255)
                                                             tGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCa                               >  1:3281043/1‑91 (MQ=255)
                                                              gTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCaa                              <  2:598268/91‑1 (MQ=255)
                                                               tttAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAAt                             >  2:1546655/1‑91 (MQ=255)
                                                               tttAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAAt                             <  2:2746244/91‑1 (MQ=255)
                                                                ttAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATa                            <  2:1849595/91‑1 (MQ=255)
                                                                ttAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATa                            >  2:814144/1‑91 (MQ=255)
                                                                 tAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAg                           <  2:3245013/91‑1 (MQ=255)
                                                                   aTGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCt                         <  2:2563015/91‑1 (MQ=255)
                                                                   aTGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCt                         >  2:196045/1‑91 (MQ=255)
                                                                      aTTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTaa                      <  1:2568881/91‑1 (MQ=255)
                                                                      aTTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTaa                      <  2:348018/91‑1 (MQ=255)
                                                                      aTTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTaa                      <  1:2436170/91‑1 (MQ=255)
                                                                      aTTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTGATAACCAATAGCTTaa                      >  2:2982371/1‑91 (MQ=255)
                                                                        tttGATACGTGTTTTTTTATTTTATACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGa                    >  2:3868702/1‑91 (MQ=255)
                                                                          tGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGAtt                  <  1:1835006/91‑1 (MQ=255)
                                                                          tGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGAtt                  <  1:861123/91‑1 (MQ=255)
                                                                           gATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGAttt                 >  2:1574153/1‑91 (MQ=255)
                                                                           gATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGAttt                 >  1:2763757/1‑91 (MQ=255)
                                                                            aTACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTa                >  2:3523924/1‑91 (MQ=255)
                                                                             tACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTaa               <  1:1393688/91‑1 (MQ=255)
                                                                              aCGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTaaa              >  1:1055676/1‑91 (MQ=255)
                                                                              aCGTGATTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTaaa              <  2:3579342/91‑1 (MQ=255)
                                                                               cGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTaaaa             <  2:2116082/91‑1 (MQ=255)
                                                                                 tgtTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACt           <  2:2295029/91‑1 (MQ=255)
                                                                                 tgtTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACt           <  1:3757777/91‑1 (MQ=255)
                                                                                    ttttATTATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACTAAt        >  1:16028/1‑91 (MQ=255)
                                                                                    ttttATTATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACTAAt        >  1:3008200/1‑91 (MQ=255)
                                                                                    ttttATAATAAAAACATATCGAACATAGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACTAAt        >  2:2208812/1‑91 (MQ=255)
                                                                                          |                                                                                          
AATTTGAAATAATCAATGATATATATAGAATGAAAGCTTAAGAAGCGGTTTAATAATCCCATGTTTAATGATTTTGATACGTGTTTTTATAATAAAAACATATCGAACATTGACTACGTTATTAAGCTGCTTTTTTGTACACTTTATAACCAATAGCTTAAGATTTAAAACTAATCGGAAA  >  USA300TCH1516_ALE/454989‑455169

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: