Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,042,796 G→A 16.9% intergenic (+26/+63) 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,7960GA16.9% 332.1 / 54.8 136intergenic (+26/+63)tagE_3/catDPoly(glycerol‑phosphate) alpha‑glucosyltransferase/Putative oxidoreductase CatD
Reads supporting (aligned to +/- strand):  ref base G (46/67);  new base A (16/7);  total (62/74)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 137 total aligned reads displayed.

CATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042707‑1042886
                                                                                         |                                                                                          
cATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGt                                                                                           >  1:528094/1‑91 (MQ=255)
 aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTa                                                                                          >  1:3791866/1‑91 (MQ=255)
  tcataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         >  1:1891896/3‑91 (MQ=255)
    ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCt                                                                                       >  1:2700011/1‑91 (MQ=255)
    ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCt                                                                                       <  1:3714556/91‑1 (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:704906/1‑91 (MQ=255)
      aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAg                                                                                     <  1:2577931/91‑1 (MQ=255)
      aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAg                                                                                     <  1:1697526/91‑1 (MQ=255)
         aaaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAGGAt                                                                                  <  1:3131360/91‑1 (MQ=255)
          aaaTATTCAACGGAATCAATATTAGAAAAATTGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtt                                                                                 >  1:3434046/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:650750/1‑91 (MQ=255)
                   aCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCAt                                                                        >  1:1378262/1‑91 (MQ=255)
                     ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                                      <  1:1740561/91‑1 (MQ=255)
                     ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                                      <  1:403759/91‑1 (MQ=255)
                     ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                                      <  1:1483087/91‑1 (MQ=255)
                       aaTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATa                                                                    <  1:663873/91‑1 (MQ=255)
                        aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                                   >  1:1826283/1‑91 (MQ=255)
                        aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                                   >  1:863777/1‑91 (MQ=255)
                         tCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaaa                                                                  <  1:2890135/91‑1 (MQ=255)
                         tCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaaa                                                                  <  1:1144692/91‑1 (MQ=255)
                          cAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAt                                                                 <  1:3661973/91‑1 (MQ=255)
                            atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                               <  1:2745246/91‑1 (MQ=255)
                            atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                               <  1:2052571/91‑1 (MQ=255)
                            atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                               <  1:2802448/91‑1 (MQ=255)
                               ttAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaa                                                            <  1:1958583/91‑1 (MQ=255)
                               ttAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaa                                                            >  1:4137022/1‑91 (MQ=255)
                                tAGAAAAATGGTGAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaa                                                           >  1:3175551/1‑91 (MQ=255)
                                 aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          <  1:1699423/91‑1 (MQ=255)
                                  gAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaaa                                                         <  1:1810952/91‑1 (MQ=255)
                                       aTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagaga                                                    >  1:793145/1‑91 (MQ=255)
                                          gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 <  1:3544883/91‑1 (MQ=255)
                                          gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 >  1:143617/1‑91 (MQ=255)
                                            tAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATa                                               >  1:1016961/1‑91 (MQ=255)
                                            tAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATa                                               >  1:1019042/1‑91 (MQ=255)
                                             aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                              <  1:2277584/91‑1 (MQ=255)
                                             aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                              <  1:2278182/91‑1 (MQ=255)
                                             aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                              <  1:2278222/91‑1 (MQ=255)
                                              aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                             <  1:3260386/91‑1 (MQ=255)
                                              aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                             <  1:3072436/91‑1 (MQ=255)
                                              aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGt                                             <  1:3260805/91‑1 (MQ=255)
                                               aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  1:97942/91‑1 (MQ=255)
                                               aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  1:98024/91‑1 (MQ=255)
                                               aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  1:1625607/91‑1 (MQ=255)
                                                 cTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCCTGATAAATACAAAAAACGAGAGGATAGTCct                                          >  1:921297/1‑91 (MQ=255)
                                                   tttCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctct                                        <  1:3893780/91‑1 (MQ=255)
                                                   tttCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctct                                        <  1:3642115/91‑1 (MQ=255)
                                                   tttCAATAGCTAAATGAAAAAATCCTAGGTATTAATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctct                                        <  1:3391995/91‑1 (MQ=255)
                                                       aaTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaaa                                    <  1:4101994/91‑1 (MQ=255)
                                                       aaTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaaa                                    <  1:1915549/91‑1 (MQ=255)
                                                        aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGAAAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                                   >  1:783218/1‑91 (MQ=255)
                                                        aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGAAAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                                   >  1:783193/1‑91 (MQ=255)
                                                         tAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACa                                  >  1:1366649/1‑91 (MQ=255)
                                                          aGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACaa                                 <  1:2833262/91‑1 (MQ=255)
                                                           gCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAAt                                >  1:3171049/1‑91 (MQ=255)
                                                             tAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAg                              <  1:1125046/91‑1 (MQ=255)
                                                               aaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGtt                            >  1:3726116/1‑91 (MQ=255)
                                                                aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                           >  1:2189685/1‑91 (MQ=255)
                                                                aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                           <  1:1898093/91‑1 (MQ=255)
                                                                aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                           >  1:1396081/1‑91 (MQ=255)
                                                                  gAAAAAATCCTAGGTATTCATGAGTACATAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTtaat                         <  1:2951803/91‑1 (MQ=255)
                                                                     aaagTCCTACGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa                      <  1:3731737/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)
                                                                       aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGCTAGTCCTCTCAAACAATAGTTAATAAAAt                    >  1:1862586/1‑91 (MQ=255)
                                                                        aTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAtt                   <  1:1378226/91‑1 (MQ=255)
                                                                        aTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAtt                   <  1:2228764/91‑1 (MQ=255)
                                                                          ccTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGt                 >  1:645686/1‑91 (MQ=255)
                                                                          ccTAGGAATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGt                 <  1:1004669/91‑1 (MQ=255)
                                                                           cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGTTAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:444737/1‑91 (MQ=255)
                                                                           cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:2556188/1‑91 (MQ=255)
                                                                           cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:1472020/1‑91 (MQ=255)
                                                                           cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:1472019/1‑91 (MQ=255)
                                                                           cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:1301323/1‑91 (MQ=255)
                                                                             aGGTATTGATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              <  1:2073376/91‑1 (MQ=255)
                                                                             aGGTATTCATGAGTACCTAGGATTTTTACGCTTGATAAATACAAGAAAAGAGAGGATCGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              <  1:3313083/91‑1 (MQ=255)
                                                                             aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              <  1:726079/91‑1 (MQ=255)
                                                                             aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              <  1:1211810/91‑1 (MQ=255)
                                                                              ggTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGcc             >  1:2033092/1‑91 (MQ=255)
                                                                              ggTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGcc             >  1:3429931/1‑91 (MQ=255)
                                                                               gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  1:3252499/91‑1 (MQ=255)
                                                                               gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  1:1627255/91‑1 (MQ=255)
                                                                               gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  1:3431891/91‑1 (MQ=255)
                                                                               gTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  1:373430/91‑1 (MQ=255)
                                                                               gTATTCATGAGTACCTAAGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  1:3551866/91‑1 (MQ=255)
                                                                               gTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            >  1:360280/1‑91 (MQ=255)
                                                                                tactCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           >  1:1112964/4‑91 (MQ=255)
                                                                                tactCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           >  1:2823732/4‑91 (MQ=255)
                                                                                tactCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           >  1:1112748/4‑91 (MQ=255)
                                                                                 actCATGAATACCTAGGATTTTTACGCATGATCAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAtt          >  1:3102928/3‑91 (MQ=255)
                                                                                  ctCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  1:872796/2‑91 (MQ=255)
                                                                                     aTGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttctt      >  1:2171655/1‑91 (MQ=255)
                                                                                      tGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttcttc     <  1:2323589/91‑1 (MQ=255)
                                                                                         gTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAg  <  1:3575843/91‑1 (MQ=255)
                                                                                         |                                                                                          
CATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042707‑1042886

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: