Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,218,732 C→G 7.5% intergenic (+223/+40) 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,7320CG 7.5% 775.2 / 39.1 227intergenic (+223/+40)USA300HOU_RS06065/USA300HOU_RS06070hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (117/93);  new base G (7/10);  total (124/103)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 246 total aligned reads displayed.

ACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAA  >  USA300TCH1516_ALE/1218642‑1218822
                                                                                          |                                                                                           
 cTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTATCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGact                                                                                            >  2:3198227/1‑91 (MQ=255)
   tAAATGATGAAACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactac                                                                                          <  2:2674063/91‑1 (MQ=255)
    aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                         <  1:489917/91‑1 (MQ=255)
    aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                         <  1:3344668/91‑1 (MQ=255)
    aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                         <  1:2908845/91‑1 (MQ=255)
    aaaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                         <  1:2098722/91‑1 (MQ=255)
     aaTGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActt                                                                                        <  1:1060237/91‑1 (MQ=255)
      atgatgATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActtt                                                                                       <  1:691039/91‑1 (MQ=255)
       tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                      <  1:581553/91‑1 (MQ=255)
       tgatgaTACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttc                                                                                      <  1:2551420/91‑1 (MQ=255)
        gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                     >  2:635964/1‑91 (MQ=255)
        gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                     >  2:2401108/1‑91 (MQ=255)
        gatgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttct                                                                                     >  2:1458130/1‑91 (MQ=255)
         atgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttctt                                                                                    >  2:941441/1‑91 (MQ=255)
          tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTACCTAATAAATATACACCTCACCATAAGAAAGGACTActttcttt                                                                                  <  2:116995/91‑1 (MQ=255)
          tgatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTActttcttt                                                                                  >  2:649801/1‑91 (MQ=255)
           gatACGTTATTTAGAGACACGAGACTTCGTGTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAATGTAGTCCTTTCTTAtg                                                                                <  2:2803184/91‑1 (MQ=12)
           gatACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCctttcttatg                                                                                 <  1:2575593/91‑4 (MQ=21)
             tACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgt                                                                               >  2:2606153/1‑91 (MQ=255)
              aCGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                              >  1:2196706/1‑91 (MQ=255)
              aCGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTtgtg                                                                              >  1:1688173/1‑91 (MQ=255)
               cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                             >  2:2064695/1‑91 (MQ=255)
               cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                             <  2:2227039/91‑1 (MQ=255)
               cGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGa                                                                             <  1:2876537/91‑1 (MQ=255)
                gTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAg                                                                            >  1:1755352/1‑91 (MQ=255)
                 ttATTTAAAAGCGCTATACTTCGTGTCAATGAACGCATTAAATCATAAATGAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAgg                                                                           <  2:3603414/91‑1 (MQ=255)
                 ttATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCATTTGTGAgg                                                                           >  2:3502686/1‑91 (MQ=255)
                   aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                         >  1:656775/1‑91 (MQ=255)
                   aTTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtg                                                                         <  1:2605515/91‑1 (MQ=255)
                    tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGactact                                                                                         <  1:480691/75‑1 (MQ=255)
                    tttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGgtgt                                                                        <  1:2671994/91‑1 (MQ=255)
                     ttAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGtt                                                                       <  2:1236386/91‑1 (MQ=14)
                       aaaaaCACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTa                                                                     >  2:1486395/1‑91 (MQ=255)
                           acacGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTAttta                                                                 <  1:3472185/91‑2 (MQ=11)
                              cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                              >  2:1187627/1‑91 (MQ=255)
                              cGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCg                                                              <  1:281834/91‑1 (MQ=255)
                               gATACTTCGTTTCAATGAACGCATTAAATAATAGATATACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                             >  1:3098116/1‑91 (MQ=255)
                               gATACTTCGTTTCAATGAACGCATTAAATAATAAATGAACACCTCACCATGCGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                             <  2:1701818/91‑1 (MQ=255)
                               gATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGa                                                             >  1:2249437/1‑91 (MQ=255)
                                aTACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGAc                                                            >  1:1354463/1‑91 (MQ=255)
                                  aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                          >  2:1992757/1‑91 (MQ=255)
                                  aCTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACaa                                                          >  1:751503/1‑91 (MQ=255)
                                    ttCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACt                                                        <  2:560727/91‑1 (MQ=255)
                                     tCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACtt                                                       >  1:1002745/1‑91 (MQ=255)
                                      cGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTg                                                      >  2:1495819/1‑91 (MQ=255)
                                        tttCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGAt                                                    >  2:325533/1‑91 (MQ=255)
                                           cAATGAACGCATTCAAGAATAAATGACCACCTCACCAGATGAGAGGACTAGTTTCTTTTGTGAGGTGTTTATTTGTCGGCAACTTtattat                                                 <  2:3175290/91‑1 (MQ=255)
                                           cAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTAt                                                 >  1:1617134/1‑91 (MQ=255)
                                                 aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                           <  1:3275020/91‑1 (MQ=255)
                                                 aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                           <  1:982641/91‑1 (MQ=255)
                                                 aCGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGt                                           >  2:2857027/1‑91 (MQ=255)
                                                  cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:3198227/91‑1 (MQ=255)
                                                  cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:1110520/91‑1 (MQ=255)
                                                  cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:3569421/91‑1 (MQ=255)
                                                  cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  1:1458130/91‑1 (MQ=255)
                                                  cGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGtt                                          <  2:2973184/91‑1 (MQ=255)
                                                   gCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGttt                                         <  2:3799854/91‑1 (MQ=255)
                                                    cATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTc                                        <  2:3637281/91‑1 (MQ=255)
                                                     aTTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCATCTGTGAGGAGTTTATTAGTCGACAACTTGATTATGATTGTTTCa                                       >  2:1811735/1‑91 (MQ=255)
                                                      ttAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAt                                      >  2:2297962/1‑91 (MQ=255)
                                                       taaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtt                                     >  2:142639/1‑91 (MQ=255)
                                                       taaataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtt                                     >  1:599926/1‑91 (MQ=255)
                                                         aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                   <  1:719202/91‑1 (MQ=255)
                                                         aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                   <  1:1415530/91‑1 (MQ=255)
                                                         aataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAtttt                                   >  1:1615408/1‑91 (MQ=255)
                                                          ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGttt                                         >  1:3288891/1‑84 (MQ=255)
                                                          ataataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCAttttt                                  >  2:2682939/1‑91 (MQ=255)
                                                            aataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTg                                <  1:2064695/91‑1 (MQ=255)
                                                             ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                               <  2:2947183/91‑1 (MQ=255)
                                                             ataaataaACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGa                               >  2:1733639/1‑91 (MQ=255)
                                                              taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                              >  1:859416/1‑91 (MQ=255)
                                                              taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                              >  2:3040105/1‑91 (MQ=255)
                                                              taaataaaCACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGaa                              >  2:2569645/1‑91 (MQ=255)
                                                                     acacCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTc                       >  1:2761356/1‑91 (MQ=255)
                                                                      cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                      >  1:2134192/1‑91 (MQ=17)
                                                                      cacCTCACAAAAGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCa                      >  2:1542798/1‑91 (MQ=17)
                                                                       acCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCAt                     >  2:478306/1‑91 (MQ=255)
                                                                        ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                    >  1:3734566/1‑91 (MQ=255)
                                                                        ccTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATa                    <  1:3836432/91‑1 (MQ=255)
                                                                         cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                   >  2:3166741/1‑91 (MQ=255)
                                                                         cTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATaa                   <  2:2184102/91‑1 (MQ=255)
                                                                            aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                >  2:3876511/1‑91 (MQ=255)
                                                                            aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                >  1:266435/1‑91 (MQ=255)
                                                                            aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                <  1:956315/91‑1 (MQ=255)
                                                                            aCCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCa                <  2:2244574/91‑1 (MQ=255)
                                                                              caaaaGAAAGTAGTCCTTTCTTATGGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAg               >  2:3729774/4‑91 (MQ=21)
                                                                              cATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAgg              <  2:878955/91‑1 (MQ=255)
                                                                               aTAAGAGAGGACTACTTTCTTTTGTGAGGTGTTTATCTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgt             <  2:3835192/91‑1 (MQ=255)
                                                                                tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg            >  1:1241385/1‑91 (MQ=255)
                                                                                tAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGgtg            <  2:577747/91‑1 (MQ=255)
                                                                                  aGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTa          <  2:3750531/91‑1 (MQ=255)
                                                                                  aGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTa          >  1:3211743/1‑91 (MQ=255)
                                                                                   gAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAg         >  2:2158132/1‑91 (MQ=255)
                                                                                       ggACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATa     >  1:2956183/1‑91 (MQ=255)
                                                                                        gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc    >  2:3693712/1‑91 (MQ=255)
                                                                                        gACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATAc    <  1:2857027/91‑1 (MQ=255)
                                                                                         actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa   >  2:2026196/1‑91 (MQ=255)
                                                                                         actactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACa   <  1:2401108/91‑1 (MQ=255)
                                                                                         actactTTCTTTTGTGAGGTGTTATTTGTCGACAAGTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa  <  2:1052960/91‑1 (MQ=255)
                                                                                          ctactTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACaa  >  2:497288/1‑91 (MQ=255)
                                                                                          |                                                                                           
ACTTAAATGATGATACGTTATTTAAAAACACGATACTTCGTTTCAATGAACGCATTAAATAATAAATAAACACCTCACCATAAGAAAGGACTACTTTCTTTTGTGAGGTGTTTATTTGTCGACAACTTGATTATGATTGTTTCATTTTTTGAATAAGTTCATAATCAGGTGTAGCATACAA  >  USA300TCH1516_ALE/1218642‑1218822

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