Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,929,570 T→G 6.1% E32A (GAG→GCG)  USA300TCH1516_01793 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,929,5700TG 6.1% 993.9 / 44.6 303E32A (GAG→GCG) USA300TCH1516_01793hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (136/144);  new base G (16/3);  total (156/147)
Fisher's exact test for biased strand distribution p-value = 3.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

ATTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTG  >  USA300TCH1516_ALE/1929480‑1929660
                                                                                          |                                                                                          
aTTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACt                                                                                            <  1:1170506/91‑1 (MQ=255)
 ttGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTc                                                                                           >  1:3805135/1‑91 (MQ=255)
   gAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgc                                                                                         >  1:215147/1‑91 (MQ=255)
   gAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTcgc                                                                                         >  2:1784278/1‑91 (MQ=255)
    aaTGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATAGCcgcg                                                                                        >  2:3646804/1‑91 (MQ=255)
      tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAGACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                      >  2:3088486/1‑91 (MQ=255)
      tGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGta                                                                                      >  2:3781229/1‑91 (MQ=255)
       gTATCCGTGCTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtat                                                                                     <  1:3553495/91‑1 (MQ=255)
       gTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATTATACGCGAGtat                                                                                     <  2:141441/91‑1 (MQ=255)
       gTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGtat                                                                                     <  1:3843941/91‑1 (MQ=255)
         aCCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATaa                                                                                   <  2:753362/91‑1 (MQ=255)
          cccGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAt                                                                                  <  1:3713162/91‑1 (MQ=255)
            cGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAAttt                                                                                <  1:2516005/91‑1 (MQ=255)
            cGTGTTTCGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGGGTATAAttt                                                                                <  1:1971275/91‑1 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               >  1:1930698/1‑91 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               >  1:1836514/1‑91 (MQ=255)
             gtgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTc                                                                               <  1:2801540/91‑1 (MQ=255)
              tgtTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCa                                                                              >  2:1003515/1‑91 (MQ=255)
                ttttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCAtt                                                                            <  2:589/91‑1 (MQ=255)
                  ttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATCATTTCATTcc                                                                          >  1:2702915/1‑91 (MQ=255)
                  ttGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTcc                                                                          >  2:2585066/1‑91 (MQ=255)
                    gATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTa                                                                        >  2:2004289/1‑91 (MQ=255)
                      tGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAgg                                                                      >  1:535743/1‑91 (MQ=255)
                       gTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGc                                                                     <  2:568641/91‑1 (MQ=255)
                        ttGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTATGAAATCATACTCGCGTATAATTTCATTCCTAGGCt                                                                    >  2:3567225/1‑91 (MQ=255)
                         tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                   >  1:492013/1‑91 (MQ=255)
                         tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                   <  1:25951/91‑1 (MQ=255)
                         tGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCtt                                                                   <  1:168135/91‑1 (MQ=255)
                         tGCTGTATCCATTTTCGTAGATTTTGTAACCGCCGTTTCCGGACGTCGAAGGTGGGAAATCATACTCGCGTATGATTTCATTCCGAGGCtt                                                                   <  2:513495/91‑1 (MQ=255)
                           cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                 <  2:2864498/91‑1 (MQ=255)
                           cTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAc                                                                 >  2:2843356/1‑91 (MQ=255)
                            tttATCCATTTTCCTACATTTTATAACCGGCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAcc                                                                <  1:3166162/91‑1 (MQ=255)
                            tttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTAcc                                                                <  1:419532/91‑1 (MQ=255)
                             ttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCa                                                               <  1:1311260/91‑1 (MQ=255)
                             ttATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCa                                                               >  1:1908874/1‑91 (MQ=255)
                              tATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAt                                                              >  1:2288160/1‑91 (MQ=255)
                              tATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAt                                                              >  1:340571/1‑91 (MQ=255)
                               aTCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtt                                                             <  2:2164778/91‑1 (MQ=255)
                                tCCATTTTCCTACATTTTATAGCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttt                                                            <  2:2527183/91‑1 (MQ=255)
                                 ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                           >  1:742943/1‑91 (MQ=255)
                                 ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                           >  2:1064638/1‑91 (MQ=255)
                                 ccATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAtttt                                                           >  2:902521/1‑91 (MQ=255)
                                  cATTTTCCTACATTTTATAACCGCCATTTACAGACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCAttttt                                                          <  2:12926/91‑1 (MQ=255)
                                      ttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatat                                                      <  2:2995734/91‑1 (MQ=255)
                                      ttCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatat                                                      <  1:133392/91‑1 (MQ=255)
                                       tCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTtatata                                                     >  2:3342015/1‑91 (MQ=255)
                                        ccTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATaa                                                    >  1:137925/1‑91 (MQ=255)
                                          tACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATc                                                  <  2:2387918/91‑1 (MQ=255)
                                          tACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATc                                                  <  1:448141/91‑1 (MQ=255)
                                           aCATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCt                                                 >  2:3378648/1‑91 (MQ=255)
                                            cATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTa                                                <  1:3045051/91‑1 (MQ=255)
                                            cATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTa                                                >  1:203759/1‑91 (MQ=255)
                                             aTTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTaa                                               <  1:2502457/91‑1 (MQ=255)
                                              ttttATAACCGCCATTTACAAACCTCGAATGTGTGAAATCATACTCGCGCATAATTTCATTCCTAGGCGTACCATCTTTATATTATCTAAc                                              >  2:156132/1‑91 (MQ=255)
                                                 tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                           <  1:3483734/91‑1 (MQ=255)
                                                 tataACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAt                                           <  1:1198498/91‑1 (MQ=255)
                                                   taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         >  1:2100055/1‑91 (MQ=255)
                                                   taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         >  2:2276434/1‑91 (MQ=255)
                                                   taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         <  2:1111382/91‑1 (MQ=255)
                                                   taACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCAttt                                         >  1:2892575/1‑91 (MQ=255)
                                                    aaCCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTg                                        <  1:1595264/91‑1 (MQ=255)
                                                      ccGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTa                                      <  2:138913/91‑1 (MQ=255)
                                                        gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                    >  1:217684/1‑91 (MQ=255)
                                                        gCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAAc                                    >  1:3327480/1‑91 (MQ=255)
                                                         ccATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACt                                   >  2:2544074/1‑91 (MQ=255)
                                                           aTTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTa                                 >  1:2411425/1‑91 (MQ=255)
                                                             ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                               >  2:1603925/1‑91 (MQ=255)
                                                             ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                               >  2:2202960/1‑91 (MQ=255)
                                                             ttACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTaaa                               <  1:3281854/91‑1 (MQ=255)
                                                               aCAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAAct                             >  2:726830/1‑91 (MQ=255)
                                                                 aaaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAActct                           <  1:3760238/91‑1 (MQ=255)
                                                                  aaCGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTg                          >  1:103081/1‑91 (MQ=255)
                                                                    cGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATGTAATCTAACCATTTGTAACTTAAACTCTGaa                        >  1:3370760/1‑91 (MQ=255)
                                                                        gAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaa                    >  1:2988010/1‑91 (MQ=255)
                                                                         aaGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTaaa                   >  2:1860100/1‑91 (MQ=255)
                                                                          aGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAAt                  >  2:3181956/1‑91 (MQ=255)
                                                                          aGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAAt                  >  2:3283274/1‑91 (MQ=255)
                                                                           ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                 <  1:3609804/91‑1 (MQ=255)
                                                                           ggTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATg                 <  1:1505871/91‑1 (MQ=255)
                                                                            gtgtGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGa                <  1:635341/91‑1 (MQ=255)
                                                                            gtgtGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGa                <  1:1639295/91‑1 (MQ=255)
                                                                              gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc              >  1:2942677/1‑91 (MQ=255)
                                                                              gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc              <  1:697400/91‑1 (MQ=255)
                                                                              gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc              <  2:412038/91‑1 (MQ=255)
                                                                              gtgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtc              <  2:446503/91‑1 (MQ=255)
                                                                               tgAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtct             >  1:807454/1‑91 (MQ=255)
                                                                                gAAATTATACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt            >  1:2854533/1‑91 (MQ=17)
                                                                                gAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtctt            <  2:1539484/91‑1 (MQ=255)
                                                                                  aaTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGAtcttct          <  1:1713736/91‑1 (MQ=255)
                                                                                   aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt         <  2:2627809/91‑1 (MQ=255)
                                                                                   aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt         >  2:2468370/1‑91 (MQ=255)
                                                                                   aTCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATcttctt         <  1:2355908/91‑1 (MQ=255)
                                                                                    tCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttcttt        <  2:415866/91‑1 (MQ=255)
                                                                                     taTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       >  2:275431/2‑91 (MQ=14)
                                                                                     taTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       >  1:2431519/2‑91 (MQ=14)
                                                                                     cATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttc       <  1:2563323/91‑1 (MQ=255)
                                                                                      aTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTTAACTAAATGATCttctttct      >  2:2942752/1‑91 (MQ=14)
                                                                                      aTACGCGAGTATGATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCttctttct      >  2:1690920/1‑91 (MQ=14)
                                                                                       tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc     >  1:2120692/1‑91 (MQ=255)
                                                                                       tACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTtctc     >  1:134852/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTGAATGTACCCGTGTTTTGATGTTGCTTTATCCATTTTCCTACATTTTATAACCGCCATTTACAAACGTCGAAGGTGTGAAATCATACTCGCGTATAATTTCATTCCTAGGCTTACCATTTTTATATAATCTAACCATTTGTAACTTAAACTCTGAACTAAATGATCTTCTTTCTCTTG  >  USA300TCH1516_ALE/1929480‑1929660

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