Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,155,991 G→C 8.0% intergenic (+167/‑6) mutS2 → / → trxA_2 Endonuclease MutS2/Thioredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,155,9910GC 8.0% 808.8 / 49.0 237intergenic (+167/‑6)mutS2/trxA_2Endonuclease MutS2/Thioredoxin
Reads supporting (aligned to +/- strand):  ref base G (112/106);  new base C (2/17);  total (114/123)
Fisher's exact test for biased strand distribution p-value = 5.30e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 247 total aligned reads displayed.

GAAATAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACATGGT  >  USA300TCH1516_ALE/1155901‑1156081
                                                                                          |                                                                                          
gAAATAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTgg                                                                                            <  1:638927/91‑1 (MQ=255)
gAAATAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTgg                                                                                            <  2:3063646/91‑1 (MQ=255)
gAAATAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTgg                                                                                            >  1:3703020/1‑91 (MQ=255)
 aaaTAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGc                                                                                           >  1:2866345/1‑91 (MQ=255)
   ttAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAt                                                                                         <  2:3843088/90‑1 (MQ=255)
   aTAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCAt                                                                                         <  1:3494663/91‑1 (MQ=255)
    tAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCAtt                                                                                        <  2:1900409/91‑1 (MQ=255)
     aaGCATAAAGATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaa                                                                                       <  2:2672994/91‑3 (MQ=255)
     aaGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCAttt                                                                                       <  2:1372662/91‑1 (MQ=255)
     aaGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCAttt                                                                                       >  1:177566/1‑91 (MQ=255)
     aaGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaa                                                                                       <  1:2180428/91‑3 (MQ=255)
      aGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTa                                                                                      >  1:2420106/1‑91 (MQ=255)
      aGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTa                                                                                      >  1:588792/1‑91 (MQ=255)
      aGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTa                                                                                      <  1:1309009/91‑1 (MQ=255)
      aGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaaa                                                                                      <  1:1914499/91‑4 (MQ=255)
       gCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTAt                                                                                     <  2:193304/91‑1 (MQ=255)
       gCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTAt                                                                                     >  2:2820327/1‑91 (MQ=255)
       gCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaaat                                                                                     <  1:2072839/91‑5 (MQ=255)
         aTAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATgg                                                                                   >  2:3356580/1‑91 (MQ=255)
          tAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGc                                                                                  <  1:2585039/91‑1 (MQ=255)
          tAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGc                                                                                  >  2:1532518/1‑91 (MQ=255)
           aaaaaTTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCa                                                                                 >  1:786949/1‑91 (MQ=255)
                ttGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGt                                                                            >  2:1044051/1‑91 (MQ=255)
                ttGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGt                                                                            <  1:3033404/91‑1 (MQ=255)
                  gAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTaa                                                                          <  1:383380/91‑1 (MQ=255)
                  gAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTaa                                                                          <  2:2093122/91‑1 (MQ=255)
                    gCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTaaaa                                                                        >  1:3554019/1‑91 (MQ=255)
                    gCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTaaaa                                                                        >  2:797165/1‑91 (MQ=255)
                    gCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTaaaa                                                                        >  1:2002623/1‑91 (MQ=255)
                     cAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAg                                                                       <  1:1363431/91‑1 (MQ=255)
                             aaTACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGAt                                                               >  1:2144505/1‑91 (MQ=255)
                              aTACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATg                                                              >  1:2479919/1‑91 (MQ=255)
                              aTACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATg                                                              >  1:3314490/1‑91 (MQ=255)
                               tACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGc                                                             >  1:1627393/1‑91 (MQ=255)
                               tACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaa                                                                                       <  1:1690954/65‑3 (MQ=255)
                               tACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGCCAtaa                                                                                       >  2:1690954/1‑63 (MQ=255)
                                aCATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCa                                                            >  2:2430726/1‑91 (MQ=255)
                                aCATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCa                                                            >  2:1854211/1‑91 (MQ=255)
                                 cATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAg                                                           <  2:887368/91‑1 (MQ=255)
                                 cATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAg                                                           >  1:812896/1‑91 (MQ=255)
                                 cATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAg                                                           <  1:1612793/91‑1 (MQ=255)
                                   tGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGAt                                                         >  2:2297602/1‑91 (MQ=255)
                                     aaGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGAttt                                                       >  2:2074131/1‑91 (MQ=255)
                                     aaGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGAttt                                                       <  1:2391474/91‑1 (MQ=255)
                                      aGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGAtttt                                                      <  1:2407118/91‑1 (MQ=255)
                                       gCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTg                                                     >  1:3436795/1‑91 (MQ=255)
                                       gCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTg                                                     <  1:3576492/91‑1 (MQ=255)
                                       gCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTg                                                     <  2:498074/91‑1 (MQ=255)
                                        cATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGa                                                    >  2:1228485/1‑91 (MQ=255)
                                         aTGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGAt                                                   <  2:1965872/91‑1 (MQ=255)
                                         aTGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGAt                                                   >  2:1200483/1‑91 (MQ=255)
                                           gTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTc                                                 >  2:1092369/1‑91 (MQ=255)
                                             tATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCaa                                               >  1:2076986/1‑91 (MQ=255)
                                             tATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCaa                                               <  2:3877359/91‑1 (MQ=255)
                                             tATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCaa                                               <  1:966759/91‑1 (MQ=255)
                                              aTCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGAt                                                   >  1:1157921/1‑86 (MQ=255)
                                               tCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCaaaa                                             <  2:3779183/91‑1 (MQ=255)
                                                cTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAg                                            >  1:1622799/1‑91 (MQ=255)
                                                 tGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGt                                           >  1:2451433/1‑91 (MQ=255)
                                                  gATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTa                                          <  2:2650388/91‑1 (MQ=255)
                                                  gATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTa                                          <  2:2649636/91‑1 (MQ=255)
                                                  gATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTa                                          <  2:3297540/91‑1 (MQ=255)
                                                  gATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTa                                          <  1:2820327/91‑1 (MQ=255)
                                                  gATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTa                                          >  1:1582167/1‑91 (MQ=255)
                                                   atatAATTTGAACATCATAATAATAATTAAGGAGGATTAGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAg                                         <  1:2791111/91‑1 (MQ=255)
                                                    tataATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGa                                        <  1:721708/91‑1 (MQ=255)
                                                      taATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAAt                                      <  2:1807010/91‑1 (MQ=255)
                                                       aaTTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATc                                     >  1:2750212/1‑91 (MQ=255)
                                                        aTTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCt                                    >  2:269895/1‑91 (MQ=255)
                                                        aTTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCt                                    >  1:929428/1‑91 (MQ=255)
                                                         tttGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTg                                   <  2:2998159/91‑1 (MQ=255)
                                                          ttGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTgg                                  >  1:2652322/1‑91 (MQ=255)
                                                            gAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGgtg                                >  2:1626986/1‑91 (MQ=255)
                                                            gAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGgtg                                >  2:291230/1‑91 (MQ=255)
                                                            gAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGgtg                                >  2:2494877/1‑91 (MQ=255)
                                                              aCATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTa                              >  2:3521803/1‑91 (MQ=255)
                                                               catcatAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTAc                             >  2:479275/1‑91 (MQ=255)
                                                               catcatAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTAc                             <  2:117329/91‑1 (MQ=255)
                                                                atcatAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACa                            <  2:3209967/91‑1 (MQ=255)
                                                                atcatAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACa                            <  1:2515017/91‑1 (MQ=255)
                                                                 tcatAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACaa                           >  2:2346082/1‑91 (MQ=255)
                                                                  catAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGt                                                                      <  2:2126701/47‑1 (MQ=255)
                                                                  catAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAAc                          >  1:2235589/1‑91 (MQ=255)
                                                                  catAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAAc                          >  1:1755131/1‑91 (MQ=255)
                                                                     aataataatTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACtag                       <  2:3486375/91‑1 (MQ=255)
                                                                     aataataatTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACtag                       <  1:1258681/91‑1 (MQ=255)
                                                                     aataataatTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACtag                       <  1:986293/91‑1 (MQ=255)
                                                                      ataataatTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACtagt                      <  1:309358/91‑1 (MQ=255)
                                                                        aataatTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACtagtag                    <  2:2144505/91‑1 (MQ=255)
                                                                             ttAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGAtttt               >  1:3255241/1‑91 (MQ=255)
                                                                              tAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGgtg                                <  1:272943/73‑1 (MQ=255)
                                                                               aaGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTg             <  1:480041/91‑1 (MQ=255)
                                                                               aaGGAGGATTGGCATTTATGGCAATCATAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTg             <  1:2261076/91‑1 (MQ=255)
                                                                                  gaggaTTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGCATCTGGTGTACAACTAGCAGATTCTTGGGc          >  1:1122194/1‑91 (MQ=255)
                                                                                  gaggaTTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGc          <  2:3105987/91‑1 (MQ=255)
                                                                                  gaggaTTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGc          <  2:896936/91‑1 (MQ=255)
                                                                                      aTTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACa      <  2:786949/91‑1 (MQ=255)
                                                                                        tGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACATg    <  1:1491493/91‑1 (MQ=255)
                                                                                         ggCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACATgg   <  1:1933080/91‑1 (MQ=255)
                                                                                          gCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACATGgt  <  2:1352586/91‑1 (MQ=255)
                                                                                          |                                                                                          
GAAATAAGCATAAAAATTGAGCAAATAGAAATACATGAAGCATGTTATCTGATATAATTTGAACATCATAATAATAATTAAGGAGGATTGGCATTTATGGCAATCGTAAAAGTAACAGATGCAGATTTTGATTCAAAAGTAGAATCTGGTGTACAACTAGTAGATTTTTGGGCAACATGGT  >  USA300TCH1516_ALE/1155901‑1156081

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: