Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 24,137 G→A 5.9% intergenic (+389/‑41) purA → / → USA300TCH1516_00018 Adenylosuccinate synthetase/tRNA‑Glu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE24,1370GA 5.9% 815.9 / 25.8 254intergenic (+389/‑41)purA/USA300TCH1516_00018Adenylosuccinate synthetase/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base G (102/136);  new base A (12/3);  total (114/140)
Fisher's exact test for biased strand distribution p-value = 6.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 267 total aligned reads displayed.

AATGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAACACG  >  USA300TCH1516_ALE/24047‑24227
                                                                                          |                                                                                          
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                            <  2:510522/91‑1 (MQ=255)
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                            <  2:355333/91‑1 (MQ=255)
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                            <  2:2437099/91‑1 (MQ=255)
aaTGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATgg                                                                                            <  2:1180355/91‑1 (MQ=255)
  tGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCa                                                                                          <  2:1128392/91‑1 (MQ=255)
   gAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCaa                                                                                         <  2:1465939/91‑1 (MQ=255)
    aaTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAg                                                                                        <  2:1418177/91‑1 (MQ=255)
     aTAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAAGAATAGCAAAATTATGTAAAAAATATCTTGCCATGACAAGa                                                                                       <  2:3187739/91‑1 (MQ=255)
      tAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGAt                                                                                      <  2:3735307/91‑1 (MQ=255)
      tAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGAt                                                                                      >  2:1328522/1‑91 (MQ=255)
         accgCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTgg                                                                                   >  2:2252993/4‑91 (MQ=255)
           tGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCt                                                                                 >  1:2167387/1‑91 (MQ=255)
             ccAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgt                                                                               <  2:1743369/91‑1 (MQ=255)
             ccAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgt                                                                               >  2:3033027/1‑91 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                              <  2:2651159/91‑1 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                              >  1:350672/1‑91 (MQ=255)
              cAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCtgtg                                                                              >  1:349462/1‑91 (MQ=255)
                aTAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCt                                                                            <  2:195468/91‑1 (MQ=255)
                aTAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCt                                                                            >  1:657805/1‑91 (MQ=255)
                  aaaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtat                                                                          <  2:948156/91‑1 (MQ=255)
                   aaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtata                                                                         <  2:3821159/91‑1 (MQ=255)
                   aaaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCtata                                                                         >  1:2113066/1‑91 (MQ=255)
                    aaTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATaa                                                                        <  1:512850/91‑1 (MQ=255)
                     aTGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAat                                                                       >  1:2375755/1‑91 (MQ=255)
                      tGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatc                                                                      >  2:3031290/1‑91 (MQ=255)
                      tGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatc                                                                      >  2:2538665/1‑91 (MQ=255)
                       gCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatct                                                                     <  1:1984965/91‑1 (MQ=255)
                        cAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatcta                                                                    <  1:1747875/91‑1 (MQ=255)
                        cAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatcta                                                                    >  1:793380/1‑91 (MQ=255)
                         aaTTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctat                                                                   >  2:450680/1‑91 (MQ=255)
                          aTTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAatctatc                                                                  >  2:2871659/1‑91 (MQ=255)
                              ttCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTtgt                                                              >  1:1009349/1‑91 (MQ=255)
                                cACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGc                                                            >  2:729880/1‑91 (MQ=255)
                                  cTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCtt                                                          >  1:2447647/1‑91 (MQ=255)
                                     ataACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAg                                                       <  2:2236320/91‑1 (MQ=255)
                                     ataACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAg                                                       <  2:192081/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                      <  1:330492/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                      <  1:1517233/91‑1 (MQ=255)
                                      taACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGa                                                      <  1:3029783/91‑1 (MQ=255)
                                       aaCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTGTCTTGTGCTTAAGaa                                                     <  1:3314398/91‑1 (MQ=255)
                                       aaCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGaa                                                     >  2:750266/1‑91 (MQ=255)
                                        aCATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAAc                                                    <  1:1167452/91‑1 (MQ=255)
                                        aCATTCTTCAAAAAATAATAGCAAAATTATATAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAAc                                                    <  2:396721/91‑1 (MQ=255)
                                         cATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTTATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACg                                                   >  2:2713609/1‑91 (MQ=255)
                                         cATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACg                                                   >  2:2399695/1‑91 (MQ=255)
                                         cATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACg                                                   <  2:2224919/91‑1 (MQ=255)
                                           ttcttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGc                                                 <  1:2924903/91‑1 (MQ=255)
                                           ttcttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGc                                                 >  2:3403655/1‑91 (MQ=255)
                                           ttcttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGc                                                 <  1:353538/91‑1 (MQ=255)
                                            tcttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCt                                                >  1:1430623/1‑91 (MQ=255)
                                             cttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTc                                               >  2:625715/1‑91 (MQ=255)
                                              ttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTcc                                              <  1:3719532/91‑1 (MQ=255)
                                              ttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTcc                                              <  1:2723478/91‑1 (MQ=255)
                                              ttcAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTcc                                              <  1:3096126/91‑1 (MQ=255)
                                                cAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCtt                                            >  1:95173/1‑91 (MQ=255)
                                                cAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCtt                                            >  2:3667112/1‑91 (MQ=255)
                                                  aaaaaTAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTgg                                          >  1:2063651/1‑91 (MQ=255)
                                                   aaaaTAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGt                                         >  2:657002/1‑91 (MQ=255)
                                                   aaaaTAATAGCAAAATTATGTAAAAAATATCGTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGt                                         <  2:3251786/91‑1 (MQ=255)
                                                    aaaTAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTc                                        >  2:772217/1‑91 (MQ=255)
                                                       taataGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAg                                     >  1:1903665/1‑91 (MQ=255)
                                                       taataGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAg                                     >  1:3864379/1‑91 (MQ=255)
                                                       taagaGCAAAATTATGGAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAg                                     <  1:3564060/91‑1 (MQ=255)
                                                         ataGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCg                                   <  2:1531970/91‑1 (MQ=255)
                                                          taGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCgg                                  <  2:1861677/91‑1 (MQ=255)
                                                           aGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGt                                 >  2:1271776/1‑91 (MQ=255)
                                                           aGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGt                                 >  1:59003/1‑91 (MQ=255)
                                                             cAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGtgt                               <  1:1273164/91‑3 (MQ=255)
                                                                  ttATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGaca                          <  1:1310165/91‑1 (MQ=255)
                                                                  ttATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGaca                          <  1:2752193/91‑1 (MQ=255)
                                                                  ttATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGaca                          <  2:328725/91‑1 (MQ=255)
                                                                     tGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCg                       <  2:2451838/91‑1 (MQ=255)
                                                                      gTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGc                      >  1:2681823/1‑91 (MQ=255)
                                                                      gTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTAGTGCTTAAGAACGGCGCCTTGGTCAAGCGGTTAAGACACCGc                      >  2:1372455/1‑91 (MQ=255)
                                                                       tAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGcc                     <  2:209636/91‑1 (MQ=255)
                                                                       tAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGcc                     <  2:1439249/91‑1 (MQ=255)
                                                                       tAAAAAATACCTTGTCATGGCAAGATTGGCTGTGCTATAATCTAGCTTGTGCTTAAGACCGGCTCCTTGGTCAAGCGGTTAAGACTCCGcc                     <  2:2038014/91‑1 (MQ=38)
                                                                        aataaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccc                    <  2:1755532/88‑1 (MQ=255)
                                                                        aaaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccc                    <  2:1576662/91‑1 (MQ=255)
                                                                        aaaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccc                    >  1:1129208/1‑91 (MQ=255)
                                                                         aaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCt                   <  1:2300530/91‑1 (MQ=255)
                                                                         aaaaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCt                   <  1:3293441/91‑1 (MQ=255)
                                                                           aaaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGccgtga                 <  1:929996/91‑5 (MQ=255)
                                                                            aaTATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTc                <  2:3348630/91‑1 (MQ=255)
                                                                             atatCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCa               <  1:1245906/91‑1 (MQ=255)
                                                                              tatCTTGTCATGGCAAGATTGGCTGTGCTATAATATATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAc              >  1:1621985/1‑91 (MQ=255)
                                                                                tCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcgg            >  2:2137881/1‑91 (MQ=255)
                                                                                 cTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggc           >  1:1483460/1‑91 (MQ=255)
                                                                                  ttgccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcg          >  1:203857/5‑91 (MQ=255)
                                                                                   tgccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcgg         >  1:906385/4‑91 (MQ=255)
                                                                                   tGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCAcggcgg         >  1:3152906/1‑91 (MQ=255)
                                                                                    gTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGt        >  1:2505275/1‑91 (MQ=255)
                                                                                    gTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGt        >  2:894934/1‑91 (MQ=255)
                                                                                     ccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTa       >  2:1014220/2‑91 (MQ=255)
                                                                                     ccATGACAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTa       >  2:1655784/2‑91 (MQ=255)
                                                                                      cATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTaa      <  1:3401451/91‑1 (MQ=255)
                                                                                      cATGACAAGATTGGCTGTGCTATCATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGACCTTTCACTGCGGTaa      >  2:3440074/1‑91 (MQ=25)
                                                                                       aTGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAac     <  1:2044803/91‑1 (MQ=255)
                                                                                       aTGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAac     >  1:3794633/1‑91 (MQ=255)
                                                                                          gCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAACACg  >  1:1818166/1‑91 (MQ=255)
                                                                                          |                                                                                          
AATGAATAGACTGCCAATAAAATGCAATTTTTCACTTATAACATTCTTCAAAAAATAATAGCAAAATTATGTAAAAAATATCTTGTCATGGCAAGATTGGCTGTGCTATAATCTATCTTGTGCTTAAGAACGGCTCCTTGGTCAAGCGGTTAAGACACCGCCCTTTCACGGCGGTAACACG  >  USA300TCH1516_ALE/24047‑24227

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