Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,743 T→G 5.1% intergenic (+234/+29) USA300HOU_RS06065 → / ← USA300HOU_RS06070 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,218,7430TG 5.1% 717.4 / 21.7 215intergenic (+234/+29)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (107/96);  new base G (7/4);  total (115/100)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 226 total aligned reads displayed.

GATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGTTTTTTGAT  >  USA300TCH1516_ALE/1218653‑1218833
                                                                                           |                                                                                          
gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                                            <  1:3486335/91‑1 (MQ=255)
    cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                                        <  1:2374706/91‑1 (MQ=255)
    cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTAtggtga                                                                                       <  2:3889438/91‑5 (MQ=18)
    cGTTATTCAACAACACGATGCTTCTTTCCACAGGACGCATTAAAGAATAAATAAACTCCTCACCATAAGAAAGAAATACTTTCTTTTGTGa                                                                                        <  2:3078643/91‑1 (MQ=255)
      ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                                      >  2:745898/1‑91 (MQ=255)
       tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                                     <  2:2922554/91‑1 (MQ=17)
       tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                                     <  1:3837020/91‑1 (MQ=17)
        aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                                    <  2:2006080/91‑1 (MQ=255)
        aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                                    >  1:3919822/1‑91 (MQ=255)
         tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                                   <  1:3225306/91‑1 (MQ=255)
           tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                                 <  1:3349698/91‑1 (MQ=255)
           tAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGttt                                                                                 <  1:2772467/91‑1 (MQ=255)
             aaaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTAATTTCTTTTGTGAGGTGTTTAt                                                                               >  2:3237558/1‑91 (MQ=255)
                acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGGGGTGTTTATTTg                                                                            <  1:3054373/91‑1 (MQ=255)
                acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                                            <  2:3654933/91‑1 (MQ=255)
                 cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                                           <  2:1184768/91‑1 (MQ=255)
                  acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                                          >  1:1895818/1‑91 (MQ=255)
                   cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                                         >  1:642275/1‑91 (MQ=255)
                   cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                                         <  1:433252/91‑1 (MQ=255)
                    gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                                        >  2:2592315/1‑91 (MQ=255)
                     aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                                       >  1:3612745/1‑91 (MQ=255)
                      tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGCAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                                      >  2:266419/1‑91 (MQ=255)
                      tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                                      >  2:2798041/1‑91 (MQ=255)
                      tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                                      <  2:3172885/91‑1 (MQ=255)
                       aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                                     >  1:1698586/1‑91 (MQ=255)
                          tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                                  >  2:1650738/1‑91 (MQ=255)
                          tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                                  <  2:2079422/91‑1 (MQ=255)
                           cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                                 <  2:3580994/91‑1 (MQ=255)
                              ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                              >  1:1010180/1‑91 (MQ=255)
                                cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                                            >  2:1783444/1‑91 (MQ=255)
                                 aaTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATg                                                           >  2:2278476/1‑91 (MQ=255)
                                 aaTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATg                                                           >  1:3939846/1‑91 (MQ=255)
                                  aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAATGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                          >  1:3423208/1‑91 (MQ=255)
                                  aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                          >  1:3010014/1‑91 (MQ=255)
                                  aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                                          >  2:1836425/1‑91 (MQ=255)
                                    gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                                        <  2:3851415/91‑1 (MQ=255)
                                    gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                                        <  2:1682575/91‑1 (MQ=255)
                                    gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                                        <  2:3839369/91‑1 (MQ=255)
                                     aaCGCATTAAATAATAAATAAGCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAAGTGTTTATTTGACGACAACTTGATTATGATTg                                                       <  2:483669/91‑1 (MQ=255)
                                      aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                                      >  1:2504358/1‑91 (MQ=255)
                                      aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                                      <  2:2728711/91‑1 (MQ=255)
                                       cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                     <  1:3099512/91‑1 (MQ=255)
                                       cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                     <  2:1027553/91‑1 (MQ=255)
                                       cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                     <  2:60259/91‑1 (MQ=255)
                                       cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                                     <  1:2280618/91‑1 (MQ=255)
                                         cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                                   <  2:2033507/91‑1 (MQ=255)
                                             aaataataaataaACACCTCACCATAAGAAAGGCCTACTTTCTTTTCTGAGGTGTTTATCTGTCGACAACTTGATTATGATTGTTTCAttt                                               >  1:3157127/1‑91 (MQ=255)
                                             aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                               <  1:1197190/91‑1 (MQ=255)
                                             aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                                               <  2:1365437/91‑1 (MQ=255)
                                               ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                             >  1:357499/1‑91 (MQ=255)
                                               ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                             >  1:1933328/1‑91 (MQ=255)
                                                taataaataaACTCCTCACCATGAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                                            <  2:2131532/91‑1 (MQ=255)
                                                  ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                                          <  1:1612640/91‑1 (MQ=255)
                                                   taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                                         <  2:2626018/91‑1 (MQ=255)
                                                   taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                                         <  2:3382182/91‑1 (MQ=255)
                                                      ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                                      <  1:2930624/91‑1 (MQ=255)
                                                       taaaCACCTTACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGATTAAg                                     >  1:1121913/1‑91 (MQ=255)
                                                         aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                                   <  1:701587/91‑1 (MQ=255)
                                                           cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                                 >  2:1581888/1‑91 (MQ=17)
                                                            acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                                <  1:2638266/91‑1 (MQ=255)
                                                             ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                               <  2:1144837/91‑1 (MQ=255)
                                                              cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                              >  1:1275346/1‑91 (MQ=255)
                                                               tCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAAt                             <  1:1061682/91‑1 (MQ=255)
                                                                cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc                            >  1:255728/1‑91 (MQ=255)
                                                                 aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                           <  2:1531724/91‑1 (MQ=255)
                                                                 aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                           >  1:3302163/1‑91 (MQ=255)
                                                                 aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                           >  2:2489810/1‑91 (MQ=255)
                                                                  ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                          <  1:470464/91‑1 (MQ=255)
                                                                  ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCCTAATCAg                          >  2:2408818/1‑91 (MQ=255)
                                                                  ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                          >  1:3709309/1‑91 (MQ=255)
                                                                   caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg                          >  2:422806/4‑91 (MQ=21)
                                                                      aaGAAAGTGCTACTTTCTTTTGTGAGGTTTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                      <  2:938597/91‑1 (MQ=255)
                                                                      aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                      <  2:2906614/91‑1 (MQ=255)
                                                                      aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt                      >  1:987422/1‑91 (MQ=255)
                                                                       aGAAAGGACTACTTTCTCTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATCGTTTCATTTCTTGAATAAGTTCATAATCAGGTGTa                     >  2:1918019/1‑91 (MQ=255)
                                                                        gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg                    >  2:107104/1‑91 (MQ=255)
                                                                        gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg                    <  2:3733675/91‑1 (MQ=255)
                                                                          aaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCa                  >  2:163644/1‑91 (MQ=255)
                                                                            ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa                >  2:51291/1‑91 (MQ=255)
                                                                            ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa                >  1:3088054/1‑91 (MQ=255)
                                                                              actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGCTTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa              >  2:888746/1‑91 (MQ=255)
                                                                              actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa              >  1:2740194/1‑91 (MQ=255)
                                                                              actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa              <  1:966470/91‑1 (MQ=255)
                                                                               gtcctttctttTGTGAGGTGTTTATTTGTCGACCACTAGATTATGACTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa             >  1:3150141/4‑91 (MQ=255)
                                                                               gtcctttcttaTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa             >  2:2384605/4‑91 (MQ=255)
                                                                                tactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaaa            <  2:3631500/91‑1 (MQ=255)
                                                                                 actTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTTATTATGTTCATAATTAGTTGTAGCATCCAAAg           >  2:3702509/1‑91 (MQ=255)
                                                                                  ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt          <  1:3724483/91‑1 (MQ=255)
                                                                                  ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt          <  1:521358/91‑1 (MQ=255)
                                                                                  ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt          <  1:747902/91‑1 (MQ=255)
                                                                                  ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt          <  1:1400499/91‑1 (MQ=255)
                                                                                  ctttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGt          >  1:2035825/1‑91 (MQ=255)
                                                                                   tttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGtt         >  2:859695/1‑91 (MQ=255)
                                                                                   tttctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGtt         >  1:2589427/1‑91 (MQ=255)
                                                                                     tctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGtttt       <  1:12224/91‑1 (MQ=255)
                                                                                      ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt      <  1:2184910/91‑1 (MQ=255)
                                                                                      ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt      <  2:3036007/91‑1 (MQ=255)
                                                                                      ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt      <  2:3035700/91‑1 (MQ=255)
                                                                                      ctttTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGttttt      >  2:2680315/1‑91 (MQ=255)
                                                                                          ttGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGTTTTTTGa   >  1:1772107/1‑91 (MQ=255)
                                                                                           |                                                                                          
GATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAAAGTTTTTTGAT  >  USA300TCH1516_ALE/1218653‑1218833

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: