Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,042,796 G→A 16.6% 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.6% 587.0 / 89.9 239intergenic (+26/+63)tagE_3/catDPoly(glycerol‑phosphate) alpha‑glucosyltransferase/Putative oxidoreductase CatD
Reads supporting (aligned to +/- strand):  ref base G (90/109);  new base A (27/13);  total (117/122)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 242 total aligned reads displayed.

ACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042706‑1042886
                                                                                          |                                                                                          
aCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAg                                                                                            <  2:421383/91‑1 (MQ=255)
  aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTa                                                                                          >  2:645199/1‑91 (MQ=255)
   ttATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         <  1:2372371/91‑1 (MQ=255)
   ttATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         >  2:248720/1‑91 (MQ=255)
    tataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAcc                                                                                        >  2:2736022/1‑91 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCt                                                                                       <  1:3884818/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTa                                                                                      >  2:1616354/1‑91 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTa                                                                                      >  1:1551984/1‑91 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATActag                                                                                      <  2:78444/91‑4 (MQ=39)
       aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAg                                                                                     <  2:2151177/91‑1 (MQ=255)
       aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAg                                                                                     <  1:615977/91‑1 (MQ=255)
       aGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAg                                                                                     <  2:615626/91‑1 (MQ=255)
        gAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAgg                                                                                    <  1:2864901/91‑1 (MQ=255)
        gAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAgg                                                                                    <  2:1168808/91‑1 (MQ=255)
         aaaaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTAGGa                                                                                   <  2:447918/91‑1 (MQ=255)
           aaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtt                                                                                 <  2:3077509/91‑1 (MQ=255)
             atatTCAGCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGTAAAAATCCTAGGTACTCATGAATACCTAGGAtttt                                                                               <  1:1950384/91‑1 (MQ=25)
             atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGTAAAAATCCTAGGTACTCATGAATACCTAGGAtttt                                                                               <  1:2104911/91‑1 (MQ=37)
             atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGTAAAAATCCTAGGTACTCATGAATACCTAGGAtttt                                                                               <  1:3441917/91‑1 (MQ=37)
              tatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAttttt                                                                              <  2:718875/91‑1 (MQ=255)
                ttCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTAc                                                                            >  2:2650335/1‑91 (MQ=255)
                  cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAttttttacg                                                                          >  2:2050887/1‑87 (MQ=255)
                  cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                          >  1:2404764/1‑91 (MQ=255)
                  cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                          >  2:1270124/1‑91 (MQ=255)
                   aaCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCa                                                                         <  1:349440/91‑1 (MQ=255)
                   aaCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCa                                                                         <  1:1605978/91‑1 (MQ=255)
                    aCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCAt                                                                        >  2:3737736/1‑91 (MQ=255)
                       gAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGCTTTTTACGCATGAt                                                                     >  1:1397446/1‑91 (MQ=255)
                         aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                                   <  2:529864/91‑1 (MQ=255)
                           cAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAt                                                                 <  2:1744995/91‑1 (MQ=255)
                            aaTATTAGAAAAATGGTTAACTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATa                                                                >  2:340490/1‑91 (MQ=255)
                            aaTATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATa                                                                >  1:1095112/1‑91 (MQ=255)
                            aaTATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATa                                                                <  1:2332036/91‑1 (MQ=255)
                              tatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACa                                                              >  1:55528/1‑91 (MQ=255)
                               atTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaa                                                             >  2:3783462/1‑91 (MQ=255)
                               atTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaa                                                             >  2:2571773/1‑91 (MQ=255)
                                ttAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaa                                                            <  2:3516320/91‑1 (MQ=255)
                                 tAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaa                                                           >  1:1132028/1‑91 (MQ=255)
                                 tAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAAGCCTAGGTATTCATGAGTACCTAGGATTTTTCCGCATGATAAATACaaaa                                                           >  2:1555206/1‑91 (MQ=255)
                                  aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          >  1:3773106/1‑91 (MQ=255)
                                  aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          >  1:2201821/1‑91 (MQ=255)
                                  aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          >  2:1926953/1‑91 (MQ=255)
                                  aGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          <  1:2445197/91‑1 (MQ=255)
                                  aGAAAAATGGTGAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaaa                                                          <  2:1875712/91‑1 (MQ=255)
                                     aaaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAag                                                       <  1:630021/91‑1 (MQ=255)
                                      aaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAaga                                                      >  2:2259433/1‑91 (MQ=255)
                                      aaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAaga                                                      >  1:1050530/1‑91 (MQ=255)
                                      aaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAaga                                                      >  2:489878/1‑91 (MQ=255)
                                       aaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagag                                                     <  1:1775674/91‑1 (MQ=255)
                                         tGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagagag                                                   <  1:711558/91‑1 (MQ=255)
                                         tGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagagag                                                   >  1:3827921/1‑91 (MQ=255)
                                         tGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagagag                                                   >  1:3665526/1‑91 (MQ=255)
                                          ggTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATCTTTACGCATGATAAATACAAAAAAAGAGAgg                                                  <  2:3659845/91‑1 (MQ=255)
                                           gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 <  2:2133673/91‑1 (MQ=255)
                                           gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 <  2:2544963/91‑1 (MQ=255)
                                           gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 <  1:1375863/91‑1 (MQ=255)
                                             tAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATa                                               >  2:428063/1‑91 (MQ=255)
                                               attCTTTTCAATTGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTATGATGTTTACGCATGATAAAAACAAAAAAAGAGAGGATAGt                                             <  2:3232914/89‑1 (MQ=255)
                                               aaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGactat                                             <  1:3841752/91‑5 (MQ=255)
                                                aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  2:1476846/91‑1 (MQ=255)
                                                       cAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaa                                     <  1:1190451/91‑1 (MQ=255)
                                                       cAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCaa                                     <  2:1847706/91‑1 (MQ=255)
                                                         aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                                   <  1:491273/91‑1 (MQ=255)
                                                         aTAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAAc                                   <  1:87303/91‑1 (MQ=255)
                                                            gCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAAt                                <  2:703643/91‑1 (MQ=255)
                                                                 gcgtaaaaaTCCTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                           >  1:973629/5‑91 (MQ=25)
                                                                      aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa                      <  2:1731130/91‑1 (MQ=255)
                                                                      aaaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaa                      >  1:2568741/1‑91 (MQ=255)
                                                                       aaaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATaaaa                     <  2:874991/91‑1 (MQ=255)
                                                                        aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAt                    >  1:2712761/1‑91 (MQ=255)
                                                                          tCCTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTg                  >  1:2551376/1‑91 (MQ=255)
                                                                          tCCTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTg                  >  1:1276861/1‑91 (MQ=255)
                                                                           ccTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGt                 <  1:1330995/91‑1 (MQ=255)
                                                                           ccTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGt                 >  2:2850747/1‑91 (MQ=255)
                                                                            cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:3814260/1‑91 (MQ=255)
                                                                            cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAAAAAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                <  2:3569209/91‑1 (MQ=255)
                                                                            cTAGGTATTCAAGAATACCCAGGATTTTTACGCACGATAAATACAAAAAAAGAGAGGATAGTCCTCTCACACAATAGATAATAAAATTGTa                >  2:141458/1‑91 (MQ=255)
                                                                            cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  2:2375409/1‑91 (MQ=255)
                                                                            cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  2:2375905/1‑91 (MQ=255)
                                                                            cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:3005897/1‑91 (MQ=255)
                                                                            cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:3182685/1‑91 (MQ=255)
                                                                            cTAGGTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:889088/1‑91 (MQ=255)
                                                                              aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              >  1:1092835/1‑91 (MQ=255)
                                                                              aGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              >  1:2983470/1‑91 (MQ=255)
                                                                               ggTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGcc             >  2:1805752/1‑91 (MQ=255)
                                                                                gTACTCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            >  1:2510103/1‑91 (MQ=255)
                                                                                 tATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           <  1:1326227/91‑1 (MQ=255)
                                                                                 tATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           <  1:1723395/91‑1 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  1:3056071/1‑91 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  1:701883/1‑91 (MQ=255)
                                                                                    tCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttc        <  1:2650335/91‑1 (MQ=255)
                                                                                     cATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttct       >  2:307113/1‑91 (MQ=255)
                                                                                     cATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttct       >  1:806944/1‑91 (MQ=255)
                                                                                       tGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttcttc     <  1:558666/91‑1 (MQ=255)
                                                                                       tGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttcttc     <  1:2953353/91‑1 (MQ=255)
                                                                                       tGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttcttc     <  1:2025956/91‑1 (MQ=255)
                                                                                        gAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGCTAATAAAATTGTAGCCATTTCTTcc    >  2:3859624/1‑91 (MQ=255)
                                                                                        aaGTGCCTTGGATTTTTACGAATTATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTcc    <  1:3022554/90‑1 (MQ=255)
                                                                                         aGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCa   >  2:2854867/1‑91 (MQ=255)
                                                                                          gTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAg  <  1:3571687/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042706‑1042886

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: