Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,042,789 T→C 12.9% intergenic (+19/+70) tagE_3 → / ← catD Poly(glycerol‑phosphate) alpha‑glucosyltransferase/Putative oxidoreductase CatD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,042,7890TC12.9% 370.5 / 39.8 138intergenic (+19/+70)tagE_3/catDPoly(glycerol‑phosphate) alpha‑glucosyltransferase/Putative oxidoreductase CatD
Reads supporting (aligned to +/- strand):  ref base T (49/71);  new base C (6/12);  total (55/83)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 141 total aligned reads displayed.

AGAGAAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTT  >  USA300TCH1516_ALE/1042699‑1042879
                                                                                          |                                                                                          
 gagaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTAtt                                                                                           >  1:1846883/1‑91 (MQ=255)
  agaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTc                                                                                          <  1:2657421/91‑1 (MQ=255)
  agaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTc                                                                                          <  1:46060/91‑1 (MQ=255)
  agaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTc                                                                                          >  1:528970/1‑91 (MQ=255)
   gaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCa                                                                                         <  1:1136241/91‑1 (MQ=255)
   gaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCa                                                                                         <  1:1136354/91‑1 (MQ=255)
   gaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCa                                                                                         <  1:1136464/91‑1 (MQ=255)
   gaAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCa                                                                                         >  1:1403505/1‑91 (MQ=255)
    aaaaCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCAt                                                                                        <  1:841728/91‑1 (MQ=255)
     aaaCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATg                                                                                       >  1:2669773/1‑91 (MQ=255)
     aaaCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATg                                                                                       <  1:3371373/91‑1 (MQ=255)
      aaCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGa                                                                                      <  1:156161/91‑1 (MQ=255)
       aCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGaa                                                                                     <  1:3009638/91‑2 (MQ=255)
        cATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGt                                                                                    >  1:528298/1‑91 (MQ=255)
         aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTa                                                                                   >  1:3791866/1‑91 (MQ=255)
         aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGaata                                                                                   <  1:120574/91‑4 (MQ=255)
          tcataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                  >  1:1891896/3‑91 (MQ=255)
            ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCt                                                                                <  1:1386839/91‑1 (MQ=255)
            ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCt                                                                                <  1:489912/91‑1 (MQ=255)
            ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCt                                                                                <  1:3011657/91‑1 (MQ=255)
             taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTa                                                                               >  1:1801357/1‑91 (MQ=255)
             taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                               <  1:3814015/91‑1 (MQ=255)
              aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAg                                                                              <  1:2577931/91‑1 (MQ=255)
              aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAg                                                                              <  1:1650624/91‑1 (MQ=255)
              aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAg                                                                              <  1:1697526/91‑1 (MQ=255)
                aaaaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAGGa                                                                            <  1:322523/91‑1 (MQ=255)
                  aaaTATTCAACGGAATCAATATTAGAAAAATTGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtt                                                                          >  1:3434046/1‑91 (MQ=255)
                   aaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAttt                                                                         >  1:3151314/1‑91 (MQ=255)
                    atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtttt                                                                        >  1:3648515/1‑91 (MQ=255)
                    atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtttt                                                                        >  1:2130883/1‑91 (MQ=255)
                    atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAACACCTAGGAtttt                                                                        <  1:3181582/91‑1 (MQ=38)
                         cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACga                                                                   >  1:1230325/1‑90 (MQ=255)
                         cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                   >  1:1366268/1‑91 (MQ=255)
                         cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                   >  1:650750/1‑91 (MQ=255)
                           aCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCAt                                                                 >  1:1378262/1‑91 (MQ=255)
                             ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                               <  1:403759/91‑1 (MQ=255)
                             ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                               <  1:1740561/91‑1 (MQ=255)
                             ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                               <  1:1483087/91‑1 (MQ=255)
                               aaTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATa                                                             <  1:663873/91‑1 (MQ=255)
                                aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                            >  1:863777/1‑91 (MQ=255)
                                aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                            >  1:1826283/1‑91 (MQ=255)
                                 tCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaaa                                                           <  1:2890135/91‑1 (MQ=255)
                                  cAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAt                                                          <  1:3661973/91‑1 (MQ=255)
                                    atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                        <  1:686608/91‑1 (MQ=255)
                                    atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                        <  1:2745246/91‑1 (MQ=255)
                                    atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                        <  1:2802448/91‑1 (MQ=255)
                                      atTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaa                                                      <  1:2916596/91‑1 (MQ=255)
                                       ttAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaa                                                     <  1:1958583/91‑1 (MQ=255)
                                       ttAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaa                                                     >  1:4137022/1‑91 (MQ=255)
                                         aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                   <  1:2093829/91‑1 (MQ=255)
                                               aTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagaga                                             >  1:793145/1‑91 (MQ=255)
                                                    tAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATa                                        >  1:1016961/1‑91 (MQ=255)
                                                    tAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATa                                        >  1:1019042/1‑91 (MQ=255)
                                                     aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                       <  1:2278182/91‑1 (MQ=255)
                                                     aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                       <  1:2278146/91‑1 (MQ=255)
                                                      aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                      <  1:3260386/91‑1 (MQ=255)
                                                      aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                      <  1:3260805/91‑1 (MQ=255)
                                                      aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                      <  1:3072436/91‑1 (MQ=255)
                                                       aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                     <  1:98024/91‑1 (MQ=255)
                                                       aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                     <  1:97942/91‑1 (MQ=255)
                                                         cTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCCTGATAAATACAAAAAACGAGAGGATAGTCct                                   >  1:921297/1‑91 (MQ=255)
                                                           tttCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctct                                 <  1:3642115/91‑1 (MQ=255)
                                                           tttCAATAGCTAAATGAAAAAATCCTAGGTATTAATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctct                                 <  1:3391995/91‑1 (MQ=255)
                                                            ttCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCtctc                                <  1:726783/91‑1 (MQ=255)
                                                              cAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaa                              >  1:1973355/1‑91 (MQ=255)
                                                               aaTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaaa                             <  1:4101994/91‑1 (MQ=255)
                                                               aaTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaaa                             <  1:1915549/91‑1 (MQ=255)
                                                                aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                            <  1:2460508/91‑1 (MQ=255)
                                                                aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGAAAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                            >  1:783218/1‑91 (MQ=255)
                                                                  aGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACaa                          <  1:2833262/91‑1 (MQ=255)
                                                                   gCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAAt                         >  1:3171049/1‑91 (MQ=255)
                                                                       aaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGtt                     >  1:3726116/1‑91 (MQ=255)
                                                                        aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                    >  1:2189685/1‑91 (MQ=255)
                                                                        aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                    >  1:1396081/1‑91 (MQ=255)
                                                                        aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                    <  1:1898093/91‑1 (MQ=255)
                                                                         tGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTtaa                   >  1:357802/1‑91 (MQ=255)
                                                                             aaacTCCTACGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa               <  1:3731714/91‑1 (MQ=255)
                                                                             aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa               >  1:1120071/1‑91 (MQ=255)
                                                                             aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa               >  1:4019748/1‑91 (MQ=255)
                                                                             aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa               >  1:1825703/1‑91 (MQ=255)
                                                                             aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa               >  1:1120345/1‑91 (MQ=255)
                                                                               aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAGTACAAAAACAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAt             <  1:3302527/91‑1 (MQ=255)
                                                                               aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGCTAGTCCTCTCAAACAATAGTTAATAAAAt             >  1:1862586/1‑91 (MQ=255)
                                                                                aTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAtt            <  1:2228764/91‑1 (MQ=255)
                                                                                aTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAtt            <  1:1378226/91‑1 (MQ=255)
                                                                                  ccTAGGAATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGt          <  1:1004669/91‑1 (MQ=255)
                                                                                   cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGTTAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa         >  1:444737/1‑91 (MQ=255)
                                                                                   cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa         >  1:1439692/1‑91 (MQ=255)
                                                                                   cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa         >  1:2556188/1‑91 (MQ=255)
                                                                                   cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa         >  1:1301323/1‑91 (MQ=255)
                                                                                   cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa         >  1:1472019/1‑91 (MQ=255)
                                                                                     aGGTATTGATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc       <  1:2073376/91‑1 (MQ=255)
                                                                                     aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc       <  1:1211810/91‑1 (MQ=255)
                                                                                     aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc       <  1:726079/91‑1 (MQ=255)
                                                                                      ggTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGcc      >  1:2033092/1‑91 (MQ=255)
                                                                                       gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa     <  1:1627255/91‑1 (MQ=255)
                                                                                       gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa     <  1:373430/91‑1 (MQ=255)
                                                                                       gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa     <  1:3431891/91‑1 (MQ=255)
                                                                                       gTATTCATGAGTACCTAAGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa     <  1:3551866/91‑1 (MQ=255)
                                                                                          ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAATCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt  <  1:3509210/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGAGAAAACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTT  >  USA300TCH1516_ALE/1042699‑1042879

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: