Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,362 G→A 8.1% H104Y (CAC→TAC)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,3620GA 8.1% 652.9 / 30.1 207H104Y (CAC→TAC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (89/101);  new base A (7/10);  total (96/111)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

CTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAA  >  USA300TCH1516_ALE/1744272‑1744452
                                                                                          |                                                                                           
cTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTtgtg                                                                                             >  2:1509716/1‑91 (MQ=255)
 tAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGc                                                                                            >  1:3646650/1‑91 (MQ=255)
  aaTTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCa                                                                                           <  1:356660/91‑1 (MQ=255)
    ttttATTGCTGCTTCAATCGTACCCCCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTCTAATTGCATCTTTGTGCAtt                                                                                         >  2:894674/1‑91 (MQ=255)
    ttttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCAtt                                                                                         <  2:2593976/91‑1 (MQ=255)
    ttttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCAtt                                                                                         >  1:193972/1‑91 (MQ=255)
      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                       <  1:871312/91‑1 (MQ=255)
      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                       <  2:1900197/91‑1 (MQ=255)
      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                       <  1:1427544/91‑1 (MQ=255)
      ttATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                       <  1:243665/91‑1 (MQ=255)
       tATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGtt                                                                                      >  2:1725554/1‑91 (MQ=255)
        aTTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTa                                                                                    <  2:2911324/91‑1 (MQ=255)
        aTTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTa                                                                                    <  1:3836617/91‑1 (MQ=255)
         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaa                                                                                   >  2:3233073/1‑91 (MQ=255)
         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaa                                                                                   <  1:3599613/91‑1 (MQ=255)
         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaa                                                                                   >  1:1572663/1‑91 (MQ=255)
         ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaa                                                                                   <  2:339534/91‑1 (MQ=255)
          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaaa                                                                                  <  1:415558/91‑1 (MQ=255)
          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaaa                                                                                  >  2:3896309/1‑91 (MQ=255)
           gctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaaaa                                                                                 >  1:3768562/1‑91 (MQ=255)
           gctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTaaaa                                                                                 >  1:2545246/1‑91 (MQ=255)
              gctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACAt                                                                              >  1:3326786/1‑91 (MQ=255)
                 tCAATCGTACCTCCAGTAGCTAATAAATCATCTTTAATTAACGCGCGTTGGCCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTg                                                                           <  2:2507276/91‑1 (MQ=255)
                  cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGt                                                                          >  2:2378480/1‑91 (MQ=255)
                    aTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTAc                                                                        <  1:1719043/91‑1 (MQ=255)
                     tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTAcc                                                                       >  1:1169012/1‑91 (MQ=255)
                      cGTACCACCAGTAGCTAATAATTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAAACATTTGTAcc                                                                       <  2:1329511/91‑1 (MQ=255)
                      cGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCa                                                                      >  2:604440/1‑91 (MQ=255)
                       gTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCat                                                                     >  1:2750908/1‑91 (MQ=255)
                        tACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCata                                                                    >  1:2500388/1‑91 (MQ=255)
                         accaccAGTAGCTAATGAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCatat                                                                   <  1:2377816/91‑1 (MQ=255)
                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                                >  1:1066535/1‑91 (MQ=255)
                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                                <  1:921801/91‑1 (MQ=255)
                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                                <  2:3072520/91‑1 (MQ=255)
                            accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCt                                                                >  1:1733606/1‑91 (MQ=255)
                             ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTa                                                               >  2:3031508/1‑91 (MQ=255)
                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                          <  1:2786974/91‑1 (MQ=255)
                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                          <  1:2621029/91‑1 (MQ=255)
                                    cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATa                                                        >  2:2167639/1‑91 (MQ=255)
                                     taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                       <  1:1067900/91‑1 (MQ=255)
                                     taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                       <  2:1986506/91‑1 (MQ=255)
                                     taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                       <  2:1364994/91‑1 (MQ=255)
                                      aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACt                                                      <  1:2291801/91‑1 (MQ=255)
                                      aataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACt                                                      <  2:1300665/91‑1 (MQ=255)
                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                     <  1:1156747/91‑1 (MQ=255)
                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                     <  2:3620936/91‑1 (MQ=255)
                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                     <  1:158611/91‑1 (MQ=255)
                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTc                                                     <  2:2314708/91‑1 (MQ=255)
                                        taaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCa                                                    >  2:1287907/1‑91 (MQ=255)
                                        taaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCa                                                    >  1:2845465/1‑91 (MQ=255)
                                         aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCAt                                                   >  1:205419/1‑91 (MQ=255)
                                            tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAAc                                                <  1:771105/91‑1 (MQ=255)
                                             catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACg                                               <  2:3844696/91‑1 (MQ=255)
                                             catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACg                                               <  1:706598/91‑1 (MQ=255)
                                             catcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACCCATAACg                                               >  2:265701/1‑91 (MQ=255)
                                              atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGa                                              >  2:530164/1‑91 (MQ=255)
                                              atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGa                                              >  2:918819/1‑91 (MQ=255)
                                                cTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAAt                                            >  2:3701729/1‑91 (MQ=255)
                                                cTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAAt                                            >  2:277079/1‑91 (MQ=255)
                                                 tGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATg                                           <  2:1047968/91‑1 (MQ=255)
                                                    aaTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACt                                        >  2:2974847/1‑91 (MQ=255)
                                                    aaTTAACACACGTTGACCTGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACtt                                       >  2:367905/1‑91 (MQ=255)
                                                     aTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACtt                                       >  1:2149614/1‑91 (MQ=255)
                                                       tAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCa                                     >  2:2181338/1‑91 (MQ=255)
                                                       tAACACACGTTGACCTGGTTTAACTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCa                                     >  2:1582115/1‑91 (MQ=255)
                                                          cacacGGTGAGGTAGTTTAATTGGATCTTTGTGCAGTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTGCACGa                                  <  2:3590299/91‑1 (MQ=255)
                                                           acacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCGTACTCATAACGAATGACTTCACGAg                                 >  1:3685026/1‑91 (MQ=255)
                                                            cacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                <  1:1572142/91‑1 (MQ=255)
                                                            cacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                >  2:342824/1‑91 (MQ=255)
                                                              cGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTa                              >  1:1344038/1‑91 (MQ=255)
                                                               gTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTaa                             <  1:1711304/91‑1 (MQ=255)
                                                                    ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                        <  1:1683603/91‑1 (MQ=255)
                                                                    ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                        <  2:2417654/91‑1 (MQ=255)
                                                                    ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                        <  2:493510/91‑1 (MQ=255)
                                                                    ccTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                        >  2:1094814/1‑91 (MQ=255)
                                                                      tGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCct                      >  1:3675179/1‑91 (MQ=255)
                                                                       ggTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATACTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCctt                     <  1:3138826/91‑1 (MQ=255)
                                                                        gTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTGGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttc                    <  2:3055715/91‑1 (MQ=255)
                                                                         tttAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAGTTTCCcttct                   >  1:3877856/1‑91 (MQ=255)
                                                                         tttAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttct                   <  2:2692890/91‑1 (MQ=255)
                                                                            aaTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                <  2:2041770/91‑1 (MQ=255)
                                                                             aTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTc               >  1:3464454/1‑91 (MQ=255)
                                                                              ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCt              <  1:962479/91‑1 (MQ=255)
                                                                               tGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa             >  1:2419691/1‑91 (MQ=255)
                                                                               tGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTa             >  1:3576568/1‑91 (MQ=255)
                                                                                gCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTaa            >  1:510779/1‑91 (MQ=255)
                                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc           <  1:662683/91‑1 (MQ=255)
                                                                                 cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc           >  1:117837/1‑91 (MQ=255)
                                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg         <  1:494034/91‑1 (MQ=255)
                                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg         <  2:211909/91‑1 (MQ=255)
                                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg         <  2:1146970/91‑1 (MQ=255)
                                                                                   tCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg         <  1:2665096/91‑1 (MQ=255)
                                                                                   tCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAg         <  2:1399094/91‑1 (MQ=255)
                                                                                       tgtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGc     <  2:1631479/91‑1 (MQ=255)
                                                                                       tgtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCATCTTTTCTAACAGGTGc     >  1:3381618/1‑91 (MQ=255)
                                                                                        gtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTTACAGGTGCa    <  1:1976032/91‑1 (MQ=255)
                                                                                         tgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa   >  1:18875/1‑91 (MQ=255)
                                                                                          gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaaa  >  1:2793407/1‑91 (MQ=255)
                                                                                          gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaaa  >  1:342050/1‑91 (MQ=255)
                                                                                          gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaaa  >  1:2353852/1‑91 (MQ=255)
                                                                                          |                                                                                           
CTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAA  >  USA300TCH1516_ALE/1744272‑1744452

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: