Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,323,313 A→G 6.3% K29E (AAA→GAA)  mutL → DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,323,3130AG 6.3% 219.0 / 2.6 74K29E (AAA→GAA) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (23/46);  new base G (3/2);  total (26/48)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCATTGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTATTTCATA  >  USA300TCH1516_ALE/1323173‑1323451
                                                                                                                                            |                                                                                                                                          
agtcATTGGTGAAGGTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAACTAGTTGAAAGACCGAGTTCTGTTGTGa                                                                                                                                            <  2:186731/138‑1 (MQ=255)
aGGCATTGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGa                                                                                                                                            <  1:17499/141‑1 (MQ=255)
 ggCATTGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGGGTTCTGTTGTGaa                                                                                                                                           >  1:376096/1‑141 (MQ=255)
 ggCATTGGTGAAGTTAAGTGAATTACAAAATCAATGAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGaa                                                                                                                                           <  1:168369/141‑1 (MQ=255)
 ggCAGTGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGaa                                                                                                                                           <  1:203054/141‑1 (MQ=255)
     ttGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGACAGACCGAGTTCTGTTGTGAAAGaa                                                                                                                                       <  1:807693/141‑1 (MQ=255)
     ttGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGaa                                                                                                                                       <  1:896951/141‑1 (MQ=255)
     ttGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGaa                                                                                                                                       <  2:461229/141‑1 (MQ=255)
      tGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAAc                                                                                                                                      >  2:1089172/1‑141 (MQ=255)
      tGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAAc                                                                                                                                      >  1:658844/1‑141 (MQ=255)
         tgactttAAGTGAATTCAAAAAAAGATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACCCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGt                                                                                                                                   <  1:469692/136‑1 (MQ=255)
         tGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGt                                                                                                                                   <  1:283605/141‑1 (MQ=255)
         ggAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGt                                                                                                                                   <  1:919352/140‑1 (MQ=255)
          gAAGTTAAGTGAATTACAAGATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGtt                                                                                                                                  <  2:671370/141‑1 (MQ=255)
          gAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAATTCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGtt                                                                                                                                  <  2:947474/141‑1 (MQ=255)
          gAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGtt                                                                                                                                  <  2:248465/141‑1 (MQ=255)
           aaGTTAAGTGAATTACAAAATCAATTAAAAAAAAAAAAAAAAAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTg                                                                                                                                 <  2:357052/141‑1 (MQ=255)
           aaGTTAAGGGAATGACAAAATCAAATAAACTAGAGGTGTTGCAACCTGGTGATAATGAAAGAGCTCCCAACCTCATTAGCAAATAAAACCGCAGCAGGTGATGTAGTTGAAAGACCTAGTTCTGTTGTGAAAGAACTGTTg                                                                                                                                 <  2:547521/141‑1 (MQ=255)
            aggtAAGGGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTgg                                                                                                                                <  2:348950/138‑1 (MQ=255)
            aGTTAAGTGAATTACAAAATCAATTAAAATTGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTgg                                                                                                                                <  1:277774/141‑1 (MQ=255)
            aGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTgg                                                                                                                                <  2:879786/141‑1 (MQ=255)
            aGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTgg                                                                                                                                <  2:978271/141‑1 (MQ=255)
             ggtAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGa                                                                                                                               <  2:322191/139‑1 (MQ=255)
             gTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGa                                                                                                                               <  2:568295/141‑1 (MQ=255)
                aaGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGaaaa                                                                                                                            >  2:494064/1‑141 (MQ=255)
                aaGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACGGTTGGaaaa                                                                                                                            >  1:240006/1‑141 (MQ=255)
                aaGAGAATTACAAAATCAATTAAAATAAAGGTGATGCAGAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATTAAATCGCAGCAGGTGAAGTAGTTGAAAAACCGAGTTCTGTTGTGAAAGACGTGTTGGaaaa                                                                                                                            >  2:98314/1‑141 (MQ=255)
                   tGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGc                                                                                                                         <  1:772809/141‑1 (MQ=255)
                            aaaTCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTCGGAAAAAGCTAAAAATgg                                                                                                                >  2:211835/1‑140 (MQ=255)
                                 aaTTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGgcg                                                                                                           <  2:285458/141‑1 (MQ=255)
                                       aaTAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAg                                                                                                     >  1:1013379/1‑141 (MQ=255)
                                       aaTAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAg                                                                                                     >  1:1033922/1‑141 (MQ=255)
                                       aaTAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAg                                                                                                     <  1:541991/141‑1 (MQ=255)
                                           gagGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAAt                                                                                                 <  2:897155/141‑1 (MQ=255)
                                             ggTGTTGAAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATaa                                                                                               <  2:894884/141‑1 (MQ=255)
                                                gtTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCa                                                                                            >  1:607927/1‑141 (MQ=255)
                                                       aaTGGGGAAAATTAAAGCACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGt                                                                                     >  1:265758/1‑141 (MQ=255)
                                                        aTGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTa                                                                                    >  1:541286/1‑141 (MQ=255)
                                                            ggAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGGAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAgg                                                                                >  2:129638/1‑141 (MQ=255)
                                                            ggAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGGTGAAAGACCGAGTTCTGTTGTGGAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAgg                                                                                >  2:328086/1‑141 (MQ=255)
                                                               aacttAAAGCACTCCGAACCCCATTAGCAAATAGAATCTCAGCATGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAAt                                                                             <  2:484247/138‑1 (MQ=255)
                                                                   tAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTgg                                                                         <  2:964585/141‑1 (MQ=255)
                                                                     aaGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCg                                                                       >  1:1019788/1‑141 (MQ=255)
                                                                     aaGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCg                                                                       >  1:975591/1‑141 (MQ=255)
                                                                            ccAAGCCTCATTACCAAATAAAATCGAAGCAGGGGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGACATAAGAATTGAAGTAGAGGCATCTGGCGTCCAATc                                                                >  2:860123/1‑141 (MQ=255)
                                                                             cAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCt                                                               >  2:298124/1‑141 (MQ=255)
                                                                             cAAACCTCATTAGCAAATAAAAGCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGAGTCCAAtgt                                                               >  1:118884/1‑139 (MQ=255)
                                                                               aaCCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACGGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAAAAAGCATTGAAGAAGAGAAAACTGGCGTCCAATCTAt                                                             >  2:694508/1‑141 (MQ=255)
                                                                                  cgcATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTcg                                                          <  2:1075320/139‑1 (MQ=255)
                                                                                  cTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTcg                                                          <  1:1089172/141‑1 (MQ=255)
                                                                                  cTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTcg                                                          <  1:984839/141‑1 (MQ=255)
                                                                                  cTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTcg                                                          <  2:658844/141‑1 (MQ=255)
                                                                                   ccATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTcgc                                                         <  1:166470/140‑1 (MQ=255)
                                                                                      ttAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTa                                                      <  2:1013379/141‑1 (MQ=255)
                                                                                                 aaTCGCAGCAGGTGAGGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGTTGCAGGCGCTACAGAGATAAGCATTGAAGTAGATGAATCTGGCGTCCAATCTATTCGCGTAGTCGATGATgg                                           >  2:945407/1‑141 (MQ=255)
                                                                                                      cagcagGTGAAGTAGCTGAAAGACCGAGTTCTATTATGGAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCg                                      >  1:984477/1‑141 (MQ=255)
                                                                                                       agcagGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTTGATAATGGAAGCgg                                     <  1:576384/141‑1 (MQ=255)
                                                                                                           ggggAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGGTGGAAAATGCTATAGCTGCAGGCGCTACAGAAATAAGTATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAAtt                                 <  2:74114/138‑1 (MQ=255)
                                                                                                           ggTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAAtt                                 >  2:901462/1‑141 (MQ=255)
                                                                                                             tGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGa                               <  1:184870/141‑1 (MQ=255)
                                                                                                             ggAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGa                               <  1:834149/140‑1 (MQ=255)
                                                                                                              taaGTAGGTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGaa                              <  1:982112/140‑1 (MQ=255)
                                                                                                              gAAGTAGTTGAAAGACCGAGTTCTGGTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGaa                              <  1:186901/141‑1 (MQ=255)
                                                                                                              gAAGTAAGTGAAAGACCCAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATTCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTACTCGCGTAGTCGATAATGGAAGCGGAATTGaa                              <  2:1054236/141‑1 (MQ=255)
                                                                                                                       gAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATGTATGCAGGCGTTACAGAAATAATCATTGAAGTAGAGGAATCTGGCGTCCAAGCTATTCGCGTATTCGATAATTGAAGCGGAATTGAAGCGGATGAc                     <  2:388373/141‑1 (MQ=255)
                                                                                                                       gAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAAGCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGAc                     <  2:347699/141‑1 (MQ=255)
                                                                                                                       gAAAGACCGAGCACTGTTGTGAAAGAACTGTTGGACATTGGTAGAGATGCTGTGGCTATCGAAATAAGCAGTGAAGTAGATGGAACTTTCGACCAATCTATTCGCGTAGTCGTTTATTGCAGCGGAATTGAAGCGGAAGAc                     <  2:388222/141‑1 (MQ=255)
                                                                                                                              cGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGAt              <  1:34717/141‑1 (MQ=255)
                                                                                                                              cGAGTTCTGTTGGGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGAt              <  1:34604/141‑1 (MQ=255)
                                                                                                                                aagtCTGTTGTGAAAGAACGAGTTGAAAATGCTAGAGGTGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTATTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTa            <  2:113983/138‑1 (MQ=255)
                                                                                                                                aGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTAGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGCCTTAGGATTa            >  2:494839/1‑141 (MQ=255)
                                                                                                                                 gTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAg           <  1:637970/141‑1 (MQ=255)
                                                                                                                                  ttCTGTAGTGAAAGAACTGTTGGAAAATGCTATAGATACAGGAGCTACAGAAATAAGCGGTGTGAGAGAGGAATCTTGTGTCCAATCTATTCGCGTAGTGGCAAATGGAAGCGGAATTGGAGCGGAAGACAGAGGATTGGt          >  2:396581/1‑141 (MQ=255)
                                                                                                                                   tCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGGTTAGTa         >  1:589470/1‑141 (MQ=255)
                                                                                                                                   tCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTa         >  1:663954/1‑141 (MQ=255)
                                                                                                                                         gtgGAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTATTTCAt   <  1:328086/141‑1 (MQ=255)
                                                                                                                                         gtgGAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTATTTCAt   <  1:129638/141‑1 (MQ=255)
                                                                                                                                          tgAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGTAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGAGTTAGGATTAGTATTTCATa  <  1:618603/141‑1 (MQ=255)
                                                                                                                                          ggAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTATTTCATa  <  1:313059/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AGGCATTGGTGAAGTTAAGTGAATTACAAAATCAATTAAAATAGAGGTGTTGCAAAATGGGGAAAATTAAAGAACTCCAAACCTCATTAGCAAATAAAATCGCAGCAGGTGAAGTAGTTGAAAGACCGAGTTCTGTTGTGAAAGAACTGTTGGAAAATGCTATAGATGCAGGCGCTACAGAAATAAGCATTGAAGTAGAGGAATCTGGCGTCCAATCTATTCGCGTAGTCGATAATGGAAGCGGAATTGAAGCGGAAGACTTAGGATTAGTATTTCATA  >  USA300TCH1516_ALE/1323173‑1323451

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