Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,719 Δ1 bp 8.8% intergenic (+210/+53) 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,7190C. 8.8% 348.3 / 26.9 102intergenic (+210/+53)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/54);  new base . (2/7);  total (41/61)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 105 total aligned reads displayed.

AAGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATC  >  USA300TCH1516_ALE/1218629‑1218807
                                                                                          |                                                                                         
aaGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCAc                                                                                           <  1:2534916/91‑1 (MQ=255)
aaGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCAc                                                                                           >  1:1546056/1‑91 (MQ=255)
aaGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCAc                                                                                           <  1:1349248/91‑1 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                        <  1:2477115/91‑1 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                        <  1:756308/91‑1 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                        <  1:720233/91‑1 (MQ=255)
   aTAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCAt                                                                                        <  1:701513/91‑1 (MQ=255)
      aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                     >  1:1882312/1‑91 (MQ=255)
      aTTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAg                                                                                     >  1:3883554/1‑91 (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:2717950/1‑91 (MQ=255)
                  aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                         >  1:3826636/1‑91 (MQ=255)
                    tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                       >  1:2047458/1‑91 (MQ=255)
                    tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                       >  1:1210362/1‑91 (MQ=255)
                    tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                       <  1:421807/91‑1 (MQ=255)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                     <  1:1585471/91‑1 (MQ=25)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                     <  1:1585444/91‑1 (MQ=25)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttctt                                                                     <  1:3947360/91‑1 (MQ=25)
                     gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATATAAATAAACACCTCACAAAAGAAAGTAGTCctttctta                                                                    <  1:663742/91‑2 (MQ=21)
                        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:2254066/91‑1 (MQ=255)
                                 tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                         <  1:2254227/91‑1 (MQ=255)
                                 tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                         <  1:648349/91‑1 (MQ=255)
                                      aaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAtt                                                    >  1:201982/1‑91 (MQ=255)
                                       aaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTAttt                                                   >  1:3505261/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)
                                        acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTg                                                  <  1:2737147/91‑1 (MQ=255)
                                        acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCGCACAAAAGAAAGTAGGCCTTTCTTATGGTGAGGTGTTTAttta                                                  <  1:1734314/91‑2 (MQ=9)
                                          acGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTc                                                >  1:1977915/1‑91 (MQ=255)
                                           cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                               >  1:373790/1‑91 (MQ=255)
                                           cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                               <  1:2098546/91‑1 (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:593848/1‑91 (MQ=255)
                                                  tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                        >  1:1343173/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:3203564/1‑91 (MQ=255)
                                                          aTGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGa                                >  1:1879615/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:3542953/91‑1 (MQ=255)
                                                               cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                           <  1:1011621/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:3310153/91‑1 (MQ=255)
                                                                     aaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttt                     >  1:1290563/1‑91 (MQ=255)
                                                                         aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                 <  1:2646079/91‑1 (MQ=255)
                                                                         aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                 <  1:4101985/91‑1 (MQ=255)
                                                                          ataaatagACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                >  1:3616726/1‑91 (MQ=255)
                                                                          ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                <  1:1028794/91‑1 (MQ=255)
                                                                          ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                <  1:2529880/91‑1 (MQ=255)
                                                                           taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa               >  1:4050174/1‑91 (MQ=255)
                                                                             aataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATa             <  1:1340590/91‑1 (MQ=255)
                                                                              ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa            <  1:2315104/91‑1 (MQ=255)
                                                                              ataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATaa            <  1:2315416/91‑1 (MQ=255)
                                                                               taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg           <  1:2447652/91‑1 (MQ=255)
                                                                               taaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg           <  1:97389/91‑1 (MQ=255)
                                                                               taaaCACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAg           >  1:3973087/1‑91 (MQ=14)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt         <  1:1544599/91‑1 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt         <  1:3022825/91‑1 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt         <  1:1039404/91‑1 (MQ=255)
                                                                                 aaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGtt         <  1:3588558/91‑1 (MQ=255)
                                                                                  acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc        <  1:444593/91‑1 (MQ=255)
                                                                                  acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc        <  1:444594/91‑1 (MQ=255)
                                                                                   cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa       >  1:2754869/1‑91 (MQ=255)
                                                                                   cacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa       >  1:2754694/1‑91 (MQ=255)
                                                                                   cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa       >  1:1492270/1‑91 (MQ=17)
                                                                                    acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt      >  1:1643653/1‑91 (MQ=255)
                                                                                    acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt      >  1:1643654/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)
                                                                                          |                                                                                         
AAGATAATTTGTAACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATC  >  USA300TCH1516_ALE/1218629‑1218807

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: