Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,730 G→T 11.9% intergenic (+221/+42) 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,7300GT11.9% 673.4 / 56.4 208intergenic (+221/+42)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (88/94);  new base T (8/17);  total (96/112)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

AACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAC  >  USA300TCH1516_ALE/1218641‑1218820
                                                                                         |                                                                                           
aaCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGa                                                                                            <  2:3798108/91‑1 (MQ=255)
  cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                                          >  2:2974548/1‑91 (MQ=255)
       atgatgATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActtt                                                                                     >  1:3001245/1‑91 (MQ=255)
         gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                   <  1:1200044/91‑1 (MQ=255)
         gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                   <  1:3944412/91‑1 (MQ=255)
         gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                                  <  2:2913815/91‑1 (MQ=25)
          atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                 <  1:798432/91‑2 (MQ=21)
          atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                                 <  2:1249390/91‑2 (MQ=21)
           tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAGCGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttat                                                                                <  2:2038579/91‑3 (MQ=21)
            gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                               <  1:3486335/91‑1 (MQ=255)
             atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatgg                                                                              <  2:1520314/91‑5 (MQ=21)
              tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                             <  2:1649225/91‑1 (MQ=255)
                cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                           <  1:2374706/91‑1 (MQ=255)
                  ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                         >  2:745898/1‑91 (MQ=255)
                  ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                         <  1:2609939/91‑1 (MQ=255)
                   tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                        <  1:330213/91‑1 (MQ=17)
                   tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGgt                                                                        <  2:2922554/91‑1 (MQ=17)
                    aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                       <  2:2006080/91‑1 (MQ=255)
                      ttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGtt                                                                     <  1:2209285/91‑1 (MQ=255)
                          aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                                 >  2:3409846/1‑91 (MQ=255)
                           aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                >  1:3925925/1‑91 (MQ=255)
                           aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                                >  1:389622/1‑91 (MQ=255)
                            acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGGGGTGTTTATTTg                                                               <  1:3054373/91‑1 (MQ=255)
                            acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                               >  1:2651647/1‑91 (MQ=255)
                            acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                               <  2:3654933/91‑1 (MQ=255)
                             cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                              <  2:1184768/91‑1 (MQ=255)
                             cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                              <  2:1558508/91‑1 (MQ=255)
                             cacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGt                                                              >  1:11082/1‑91 (MQ=255)
                              acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                             >  1:1895818/1‑91 (MQ=255)
                               cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                            >  1:642275/1‑91 (MQ=255)
                                gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                           >  2:2592315/1‑91 (MQ=255)
                                 aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                          <  2:2909669/91‑1 (MQ=255)
                                  tACTTCGTTTCAATGAACGCATTAAATAATAAATACACACCTCACCATACGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                         >  2:3165268/1‑91 (MQ=255)
                                  tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGCAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                         >  2:266419/1‑91 (MQ=255)
                                  tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACa                                                         >  2:2798041/1‑91 (MQ=255)
                                  tACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTCGTGAGGTGTTTATTTGTCGACa                                                         >  2:395805/1‑91 (MQ=255)
                                   aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                        >  1:1698586/1‑91 (MQ=255)
                                    cTTCGTTTCAATGAACGCATTAAATGATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAAc                                                       <  1:3681273/91‑1 (MQ=255)
                                     ttCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACt                                                      <  1:716904/91‑1 (MQ=255)
                                      tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                     <  2:2079422/91‑1 (MQ=255)
                                      tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGAGAGGTGTTTATTTGTCGACAACtt                                                     >  2:1650773/1‑91 (MQ=255)
                                       cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                    <  2:3580994/91‑1 (MQ=255)
                                        gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                                   >  1:577277/1‑91 (MQ=255)
                                        gTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGa                                                   >  1:2048889/1‑91 (MQ=255)
                                          ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                 <  2:2492768/91‑1 (MQ=255)
                                          ttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAtt                                                 >  1:1010180/1‑91 (MQ=255)
                                           tCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTa                                                <  1:1831263/91‑1 (MQ=255)
                                            cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                               >  2:1783444/1‑91 (MQ=255)
                                              aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAATGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                             >  1:3423208/1‑91 (MQ=255)
                                              aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                             >  2:1836425/1‑91 (MQ=255)
                                               tGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAt                                            >  2:821556/1‑91 (MQ=255)
                                                gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                           <  2:3851415/91‑1 (MQ=255)
                                                gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                           <  2:3106492/91‑1 (MQ=255)
                                                gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                           <  2:3839369/91‑1 (MQ=255)
                                                gAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGAtt                                           <  2:1682575/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)
                                                  aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                         >  1:115873/1‑91 (MQ=255)
                                                   cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:3099512/91‑1 (MQ=255)
                                                   cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  2:1027553/91‑1 (MQ=255)
                                                   cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:3325662/91‑1 (MQ=255)
                                                   cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                        <  1:2280618/91‑1 (MQ=255)
                                                     cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                      <  2:2033507/91‑1 (MQ=255)
                                                      aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                                     >  1:2757447/1‑91 (MQ=255)
                                                       ttAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAt                                    >  1:3799129/1‑91 (MQ=255)
                                                         aaataataaataaACACCTCACCATAAGAAAGGCCTACTTTCTTTTCTGAGGTGTTTATCTGTCGACAACTTGATTATGATTGTTTCAttt                                  >  1:3157127/1‑91 (MQ=255)
                                                           ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                >  1:1933328/1‑91 (MQ=255)
                                                           ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                >  1:357499/1‑91 (MQ=255)
                                                            taataaataaACTCCTCACCATGAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                               <  2:2131532/91‑1 (MQ=255)
                                                            taataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttttt                               <  1:1379458/91‑1 (MQ=255)
                                                              ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                             <  1:1612640/91‑1 (MQ=255)
                                                               taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                            <  2:2626018/91‑1 (MQ=255)
                                                                 aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                          >  1:472804/1‑91 (MQ=255)
                                                                  ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                         <  1:2930751/91‑1 (MQ=255)
                                                                   taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTCTTTGAATacg                        >  2:3620574/1‑89 (MQ=255)
                                                                     aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                      >  1:46993/1‑91 (MQ=255)
                                                                     aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt                      <  2:1963755/91‑1 (MQ=255)
                                                                       cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                    >  2:1581888/1‑91 (MQ=17)
                                                                         ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                  >  2:434273/1‑91 (MQ=255)
                                                                         ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                  <  2:1144837/91‑1 (MQ=255)
                                                                          cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                 >  1:1275346/1‑91 (MQ=255)
                                                                            cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc               >  1:255728/1‑91 (MQ=255)
                                                                              ccATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg             >  1:3578972/1‑91 (MQ=255)
                                                                              ccATAAGAAAGGACTACTTTCTATTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCCTAATCAg             >  2:2408818/1‑91 (MQ=255)
                                                                               caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg             >  1:121043/4‑91 (MQ=21)
                                                                               cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg            >  1:2888643/1‑91 (MQ=255)
                                                                                 tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTTATTTTTTGAATAAGTTCATAATCAGgtg          >  2:1040595/1‑91 (MQ=255)
                                                                                  aaGAAAGTGCTACTTTCTTTTGTGAGGTTTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt         <  2:938597/91‑1 (MQ=255)
                                                                                  aaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg          >  2:3193226/1‑91 (MQ=21)
                                                                                  aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt         >  1:987422/1‑91 (MQ=255)
                                                                                  aaGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtgt         <  1:3371915/91‑1 (MQ=255)
                                                                                   aGAAAGGACTACGTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTa        >  1:3030287/1‑91 (MQ=255)
                                                                                    gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg       <  2:3733675/91‑1 (MQ=255)
                                                                                    gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg       >  2:107104/1‑91 (MQ=255)
                                                                                     aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc      >  2:3008547/1‑91 (MQ=255)
                                                                                     aaaGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGc      >  2:2380536/1‑91 (MQ=255)
                                                                                        ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa   >  2:51291/1‑91 (MQ=255)
                                                                                        ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa   >  1:3088054/1‑91 (MQ=255)
                                                                                         gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc  >  2:2708180/1‑91 (MQ=255)
                                                                                         |                                                                                           
AACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAC  >  USA300TCH1516_ALE/1218641‑1218820

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: