Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 962,338 T→C 7.3% intergenic (+42/‑118) spsB_2 → / → addB Signal peptidase IB/ATP‑dependent helicase/deoxyribonuclease subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE962,3380TC 7.3% 355.9 / 13.5 109intergenic (+42/‑118)spsB_2/addBSignal peptidase IB/ATP‑dependent helicase/deoxyribonuclease subunit B
Reads supporting (aligned to +/- strand):  ref base T (44/57);  new base C (5/3);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 109 total aligned reads displayed.

TTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACTTTGTGTTA  >  USA300TCH1516_ALE/962249‑962424
                                                                                         |                                                                                      
tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                                       >  1:2774868/1‑91 (MQ=255)
tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                                       >  1:900638/1‑91 (MQ=255)
tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                                       >  1:2774367/1‑91 (MQ=255)
  tAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATaaa                                                                                     <  1:2662405/91‑1 (MQ=255)
   aGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAAc                                                                                    >  1:510585/1‑91 (MQ=255)
    gTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACg                                                                                   <  1:1016794/91‑1 (MQ=255)
    gTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACg                                                                                   <  1:356379/91‑1 (MQ=255)
    gTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACg                                                                                   <  1:2241377/91‑1 (MQ=255)
     tGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGa                                                                                  <  1:1849993/91‑1 (MQ=255)
     tGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGa                                                                                  <  1:1833661/91‑1 (MQ=255)
      gAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGAt                                                                                 >  1:3266146/1‑91 (MQ=255)
      gAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTATCAGTATTAAAACCGACAAACGAt                                                                                 <  1:1811360/91‑1 (MQ=21)
         tttAAACATCATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtt                                                                              <  1:3904351/91‑1 (MQ=255)
         tttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtt                                                                              <  1:899516/91‑1 (MQ=255)
         tttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtt                                                                              <  1:708920/91‑1 (MQ=255)
         tttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtt                                                                              <  1:709640/91‑1 (MQ=255)
          ttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGttt                                                                             <  1:1232550/91‑1 (MQ=255)
          ttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGttt                                                                             <  1:445163/91‑1 (MQ=255)
           tAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtttt                                                                            <  1:2884851/91‑1 (MQ=255)
           tAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtttt                                                                            <  1:2991357/91‑1 (MQ=255)
              aCATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATCTTTTAtt                                                                         >  1:2070280/1‑91 (MQ=255)
              aCATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATATTTTAtt                                                                         >  1:2070170/1‑91 (MQ=255)
               cATAATTTCAATCCTGAAAATACTATAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTAttt                                                                        <  1:3875613/91‑1 (MQ=255)
               cATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACCACATCGTTTGTCGGTTTTAATACTGATAAGCGATGTTTTAttt                                                                        >  1:4153646/1‑91 (MQ=255)
               cATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACCACATCGTTTGTCGGTTTTAATACTGATAAGCGATGTTTTAttt                                                                        >  1:4153589/1‑91 (MQ=255)
                  aaTTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGt                                                                     >  1:3661975/1‑91 (MQ=255)
                   aTTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATAAAACATCGTTTATCAGTATTAAAACCGACAAACGATGTTGTAtttgttt                                                                    <  1:2996262/91‑5 (MQ=2)
                   aTTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATAAAACATCGTTTATCAGTATTAAAACCGACAAACGATGTTGTAtttgttt                                                                    <  1:2867024/91‑5 (MQ=2)
                     ttCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAg                                                                  <  1:4059980/91‑1 (MQ=255)
                       cAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTa                                                                >  1:700272/1‑91 (MQ=255)
                        aaTCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTAc                                                               >  1:2876056/1‑91 (MQ=255)
                         aTCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTAcc                                                              <  1:578937/91‑1 (MQ=255)
                         aTCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTAcc                                                              <  1:578854/91‑1 (MQ=255)
                         aTCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTAcc                                                              <  1:1200724/91‑1 (MQ=255)
                          tCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACca                                                             >  1:426048/1‑91 (MQ=255)
                              gAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAAt                                                         <  1:3804002/91‑1 (MQ=255)
                              gAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAAt                                                         <  1:510869/91‑1 (MQ=255)
                              gAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAAt                                                         >  1:87159/1‑91 (MQ=255)
                                aaaTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTCATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaa                                                       <  1:3596343/91‑1 (MQ=255)
                                aaaTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaa                                                       <  1:1812934/91‑1 (MQ=255)
                                  aTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaaaa                                                     >  1:1267726/1‑91 (MQ=255)
                                  aTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaaaa                                                     >  1:1271017/1‑91 (MQ=255)
                                   tACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAg                                                    <  1:2310791/91‑1 (MQ=255)
                                     cTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCt                                                  <  1:3627356/91‑1 (MQ=255)
                                      tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                                 >  1:1600577/1‑91 (MQ=255)
                                      tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                                 <  1:3818663/91‑1 (MQ=255)
                                      tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                                 >  1:1600585/1‑91 (MQ=255)
                                      tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                                 <  1:466858/91‑1 (MQ=255)
                                      tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                                 >  1:1956067/1‑91 (MQ=255)
                                         aaaTTAATATGAAACAAATACAACATCGTTTGTCTGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGt                                              <  1:756104/91‑1 (MQ=255)
                                         aaaTTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGt                                              <  1:756073/91‑1 (MQ=255)
                                          aaTTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGtt                                             <  1:1660082/91‑1 (MQ=255)
                                          aaTTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGtt                                             <  1:3358280/91‑1 (MQ=255)
                                           aTTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTc                                            <  1:1140391/91‑1 (MQ=255)
                                           aTTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTc                                            <  1:1142369/91‑1 (MQ=255)
                                              aaTATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                                         >  1:1441699/1‑91 (MQ=255)
                                              aaTATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                                         >  1:1441625/1‑91 (MQ=255)
                                              aaTATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                                         <  1:1299908/91‑1 (MQ=255)
                                               atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTcattt                                        <  1:4055642/91‑5 (MQ=255)
                                               atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTcattt                                        <  1:4055556/91‑5 (MQ=255)
                                               atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaaa                                        >  1:2095850/1‑91 (MQ=255)
                                               atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaaa                                        <  1:1027944/91‑1 (MQ=255)
                                                tatGAAACACATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAAt                                       <  1:3964368/91‑1 (MQ=255)
                                                tatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAAt                                       >  1:524220/1‑91 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      >  1:943777/1‑91 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      >  1:3807276/1‑91 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      >  1:3806860/1‑91 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      >  1:3807288/1‑91 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      <  1:3571518/91‑1 (MQ=255)
                                                 atGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATg                                      >  1:3244578/1‑91 (MQ=255)
                                                  tGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGa                                     >  1:1729415/1‑91 (MQ=255)
                                                      caaaaTAAAACATCGTTTATCAGTATTAAAACCGACAAACGATGTTGTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACtt                                 >  1:697530/3‑91 (MQ=2)
                                                      aCAAATGCAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACtt                                 >  1:1013211/1‑91 (MQ=255)
                                                       cAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTta                                <  1:753074/91‑1 (MQ=255)
                                                         aaTAAAACATCGTTTATCAGTATTAAAACCGACAAACGATGTTGTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTtata                              >  1:778764/1‑91 (MQ=14)
                                                          aTACAACATCGTTTGTCGGTTTTAATTCTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTataa                             <  1:1852879/91‑1 (MQ=255)
                                                          aTACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTataa                             <  1:2110339/91‑1 (MQ=255)
                                                           tACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTataat                            >  1:580379/1‑91 (MQ=255)
                                                            acaacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTataata                           <  1:829996/91‑1 (MQ=255)
                                                            acaacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTataata                           <  1:1857396/91‑1 (MQ=255)
                                                             caacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa                          >  1:3428707/1‑91 (MQ=255)
                                                             caacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa                          >  1:3189508/1‑91 (MQ=255)
                                                             aaacaTCGTTTATCAGTATTAAAACCGACAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa                          >  1:3140967/2‑91 (MQ=18)
                                                             aaacaTCGTTTATCAGTATTAAAACCGACAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa                          >  1:3141644/2‑91 (MQ=18)
                                                              aacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATaaa                         <  1:2730822/91‑1 (MQ=255)
                                                              aacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATaaa                         >  1:1853407/1‑91 (MQ=255)
                                                                   cGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAaatcaat                    <  1:2117156/91‑1 (MQ=255)
                                                                      ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac                 <  1:113674/91‑1 (MQ=255)
                                                                      ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac                 <  1:113692/91‑1 (MQ=255)
                                                                      ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac                 <  1:3689991/91‑1 (MQ=255)
                                                                       tGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcaca                >  1:2580247/1‑91 (MQ=255)
                                                                       tGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcaca                >  1:2580117/1‑91 (MQ=255)
                                                                          cGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcacaatc             >  1:1034843/1‑91 (MQ=255)
                                                                             ttttAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACt          <  1:2994385/91‑1 (MQ=255)
                                                                             ttttAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACt          <  1:2995319/91‑1 (MQ=255)
                                                                              tttAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACtt         >  1:4049841/1‑91 (MQ=255)
                                                                              tttAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACtt         >  1:2657951/1‑91 (MQ=255)
                                                                                    aCTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACTTTGTGtt   >  1:2523386/1‑91 (MQ=255)
                                                                                    aCTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACTTTGTGtt   <  1:3892838/91‑1 (MQ=255)
                                                                                     cTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACTTTGTGTTa  <  1:4013059/91‑1 (MQ=255)
                                                                                         |                                                                                      
TTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACAATCACTTTGTGTTA  >  USA300TCH1516_ALE/962249‑962424

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: