Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,078,152 A→T 13.7% Q510L (CAG→CTG)  purL → Phosphoribosylformylglycinamidine synthase subunit PurL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,078,1520AT13.7% 692.3 / 87.8 248Q510L (CAG→CTG) purLPhosphoribosylformylglycinamidine synthase subunit PurL
Reads supporting (aligned to +/- strand):  ref base A (108/104);  new base T (8/26);  total (117/131)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 251 total aligned reads displayed.

ATCGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATAACGAAAC  >  USA300TCH1516_ALE/1078062‑1078242
                                                                                          |                                                                                          
aTCGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACa                                                                                            >  2:3722515/1‑91 (MQ=255)
  ctaaaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGt                                                                                          <  1:1301786/89‑1 (MQ=255)
   gAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGtt                                                                                         <  2:1534232/91‑1 (MQ=255)
   gAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTAGGGTTCTACTGAAAAGAAAGAAATCTATCAACTGtt                                                                                         <  2:2723104/91‑1 (MQ=255)
    aaaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTg                                                                                        >  2:826470/1‑91 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGa                                                                                       <  2:3554722/91‑1 (MQ=255)
     aaTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGa                                                                                       <  2:1250598/91‑1 (MQ=255)
      agttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGAt                                                                                      <  2:3444258/89‑1 (MQ=255)
      aTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGAt                                                                                      <  1:1881194/91‑1 (MQ=255)
      aTTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGAt                                                                                      >  2:3061077/1‑91 (MQ=255)
      aTTTAATTGCCGTGGGTGCAACACCATTAGCAAAGACCGATTGTTTCAATTATGGTTCTGCGGAAAAGAAAGAAATCTATCAACAGTTGAt                                                                                      >  1:3162944/1‑91 (MQ=255)
       tttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATa                                                                                     <  1:1232838/91‑1 (MQ=255)
        ttAATTGCCGTGGGTGCAGCACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAg                                                                                    <  2:2279959/91‑1 (MQ=255)
        ttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAg                                                                                    <  1:1497159/91‑1 (MQ=255)
        ttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAg                                                                                    <  1:3860731/91‑1 (MQ=255)
        ttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAg                                                                                    <  1:1935503/91‑1 (MQ=255)
        ttAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAg                                                                                    <  2:3592612/91‑1 (MQ=255)
         tAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGa                                                                                   <  2:3254300/91‑1 (MQ=255)
         tAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGa                                                                                   >  1:3693459/1‑91 (MQ=255)
         tAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGa                                                                                   <  1:589699/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAGAGAAAGAAATCTATCAACAGTTGATAGAt                                                                                  <  2:2507863/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGAt                                                                                  <  1:2842277/91‑1 (MQ=255)
          aaTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAt                                                                                  <  2:3800806/91‑1 (MQ=255)
           aTTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGAtt                                                                                 <  1:3012567/91‑1 (MQ=255)
             tGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGAtttc                                                                               <  1:2233853/91‑3 (MQ=255)
             tGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGAtttc                                                                               <  2:3208056/91‑3 (MQ=255)
               ccGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACTGTTGATAGAtttctt                                                                             <  1:2645762/91‑5 (MQ=255)
               ccGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAAc                                                                             >  1:2504415/1‑91 (MQ=255)
               ccGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAAc                                                                             >  1:2519541/1‑91 (MQ=255)
               ccGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAAc                                                                             <  1:730752/91‑1 (MQ=255)
               ccGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTACCAACAGTTGATAGATTCAAc                                                                             <  1:2727170/91‑1 (MQ=255)
                cGTGGGTGCAAGATCATTAGCAATGACAGGTTGTTTAAATTATGGTTCTCCTGTAAAGAAAGAAATCTATCAACAGTTGATATATTCAACg                                                                            <  2:3153396/91‑1 (MQ=255)
                cGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACg                                                                            <  1:1993713/91‑1 (MQ=255)
                cGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACg                                                                            >  1:3558744/1‑91 (MQ=255)
                 gTGGGTGCGACACCATTAGCAATGACAGATTGTGTAAATTATGGTTCTCCTGAGAGGAAAGAAATCTATCAACCGTTGAGAGATTCAACGa                                                                           <  2:1024424/91‑1 (MQ=255)
                 gTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGa                                                                           <  1:3638339/91‑1 (MQ=255)
                  tGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGaa                                                                          <  1:3747381/91‑1 (MQ=255)
                   gggTGCAACACCATTAGCAATGGCAGATTGCTTAATTGATGGTTCTCCTGACAAGAAAGTAATCTATCAACCGTTGATAGATTCAACGaaa                                                                         >  2:1213372/1‑91 (MQ=255)
                   gggTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGaaa                                                                         >  2:3693700/1‑91 (MQ=255)
                     gTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTGAACGAAAgg                                                                       >  1:3547751/1‑91 (MQ=255)
                     gTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAgg                                                                       <  2:487879/91‑1 (MQ=255)
                      tGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGt                                                                      >  2:2021092/1‑91 (MQ=255)
                       gCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTa                                                                     <  1:2009491/91‑1 (MQ=255)
                         aaCACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATg                                                                   >  2:3590706/1‑91 (MQ=255)
                         aaCACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATg                                                                   >  2:2442549/1‑91 (MQ=255)
                          acacCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATgg                                                                  >  2:2337151/1‑91 (MQ=255)
                          acacCATTAGCAATGACAGATCGCTTAAAATATGGTCCTCCTGAAAAGAAAGAAATCTATCAACAGCTCATAGATTCAACGAAAGGTAtcg                                                                  >  2:3264159/1‑89 (MQ=255)
                            acCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCa                                                                <  2:1129263/91‑1 (MQ=255)
                            acCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCa                                                                <  2:3627945/91‑1 (MQ=255)
                             ccATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAg                                                               >  1:273942/1‑91 (MQ=255)
                              cATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGa                                                              >  1:617565/1‑91 (MQ=255)
                                 tAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGc                                                           >  2:784323/1‑91 (MQ=255)
                                 tAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGc                                                           >  2:2166576/1‑91 (MQ=255)
                                  aGCAATGACTGACTGCTTAAATTATGGTTCTCCTGAAAAGAAAGAAATTTATCAACAGCTGATAGATTCAACGAAGGGTATGGCAGAAGCa                                                          <  1:2388245/91‑1 (MQ=255)
                                  aGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCa                                                          >  1:3224572/1‑91 (MQ=255)
                                   gCAATGACAGATTGTTTAAATTATTGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCAt                                                         <  2:3674885/91‑1 (MQ=255)
                                     aaTGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGc                                                       <  1:3911371/91‑1 (MQ=255)
                                      aTGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCg                                                      <  2:3851517/91‑1 (MQ=255)
                                      aTGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCg                                                      <  2:1653116/91‑1 (MQ=255)
                                       tGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGa                                                     <  1:1802547/91‑1 (MQ=255)
                                        gACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGAc                                                    <  2:1969642/91‑1 (MQ=255)
                                         aCAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACa                                                   <  2:2096633/91‑1 (MQ=255)
                                          gaGATGGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGAGAGATTCAACGAAAGGTATGGCAGAAGCATGCGACAt                                                  <  2:1327091/90‑1 (MQ=255)
                                          cAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACAt                                                  >  2:2866882/1‑91 (MQ=255)
                                             aTTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCt                                               >  2:1891797/1‑91 (MQ=255)
                                             aTTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCt                                               >  1:3818238/1‑91 (MQ=255)
                                              ttGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCtt                                              >  1:1051196/1‑91 (MQ=255)
                                                  ttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGa                                          >  1:937151/1‑91 (MQ=255)
                                                  ttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGa                                          >  1:1778840/1‑91 (MQ=255)
                                                  ttAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGa                                          >  1:3256118/1‑91 (MQ=255)
                                                    aaaTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGaca                                        <  2:1221681/91‑1 (MQ=255)
                                                      aTTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACAcc                                      >  2:675857/1‑91 (MQ=255)
                                                        tagggTTCTCCGGAAAAAAAAGAAATCGATCAACAGTTGATAGGTTCAACGAAGGGTATGGCAGAAGCATGCGACATTCTTAAGACACCag                                    <  2:1465691/88‑1 (MQ=255)
                                                        tATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCag                                    <  1:1116641/91‑1 (MQ=255)
                                                        tATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCag                                    <  1:3000710/91‑1 (MQ=255)
                                                         aTGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCagt                                   <  1:2314726/91‑1 (MQ=255)
                                                          tGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCagta                                  >  1:1442866/1‑91 (MQ=255)
                                                           ggTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCagtag                                 >  2:3908964/1‑91 (MQ=255)
                                                           ggTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCagtag                                 >  1:1515785/1‑91 (MQ=255)
                                                             ttCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGtt                               <  1:2098947/91‑1 (MQ=255)
                                                              tctcCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTCAAGACACCAGTAGttt                              >  2:3784883/1‑91 (MQ=255)
                                                               ctcCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTc                             >  1:3805201/1‑91 (MQ=255)
                                                                   tGAGAAGAAAGAAATCGAACAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                         <  2:3893625/91‑1 (MQ=255)
                                                                   tGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGt                         <  1:3370283/91‑1 (MQ=255)
                                                                    gAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTa                        <  1:1634544/91‑1 (MQ=255)
                                                                    gAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTa                        >  2:1164263/1‑91 (MQ=255)
                                                                    gAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTa                        <  1:1369958/91‑1 (MQ=255)
                                                                     aaaaGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTaa                       <  2:2519813/91‑1 (MQ=255)
                                                                        agaaagaaATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGt                    <  2:1777224/91‑1 (MQ=255)
                                                                           aagaaaTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATc                 <  1:2871153/91‑1 (MQ=255)
                                                                                aTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtat            <  2:3418129/91‑1 (MQ=255)
                                                                                aTCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtat            >  1:1194268/1‑91 (MQ=255)
                                                                                   tATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTtatata         >  1:2203772/1‑91 (MQ=255)
                                                                                    aTCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaa        >  2:1457013/1‑91 (MQ=255)
                                                                                    aTCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaa        >  1:2085261/1‑91 (MQ=255)
                                                                                    aTCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaa        <  2:245687/91‑1 (MQ=255)
                                                                                     tCAACTGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACACTCTTAAGACACCAGTAGTTTCTGGTAATGCATCTTTATATaac       >  1:3955386/1‑91 (MQ=255)
                                                                                       aaCAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaacga     >  2:3048562/1‑91 (MQ=255)
                                                                                        aCAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaacgaa    >  1:2608758/1‑91 (MQ=255)
                                                                                          aGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATaacgaaac  <  2:574219/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATCGAAATTTAATTGCCGTGGGTGCAACACCATTAGCAATGACAGATTGTTTAAATTATGGTTCTCCTGAAAAGAAAGAAATCTATCAACAGTTGATAGATTCAACGAAAGGTATGGCAGAAGCATGCGACATTCTTAAGACACCAGTAGTTTCTGGTAATGTATCTTTATATAACGAAAC  >  USA300TCH1516_ALE/1078062‑1078242

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: