Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,042,796 G→A 18.8% 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,7960GA18.8% 631.4 / 126.8 286intergenic (+26/+63)tagE_3/catDPoly(glycerol‑phosphate) alpha‑glucosyltransferase/Putative oxidoreductase CatD
Reads supporting (aligned to +/- strand):  ref base G (99/132);  new base A (31/23);  total (131/155)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 292 total aligned reads displayed.

ACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042706‑1042886
                                                                                          |                                                                                          
aCATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAg                                                                                            >  2:2233920/1‑91 (MQ=255)
 cATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGt                                                                                           <  1:3670513/91‑1 (MQ=255)
 cATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGt                                                                                           >  1:209388/1‑91 (MQ=255)
  aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTa                                                                                          >  1:2343866/1‑91 (MQ=255)
  aTTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTa                                                                                          >  1:1692625/1‑91 (MQ=255)
   ttATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         >  1:1085773/1‑91 (MQ=255)
   ttATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         <  1:1241876/91‑1 (MQ=255)
   ttATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAc                                                                                         <  1:2244910/91‑1 (MQ=255)
    tataGCAAAATATTCAACGGAATCAAGATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAcc                                                                                        <  2:3641861/91‑1 (MQ=255)
    tataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAcc                                                                                        >  1:3775027/1‑91 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCt                                                                                       <  1:899247/91‑1 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCt                                                                                       <  1:2493481/91‑1 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATAAGTACCt                                                                                       >  1:1958092/1‑91 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCt                                                                                       <  2:2659290/91‑1 (MQ=255)
     ataGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCt                                                                                       <  2:1151831/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTa                                                                                      >  1:3380789/1‑91 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                                      <  2:2126268/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                                      <  1:2217538/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                                      <  1:554261/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                                      <  2:2832277/91‑1 (MQ=255)
      taGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTACTCATGAATACCTa                                                                                      <  1:857336/91‑1 (MQ=255)
         aaaaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGa                                                                                   <  2:228513/91‑1 (MQ=255)
           aaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtt                                                                                 <  2:1550007/91‑1 (MQ=255)
           aaaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAtt                                                                                 <  1:1432284/91‑1 (MQ=255)
            aaTATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAttt                                                                                >  1:1159743/1‑91 (MQ=255)
             atatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAAGCCTAGGTATTCATGAGTACCTAGGAtttt                                                                               <  2:1262102/91‑1 (MQ=255)
              tatTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGAttttt                                                                              <  1:724154/91‑1 (MQ=255)
                ttCCACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGCTTTTTAc                                                                            >  2:3500275/1‑91 (MQ=255)
                ttCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTAc                                                                            >  1:3049592/1‑91 (MQ=255)
                  cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                          >  1:665981/1‑91 (MQ=255)
                  cAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGc                                                                          >  2:1949170/1‑91 (MQ=255)
                    aCGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCAt                                                                        <  2:1044375/91‑1 (MQ=255)
                      ggAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGa                                                                      <  1:925528/91‑1 (MQ=255)
                       gAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTATGCATGAt                                                                     >  2:82937/1‑91 (MQ=255)
                        aaTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATa                                                                    <  1:408850/91‑1 (MQ=255)
                        aaTCAATATTAGAAAAATGGTTAAATCTTTGCAATAGCAAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATa                                                                    <  2:1874918/91‑1 (MQ=255)
                         aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                                   >  1:3583339/1‑91 (MQ=255)
                         aTCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaa                                                                   >  1:3910800/1‑91 (MQ=255)
                          tCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATaaa                                                                  >  1:513235/1‑91 (MQ=255)
                           cAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAt                                                                 >  1:1295020/1‑91 (MQ=255)
                           cAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAt                                                                 >  1:2026241/1‑91 (MQ=255)
                            aaTATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATa                                                                <  2:901310/91‑1 (MQ=255)
                             atatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAc                                                               <  2:3380257/91‑1 (MQ=255)
                              tatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATAcc                                                              >  2:402027/1‑90 (MQ=255)
                              tatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACa                                                              >  1:3262182/1‑91 (MQ=255)
                              tatTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACa                                                              >  1:2951392/1‑91 (MQ=255)
                               atTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaa                                                             >  1:3355268/1‑91 (MQ=255)
                                atAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAAACCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAaaaaaaaa                                                            <  2:3211285/90‑1 (MQ=39)
                                 tAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaa                                                           <  1:801346/91‑1 (MQ=255)
                                 tAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACaaaa                                                           <  1:2659056/91‑1 (MQ=255)
                                     aaaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAag                                                       >  1:2775306/1‑91 (MQ=255)
                                     aaaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAag                                                       >  1:2754011/1‑91 (MQ=255)
                                      aaaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAaga                                                      >  1:3108110/1‑91 (MQ=255)
                                       aaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGTATGATAAAAAAAAAAAAagag                                                     <  2:2386504/91‑1 (MQ=255)
                                       aaTGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAagag                                                     <  1:1006923/91‑1 (MQ=255)
                                          ggTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAgg                                                  <  2:2974533/91‑1 (MQ=255)
                                          ggTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAgg                                                  <  1:1363476/91‑1 (MQ=255)
                                           gTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGa                                                 >  1:1636329/1‑91 (MQ=255)
                                            ttAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGAt                                                <  1:993343/91‑1 (MQ=255)
                                            ttAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAAAGAGGAt                                                >  1:1561834/1‑91 (MQ=255)
                                            ttAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTAACTAGGATTTTTACGCATGAAAAAAAAAAAAAAAGAGAGGAt                                                <  2:1638560/91‑1 (MQ=255)
                                              aaaTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAg                                              >  2:728326/1‑91 (MQ=255)
                                               aaTCTTTACAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGactat                                             <  1:1978864/91‑5 (MQ=255)
                                                aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  1:3647891/91‑1 (MQ=255)
                                                aTCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTc                                            <  1:1817389/91‑1 (MQ=255)
                                                 tCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTcc                                           <  1:1873843/91‑1 (MQ=255)
                                                   ttttCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCctc                                         <  2:2551039/91‑1 (MQ=255)
                                                      tCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCa                                      <  1:3858668/91‑1 (MQ=255)
                                                       cAATAGCTAAATGAAAAAATCCGAGGTATTCATGAGTGCCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGACTATCCtctctt                                     <  2:1588381/91‑3 (MQ=255)
                                                            gCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAAt                                <  2:3242214/91‑1 (MQ=255)
                                                            gCTAAATGAAAAAATCCTAGGTATTCATGAGAACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAAt                                >  1:3273061/1‑91 (MQ=255)
                                                              aaaaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAg                              <  2:1240120/90‑1 (MQ=255)
                                                               aaaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGt                             <  1:2152680/91‑1 (MQ=255)
                                                                aaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGtt                            >  1:2046599/1‑91 (MQ=255)
                                                                aaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGtt                            >  1:3753239/1‑91 (MQ=255)
                                                                aaTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAAAACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGtt                            <  2:2003525/91‑1 (MQ=255)
                                                                 aTGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTta                           <  1:3496792/91‑1 (MQ=255)
                                                                   gAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTtaat                         <  2:593655/91‑1 (MQ=255)
                                                                        aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAt                    >  1:3839273/1‑91 (MQ=255)
                                                                        aaTCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAt                    <  2:1108563/91‑1 (MQ=255)
                                                                         aTCCTAGGTATTCATGAGTACCGAGGATTTTTACGCATGAGAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAAtt                   <  2:390041/91‑1 (MQ=255)
                                                                            cTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTa                >  1:336581/1‑91 (MQ=255)
                                                                             tAGGTATTCATGAGTACCTAGGATTTTTACGCATGATATATACAAAAAAAGAGAGGATAGTCCTCTCAAATAATAGTTAATAAAATTGTAg               <  1:3442982/91‑1 (MQ=255)
                                                                             tAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCACTCAAACAATGGTTAATAAAATTGTAg               >  2:547991/1‑91 (MQ=255)
                                                                              aGGTATTCATGAGTACCTAAGATGTTTACGCATGATAAATACAAAAAGAGGGAGGAGAGTCCTCGCAAACAATATTTAATAAAATTGTAGc              <  2:1016684/91‑1 (MQ=255)
                                                                              aGGTAGTAATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGc              <  2:315751/91‑1 (MQ=255)
                                                                                gTATTCATGAGTACCTGGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCa            <  2:3787646/91‑1 (MQ=255)
                                                                                 tactCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAt           >  2:1075344/4‑91 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         <  2:1387608/91‑1 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  1:2994378/1‑91 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  2:475995/1‑91 (MQ=255)
                                                                                   ttCATGAGTACCTAGGATTTTTACGAATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAGTAAAATTGTAGCCAttt         <  2:170189/91‑1 (MQ=255)
                                                                                   ctCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAGATTGTAGCCAttt         >  1:3956405/2‑91 (MQ=255)
                                                                                   ctCATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCAttt         >  1:150960/2‑91 (MQ=255)
                                                                                     cATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttct       >  1:2679926/1‑91 (MQ=255)
                                                                                     cATGAATACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttct       >  1:2015200/1‑91 (MQ=255)
                                                                                     aaagAGTACCTAGGATTTTTACGAATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATttct       <  2:1487814/88‑1 (MQ=255)
                                                                                         aGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCa   >  2:1284973/1‑91 (MQ=255)
                                                                                          gTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAg  <  1:2276458/91‑1 (MQ=255)
                                                                                          gTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAg  <  1:2856044/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACATTATAGAAAAATATTCAACGGAATCAATATTAGAAAAATGGTTAAATCTTTTCAATAGCTAAATGAAAAAATCCTAGGTATTCATGAGTACCTAGGATTTTTACGCATGATAAATACAAAAAAAGAGAGGATAGTCCTCTCAAACAATAGTTAATAAAATTGTAGCCATTTCTTCCAG  >  USA300TCH1516_ALE/1042706‑1042886

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