Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 962,321 G→A 6.9% intergenic (+25/‑135) 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,3210GA 6.9% 387.7 / 12.6 116intergenic (+25/‑135)spsB_2/addBSignal peptidase IB/ATP‑dependent helicase/deoxyribonuclease subunit B
Reads supporting (aligned to +/- strand):  ref base G (47/61);  new base A (4/4);  total (51/65)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 124 total aligned reads displayed.

TCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACA  >  USA300TCH1516_ALE/962231‑962410
                                                                                          |                                                                                         
tCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTACAAATTAATATGAAACAAATACAACATCGTTTg                                                                                           >  1:3544521/1‑91 (MQ=255)
tCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTg                                                                                           >  1:2639784/1‑91 (MQ=255)
tCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTg                                                                                           <  1:1324339/91‑1 (MQ=255)
tCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTg                                                                                           <  1:1324221/91‑1 (MQ=255)
tCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTg                                                                                           <  1:2514143/91‑1 (MQ=255)
 cATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGt                                                                                          >  1:3940233/1‑91 (MQ=255)
  aTTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTc                                                                                         >  1:919763/1‑91 (MQ=255)
    ttAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCgg                                                                                       <  1:2280481/91‑1 (MQ=255)
      aGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGtt                                                                                     <  1:941416/91‑1 (MQ=255)
           cTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAAt                                                                                <  1:2591020/91‑1 (MQ=255)
           cTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAAt                                                                                >  1:200037/1‑91 (MQ=255)
           cTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAAt                                                                                <  1:1056523/91‑1 (MQ=255)
           cTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAAt                                                                                <  1:3806064/91‑1 (MQ=255)
              gCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACt                                                                             <  1:2637128/91‑1 (MQ=255)
              gCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACt                                                                             >  1:3876438/1‑91 (MQ=255)
              gCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACt                                                                             >  1:3876420/1‑91 (MQ=255)
               ccATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTg                                                                            >  1:101780/1‑91 (MQ=255)
               ccATTTAGTGAATTTAAACATAATGTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATAAAACATCGTTTATCAGTATTAAAAccg                                                                            <  1:1582668/91‑3 (MQ=25)
                  tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                         >  1:2774868/1‑91 (MQ=255)
                  tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                         >  1:2774367/1‑91 (MQ=255)
                  tttAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATa                                                                         >  1:900638/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)
                       tGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGa                                                                    <  1:1833661/91‑1 (MQ=255)
                       tGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGa                                                                    <  1:1849993/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:708920/91‑1 (MQ=255)
                           tttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGtt                                                                <  1:709640/91‑1 (MQ=255)
                            ttAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGttt                                                               <  1:1232550/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)
                                 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)
                                       ttCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAg                                                    <  1:4059980/91‑1 (MQ=255)
                                         cAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTa                                                  >  1:701139/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)
                                            tCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACca                                               >  1:425641/1‑91 (MQ=255)
                                            tCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACca                                               >  1:426048/1‑91 (MQ=255)
                                                gAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAAt                                           >  1:87159/1‑91 (MQ=255)
                                                gAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAAt                                           <  1:510869/91‑1 (MQ=255)
                                                  aaaTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTCATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaa                                         <  1:3596343/91‑1 (MQ=255)
                                                  aaaTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaa                                         <  1:1812934/91‑1 (MQ=255)
                                                    aTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaaaa                                       >  1:1271017/1‑91 (MQ=255)
                                                    aTACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATaaaa                                       >  1:1267726/1‑91 (MQ=255)
                                                     tACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAg                                      <  1:2310791/91‑1 (MQ=255)
                                                       cTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCt                                    <  1:3627356/91‑1 (MQ=255)
                                                        tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                   >  1:1956067/1‑91 (MQ=255)
                                                        tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                   <  1:466858/91‑1 (MQ=255)
                                                        tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                   >  1:1600585/1‑91 (MQ=255)
                                                        tAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTa                                   >  1:1600577/1‑91 (MQ=255)
                                                        tAAAAATTAATATGAAACAAATACAACATCGTTTGTCCGTTTTAATACTGACAAACGATGTTTTAGTTTGTTAGTACCACAATAACAGCTa                                   >  1:3861976/1‑91 (MQ=255)
                                                           aaaTTAATATGAAACAAATACAACATCGTTTGTCTGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGt                                <  1:756104/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)
                                                                aaTATGAATCAAATACCACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                           <  1:1299492/91‑1 (MQ=255)
                                                                aaTATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                           >  1:1441699/1‑91 (MQ=255)
                                                                aaTATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaa                           >  1:1441625/1‑91 (MQ=255)
                                                                 atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTcattt                          <  1:4055556/91‑5 (MQ=255)
                                                                 atatGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGaaa                          <  1:2507486/91‑1 (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:3244578/1‑91 (MQ=255)
                                                                      aaaCAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAAc                     <  1:34937/91‑1 (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:3189508/1‑91 (MQ=255)
                                                                               caacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa            >  1:3428707/1‑91 (MQ=255)
                                                                               aaacaTCGTTTATCAGTATTAAAACCGACAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTAtaataa            >  1:3141644/2‑91 (MQ=18)
                                                                                aacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATaaa           >  1:1853407/1‑91 (MQ=255)
                                                                                aacaTCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATaaa           <  1:2730822/91‑1 (MQ=255)
                                                                                     cGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAaatcaat      <  1:2117156/91‑1 (MQ=255)
                                                                                        ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac   <  1:113692/91‑1 (MQ=255)
                                                                                        ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac   <  1:113674/91‑1 (MQ=255)
                                                                                        ttGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcac   <  1:3689991/91‑1 (MQ=255)
                                                                                         tGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATcaatcaca  >  1:2580247/1‑91 (MQ=255)
                                                                                          |                                                                                         
TCATTTAGGTTCTGGCCATTTAGTGAATTTAAACATAATTTCAATCCTGAAAATACTAAAAATTAATATGAAACAAATACAACATCGTTTGTCGGTTTTAATACTGATAAACGATGTTTTATTTTGTTAGTACCACAATAAAAGCTAAGTTCGAAATGAACTTATAATAAATCAATCACA  >  USA300TCH1516_ALE/962231‑962410

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: