Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,722 T→A 8.9% intergenic (+213/+50) 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,7220TA 8.9% 341.6 / 20.7 101intergenic (+213/+50)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/50);  new base A (2/7);  total (44/57)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 106 total aligned reads displayed.

ATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTG  >  USA300TCH1516_ALE/1218632‑1218812
                                                                                          |                                                                                           
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  1:2477115/91‑1 (MQ=255)
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  1:701513/91‑1 (MQ=255)
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  1:720233/91‑1 (MQ=255)
aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                             <  1:756308/91‑1 (MQ=255)
   aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                          >  1:1882312/1‑91 (MQ=255)
   aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                          >  1:3883554/1‑91 (MQ=255)
        tAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAgg                                                                                     <  1:2377804/91‑1 (MQ=255)
         attttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGa                                                                                    <  1:543140/88‑1 (MQ=255)
         aaCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGa                                                                                    <  1:152962/91‑1 (MQ=255)
          aCTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGac                                                                                   <  1:3892062/91‑1 (MQ=255)
            ttAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGacta                                                                                 >  1:2207049/1‑91 (MQ=255)
             tAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactac                                                                                >  1:2093287/1‑91 (MQ=255)
               aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGTATTAAATAATAAATAAACACCTCACCATACGAAAGGACTActt                                                                              >  1:4147222/1‑91 (MQ=255)
               aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                              >  1:3826636/1‑91 (MQ=255)
               aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                              >  1:2717950/1‑91 (MQ=255)
                 tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                            <  1:421807/91‑1 (MQ=255)
                 tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                            >  1:2047458/1‑91 (MQ=255)
                 tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                            >  1:1210362/1‑91 (MQ=255)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  1:1585444/91‑1 (MQ=25)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  1:1585471/91‑1 (MQ=25)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  1:1601879/91‑1 (MQ=25)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                          <  1:3947360/91‑1 (MQ=25)
                  gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                         <  1:663742/91‑2 (MQ=21)
                    tgatACGTTATTTAAGAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                        <  1:3426939/91‑1 (MQ=255)
                     gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                       >  1:1824839/1‑91 (MQ=255)
                     gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCtttt                                                                       >  1:4019707/1‑91 (MQ=255)
                     gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatg                                                                       <  1:1884895/91‑4 (MQ=21)
                      atACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtg                                                                      <  1:694064/91‑1 (MQ=255)
                       tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                     >  1:3665857/1‑91 (MQ=255)
                       tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                     <  1:985161/91‑1 (MQ=255)
                           ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                 <  1:3816301/91‑1 (MQ=255)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                >  1:3077971/1‑91 (MQ=255)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgt                                                                <  1:2599431/91‑1 (MQ=255)
                            tATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTCTTGTGAGgt                                                                >  1:3962244/1‑91 (MQ=255)
                             aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                               <  1:3770732/91‑1 (MQ=255)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  1:648349/91‑1 (MQ=255)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  1:2254066/91‑1 (MQ=255)
                              tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                              <  1:2254227/91‑1 (MQ=255)
                                   aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                         >  1:201982/1‑91 (MQ=255)
                                    aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                        <  1:3499544/91‑1 (MQ=255)
                                    aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                        >  1:1532662/1‑91 (MQ=255)
                                    aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                        >  1:1127949/1‑91 (MQ=255)
                                    aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                        >  1:3505261/1‑91 (MQ=255)
                                     acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCGCACAAAAGAAAGTAGGCCTTTCTTATGGTGAGGTGTTTAttta                                                       <  1:1734314/91‑2 (MQ=9)
                                       acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                     >  1:1977915/1‑91 (MQ=255)
                                        cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                    >  1:373790/1‑91 (MQ=255)
                                         gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                   <  1:1231524/91‑1 (MQ=255)
                                            tcTTCGTTTCGATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                <  1:2457997/90‑1 (MQ=255)
                                            aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                <  1:2458255/91‑1 (MQ=255)
                                               tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             >  1:1343173/1‑91 (MQ=255)
                                               tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             >  1:593848/1‑91 (MQ=255)
                                               tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             >  1:2174579/1‑91 (MQ=255)
                                               tCGATTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                             <  1:3732165/91‑1 (MQ=255)
                                                cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaatcat                                            <  1:1713389/91‑5 (MQ=255)
                                                       aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                     >  1:1879615/1‑91 (MQ=255)
                                                       aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                     >  1:3203564/1‑91 (MQ=255)
                                                          aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                  <  1:2064164/91‑1 (MQ=255)
                                                          aaCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTg                                  <  1:3470896/91‑1 (MQ=255)
                                                           aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                 <  1:1320388/91‑1 (MQ=255)
                                                           aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                 <  1:3898920/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  1:1011621/91‑1 (MQ=255)
                                                            cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                <  1:3542953/91‑1 (MQ=255)
                                                             gCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGttt                               <  1:2694856/91‑1 (MQ=255)
                                                              cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                              <  1:123855/91‑1 (MQ=255)
                                                              cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                              >  1:4032216/1‑91 (MQ=255)
                                                               aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCa                             >  1:3787701/1‑91 (MQ=255)
                                                                ttAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAt                            <  1:3947531/91‑1 (MQ=255)
                                                                 taaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtt                           >  1:3094013/1‑91 (MQ=255)
                                                                  aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                          >  1:3359560/1‑91 (MQ=255)
                                                                  aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                          >  1:1290563/1‑91 (MQ=255)
                                                                  aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                          <  1:3310153/91‑1 (MQ=255)
                                                                      aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                      <  1:4101985/91‑1 (MQ=255)
                                                                      aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                      <  1:2646079/91‑1 (MQ=255)
                                                                       ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                     <  1:1028794/91‑1 (MQ=255)
                                                                       ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                     <  1:2529880/91‑1 (MQ=255)
                                                                          aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa                  <  1:1340590/91‑1 (MQ=255)
                                                                           ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                 <  1:2315416/91‑1 (MQ=255)
                                                                           ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa                 <  1:2315104/91‑1 (MQ=255)
                                                                            taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                <  1:97389/91‑1 (MQ=255)
                                                                            taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                <  1:2447652/91‑1 (MQ=255)
                                                                            taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                >  1:3971644/1‑91 (MQ=255)
                                                                            taaaCACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg                >  1:3973087/1‑91 (MQ=14)
                                                                              aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt              <  1:3588558/91‑1 (MQ=255)
                                                                              aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt              <  1:1039404/91‑1 (MQ=255)
                                                                              aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt              <  1:3022825/91‑1 (MQ=255)
                                                                               acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc             <  1:444594/91‑1 (MQ=255)
                                                                               acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc             <  1:444593/91‑1 (MQ=255)
                                                                                cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa            >  1:2754694/1‑91 (MQ=255)
                                                                                cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa            >  1:2754869/1‑91 (MQ=255)
                                                                                cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa            >  1:1492270/1‑91 (MQ=17)
                                                                                 acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt           >  1:1643654/1‑91 (MQ=255)
                                                                                 acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt           >  1:1643653/1‑91 (MQ=255)
                                                                                   cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa         >  1:34598/1‑91 (MQ=255)
                                                                                    tCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAAt        >  1:102284/1‑91 (MQ=255)
                                                                                     cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATGAGTTCATAATc       <  1:659202/91‑1 (MQ=255)
                                                                                     cACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATc       <  1:2890507/91‑1 (MQ=255)
                                                                                        cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg    >  1:1513701/1‑91 (MQ=255)
                                                                                        cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg    >  1:1166397/1‑91 (MQ=255)
                                                                                         aTAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgt   >  1:2021047/1‑91 (MQ=255)
                                                                                          tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg  >  1:3148655/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTG  >  USA300TCH1516_ALE/1218632‑1218812

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