Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,774,463 A→T 5.0% intergenic (‑93/+23) USA300HOU_RS14290 ← / ← clfB S‑formylglutathione hydrolase/Clumping factor B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,774,4630AT 5.0% 799.8 / 22.6 238intergenic (‑93/+23)USA300HOU_RS14290/clfBS‑formylglutathione hydrolase/Clumping factor B
Reads supporting (aligned to +/- strand):  ref base A (109/117);  new base T (3/9);  total (112/126)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

Only 100 of 247 total aligned reads displayed.

GTCATGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCCATC  >  USA300TCH1516_ALE/2774373‑2774552
                                                                                          |                                                                                         
gTCATGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACaaa                                                                                           >  2:1263095/1‑91 (MQ=255)
gTCATGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACaaa                                                                                           <  1:1419073/91‑1 (MQ=255)
 tCATGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAAt                                                                                          <  1:2227290/91‑1 (MQ=255)
   aTGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGt                                                                                        <  2:1446006/91‑1 (MQ=255)
   aTGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGt                                                                                        >  2:1923832/1‑91 (MQ=255)
     gACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTAt                                                                                      >  2:1775412/1‑91 (MQ=255)
      aCCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACATTTGTAtt                                                                                     <  1:1316370/91‑1 (MQ=255)
        cccTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACATTTGTATTCg                                                                                   <  1:3778085/91‑1 (MQ=255)
        cccTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACATTTGTATTCg                                                                                   <  1:3248277/91‑1 (MQ=255)
         ccTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCgg                                                                                  <  2:2506752/91‑1 (MQ=255)
         ccTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCgg                                                                                  >  1:3675583/1‑91 (MQ=255)
          cTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACATTTGTATTCggt                                                                                 <  2:2369323/91‑2 (MQ=255)
          cTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGc                                                                                 >  1:1032684/1‑91 (MQ=255)
          cTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGc                                                                                 >  1:3100541/1‑91 (MQ=255)
           tttCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGcc                                                                                >  2:3018387/1‑91 (MQ=255)
              cTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTaa                                                                             <  1:1882088/91‑1 (MQ=255)
                  ttCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTaaaaaa                                                                         >  2:2597907/1‑91 (MQ=255)
                   tCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAg                                                                        >  2:1524624/1‑91 (MQ=255)
                    cTACCTTATCAACATTTTACGCTTAATAAATTCACTTTACAATCATTTTCAACAAAAAAACCGAATACACACGTCTTCGGCCTAAAACAGt                                                                       >  2:1674856/1‑91 (MQ=255)
                      aCTTTATCAACATTTTACGCTTAATCAATTCACTTTACAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTAt                                                                     >  2:1374353/1‑91 (MQ=255)
                        tttATCAACATTTTACGCTTAATCAATTCGCTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTa                                                                      <  1:3308935/88‑1 (MQ=255)
                        tttATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTa                                                                      >  2:3308935/1‑88 (MQ=255)
                        tttATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTAttt                                                                   >  2:343229/1‑91 (MQ=255)
                        tttATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTAttt                                                                   >  2:3190899/1‑91 (MQ=255)
                         ttATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTa                                                                  >  2:1565856/1‑91 (MQ=255)
                           aTCAACATTTTACGCTTAATCAATTTACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAACGTATTTACg                                                                >  2:1387302/1‑91 (MQ=255)
                           aTCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAGACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACg                                                                <  2:1308995/91‑1 (MQ=255)
                           aTCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACg                                                                <  1:2898234/91‑1 (MQ=255)
                            tCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGc                                                               <  2:1713314/91‑1 (MQ=255)
                            tCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGc                                                               <  1:3385450/91‑1 (MQ=255)
                                cATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCtttt                                                           >  1:2106728/1‑91 (MQ=255)
                                cATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCtttt                                                           >  1:934851/1‑91 (MQ=255)
                                  ttttACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTc                                                         <  1:1418875/91‑1 (MQ=255)
                                   tttACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCt                                                        <  1:2170650/91‑1 (MQ=255)
                                    ttACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCtt                                                       <  2:948503/91‑1 (MQ=255)
                                    ttACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCtt                                                       >  2:884546/1‑91 (MQ=255)
                                     tACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCttt                                                      <  1:501427/91‑1 (MQ=255)
                                      aCGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTa                                                     <  2:960505/91‑1 (MQ=255)
                                      aCGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTa                                                     >  1:3482526/1‑91 (MQ=255)
                                      aCGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTa                                                     <  1:2905383/91‑1 (MQ=255)
                                      aCGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTa                                                     >  1:347210/1‑91 (MQ=255)
                                       cGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTCAAAAAGTATTTACGCTTTTTCTTTAt                                                    >  2:3764169/1‑91 (MQ=255)
                                       cGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTAt                                                    >  2:684385/1‑91 (MQ=255)
                                       cGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTAt                                                    >  2:3676672/1‑91 (MQ=255)
                                       cGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTAt                                                    >  2:1408148/1‑91 (MQ=255)
                                        gCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATg                                                   <  1:1381611/91‑1 (MQ=255)
                                          ttAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGAt                                                 <  2:3112734/91‑1 (MQ=255)
                                          ttAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGAt                                                 <  1:2773920/91‑1 (MQ=255)
                                           tAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATc                                                <  2:625887/91‑1 (MQ=255)
                                            aaTCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATct                                               <  2:2881638/91‑1 (MQ=255)
                                              tCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttg                                             >  2:3336164/1‑91 (MQ=255)
                                               cAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgc                                            >  1:599921/1‑91 (MQ=255)
                                               cAATTCACTTTAAAATCATTTACAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgc                                            <  2:3140033/91‑1 (MQ=255)
                                                aaTTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgct                                           <  1:1409290/91‑1 (MQ=255)
                                                aaTTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgct                                           <  1:1401804/91‑1 (MQ=255)
                                                aaTTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgct                                           <  2:1201942/91‑1 (MQ=255)
                                                 aTTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgctt                                          <  2:1138719/91‑1 (MQ=255)
                                                 aTTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATcttgctt                                          <  2:2152537/91‑1 (MQ=255)
                                                    cACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCg                                       <  1:3164019/91‑1 (MQ=255)
                                                    cACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCg                                       <  1:1600507/91‑1 (MQ=255)
                                                    aaCTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCg                                       <  2:285936/90‑1 (MQ=255)
                                                     aCTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCttgcttgcat                                      >  1:263056/1‑89 (MQ=255)
                                                      cTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGtt                                     >  2:1602948/1‑91 (MQ=255)
                                                       tttAAAATCATTTTCAACAAAAAATCCGAATACAAATGTATTCAGCCTAAAAAAGTATTTACGCTTTTTCTTTATGTTCTTGCTTGCGttt                                    <  2:544591/91‑1 (MQ=255)
                                                       tttAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAGAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGttt                                    >  2:3409796/1‑91 (MQ=255)
                                                             aTCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTaaa                              >  1:827660/1‑91 (MQ=255)
                                                              tCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAAc                             <  1:2509042/91‑1 (MQ=255)
                                                              tCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTCGATGATCTTGCTTGCGTGTTCTAAAc                             >  2:3369636/1‑91 (MQ=255)
                                                                aTTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACaa                           >  2:2164642/1‑91 (MQ=255)
                                                                aTTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACaa                           <  2:2074795/91‑1 (MQ=255)
                                                                aTTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACaa                           >  1:763375/1‑91 (MQ=255)
                                                                 ttttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAAt                          <  1:884546/91‑1 (MQ=255)
                                                                  tttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATa                         <  2:3738785/91‑1 (MQ=255)
                                                                  tttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATa                         >  1:243384/1‑91 (MQ=255)
                                                                  tttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATa                         >  1:68499/1‑91 (MQ=255)
                                                                   ttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAg                        <  1:848025/91‑1 (MQ=255)
                                                                   ttCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAg                        <  2:1746495/91‑1 (MQ=255)
                                                                    tCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGt                       >  1:2008071/1‑91 (MQ=255)
                                                                    tCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGt                       <  2:2660842/91‑1 (MQ=255)
                                                                         aaaaaaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGa                  >  1:2382385/1‑91 (MQ=255)
                                                                           aaaaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGa                  >  2:479935/1‑89 (MQ=255)
                                                                           aaaaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATc                <  2:582397/91‑1 (MQ=255)
                                                                           aaaaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATc                <  2:3580322/91‑1 (MQ=255)
                                                                           aaaaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATc                <  2:2226574/91‑1 (MQ=255)
                                                                             aaaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCt              <  2:293953/91‑1 (MQ=255)
                                                                              aaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCta             >  1:1518524/1‑91 (MQ=255)
                                                                              aaCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCta             <  2:2000857/91‑1 (MQ=255)
                                                                               aCCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCtaa            >  1:2017931/1‑91 (MQ=255)
                                                                                ccGAATACATTTGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCtaat           >  1:1093104/1‑91 (MQ=255)
                                                                                ccGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCtaat           >  1:3853374/1‑91 (MQ=255)
                                                                                 cGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCtaata          >  2:1597922/1‑91 (MQ=255)
                                                                                 cGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCtaata          >  2:1081070/1‑91 (MQ=255)
                                                                                    aTACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATg       >  1:1070051/1‑91 (MQ=255)
                                                                                    aTACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATg       <  1:3336164/91‑1 (MQ=255)
                                                                                     tACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGc      >  1:2808073/1‑91 (MQ=255)
                                                                                     tACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGc      <  1:2294423/91‑1 (MQ=255)
                                                                                       cAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCca    >  1:2388204/1‑91 (MQ=255)
                                                                                       cAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCca    <  2:3175812/91‑1 (MQ=255)
                                                                                         aaTGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCcatc  >  2:184449/1‑91 (MQ=255)
                                                                                         aaTGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCcatc  >  1:1822464/1‑91 (MQ=255)
                                                                                          |                                                                                         
GTCATGACCCCTTTCTAATTCTACTTTATCAACATTTTACGCTTAATCAATTCACTTTAAAATCATTTTCAACAAAAAAACCGAATACAAATGTATTCGGCCTAAAAAAGTATTTACGCTTTTTCTTTATGATCTTGCTTGCGTTTTCTAAACAATAGTAATGATCCTAATAATGCCATC  >  USA300TCH1516_ALE/2774373‑2774552

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