Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,752,368 G→A 5.1% G276G (GGC→GGT ruvB ← Holliday junction ATP‑dependent DNA helicase RuvB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,752,3680GA 5.1% 399.5 / 5.1 127G276G (GGC→GGTruvBHolliday junction ATP‑dependent DNA helicase RuvB
Reads supporting (aligned to +/- strand):  ref base G (61/59);  new base A (5/2);  total (66/61)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 127 total aligned reads displayed.

TAATGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTGGAAG  >  USA300TCH1516_ALE/1752233‑1752506
                                                                                                                                       |                                                                                                                                          
tAATGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCAtt                                                                                                                                       >  3:296786/1‑141 (MQ=255)
 aaTGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATta                                                                                                                                      <  2:633572/141‑1 (MQ=255)
 aaTGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATta                                                                                                                                      <  2:1029035/141‑1 (MQ=255)
  aTGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATtat                                                                                                                                     <  2:723765/141‑1 (MQ=255)
  aTGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATtat                                                                                                                                     <  1:131952/141‑1 (MQ=255)
   tGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATtata                                                                                                                                    <  2:773581/141‑1 (MQ=255)
    ggTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATAc                                                                                                                                   >  2:1044974/1‑141 (MQ=255)
      tgtTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTg                                                                                                                                 >  2:456638/1‑141 (MQ=255)
      tgtTGCTATTCTGCCACGTGGCGTACGTTCTAAAAAGCCTATCTGAATAAGAAATGGCACATAAACGTCCTCAATTGTAATACGTTCTTCACCAATAGTTACAGCAATCGTAACTAAACCAACAAGTCCGCCATTATAcag                                                                                                                                 >  2:356530/1‑139 (MQ=255)
        ttGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGGATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCt                                                                                                                               <  1:852741/141‑1 (MQ=255)
         tGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCtt                                                                                                                              >  2:907174/1‑141 (MQ=255)
         tGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATAGTATCTAAACCAACAGGTCCGCCATTATACTGCtt                                                                                                                              <  1:777233/141‑1 (MQ=255)
                gtGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAataata                                                                                                                       >  2:882716/2‑141 (MQ=255)
                cTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAataata                                                                                                                       <  1:494788/141‑1 (MQ=255)
                  gCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACa                                                                                                                     <  1:519019/141‑1 (MQ=255)
                   ccACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAg                                                                                                                    >  4:266817/1‑141 (MQ=255)
                   ccACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAg                                                                                                                    <  3:389813/141‑1 (MQ=255)
                   ccACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAg                                                                                                                    <  2:730793/141‑1 (MQ=255)
                     aCGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGtt                                                                                                                  <  3:132486/141‑1 (MQ=255)
                         ggCGTAAGTTCTAAAAAGTCTTTCTGAATAAGAAATGGCAAATAAACGTACCAAATTGTAATACGTTATTCACCAATTGTTACGGCAATCGTATTTAAAACAACAGGTCCAACATTATACTGCTTAATAATACAGTtcatc                                                                                                              <  4:5009/141‑1 (MQ=255)
                          gCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTtcatca                                                                                                             >  2:548896/1‑141 (MQ=255)
                          gCGTACGGTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGGCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTtcatca                                                                                                             >  2:750940/1‑141 (MQ=255)
                           cGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTcatcat                                                                                                            >  3:242455/1‑141 (MQ=255)
                              aCGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCAtttt                                                                                                         <  2:937931/141‑1 (MQ=255)
                               cGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTa                                                                                                        <  3:275312/141‑1 (MQ=255)
                                  tCTAAAAAGCCTTACTGAATAAGAAATGGCTCATAAACGTCCTCATTTGTAATACTTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGa                                                                                                     <  3:179259/141‑1 (MQ=255)
                                    tAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATc                                                                                                   <  2:114399/141‑1 (MQ=255)
                                     aaaaaGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCa                                                                                                  >  2:675345/1‑141 (MQ=255)
                                      aaaaGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCaa                                                                                                 <  1:218844/141‑1 (MQ=255)
                                          gCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTa                                                                                             >  1:632765/1‑141 (MQ=255)
                                              ttCTGAATAAGAAATGGCTCATCAACGTCCTCAAATGTAATACGTTCTCCACCAATTGTTACGGAAATCGTATATAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCGTCATTTTATGATCAATGTCATCt                                                                                         <  4:50554/141‑1 (MQ=255)
                                              ttCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCTGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCt                                                                                         <  2:1005336/141‑1 (MQ=255)
                                              ttCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCACTGGAATCt                                                                                         >  2:314551/1‑141 (MQ=255)
                                              ttCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCt                                                                                         <  1:405494/141‑1 (MQ=255)
                                              ttCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCt                                                                                         >  2:286504/1‑141 (MQ=255)
                                               tCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTa                                                                                        >  1:275002/1‑141 (MQ=255)
                                               tCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTa                                                                                        >  3:171931/1‑141 (MQ=255)
                                               tCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTa                                                                                        >  2:1030841/1‑141 (MQ=255)
                                                    aTAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTCATCTAGTCCg                                                                                   >  3:139306/1‑141 (MQ=255)
                                                      aaGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCgtg                                                                                 >  2:754514/1‑141 (MQ=255)
                                                      aaGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCgtg                                                                                 <  4:29114/141‑1 (MQ=255)
                                                      aaGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCgtg                                                                                 >  1:824897/1‑141 (MQ=255)
                                                       aGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTTACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTACgtgt                                                                                >  4:59783/1‑141 (MQ=255)
                                                        gAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGtt                                                                               <  1:626503/141‑1 (MQ=255)
                                                          aaTGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGa                                                                             <  1:692383/141‑1 (MQ=255)
                                                           aTGGCTCATACACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGAt                                                                            <  2:445015/141‑1 (MQ=255)
                                                           aTGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGAt                                                                            <  1:240420/141‑1 (MQ=255)
                                                           aTGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGAt                                                                            <  1:750940/141‑1 (MQ=255)
                                                           aTGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGAt                                                                            <  2:9948/141‑1 (MQ=255)
                                                                 cATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACtt                                                                      <  4:144112/141‑1 (MQ=255)
                                                                  aTAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTg                                                                     <  2:289376/141‑1 (MQ=255)
                                                                  aTAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTg                                                                     >  4:256126/1‑141 (MQ=255)
                                                                    aaaCGTCCTCAAGTGTAATACGTTCTTCACCAATTGTCACGGCCATCGTATCTGAACGAACAGGACCACCATTATACAGCCTAATAATACAGTTCAGCATTTTATGATCAATGTAATCTAGTCCGTGTTGCTTAACTTGaa                                                                   >  4:220646/1‑141 (MQ=255)
                                                                      aCGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGt                                                                 <  1:1029829/141‑1 (MQ=255)
                                                                       cGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATATTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTa                                                                >  1:852541/1‑141 (MQ=255)
                                                                        gTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTaa                                                               >  4:3295/1‑141 (MQ=255)
                                                                          ccTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAAc                                                             >  1:1072372/1‑141 (MQ=255)
                                                                               aTTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAAt                                                        <  1:67623/141‑1 (MQ=255)
                                                                               aTTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAAt                                                        >  1:262174/1‑141 (MQ=255)
                                                                                     aTACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGc                                                  >  2:366804/1‑141 (MQ=255)
                                                                                     aTACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGc                                                  >  2:366798/1‑141 (MQ=255)
                                                                                     aTACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGc                                                  >  2:299897/1‑141 (MQ=255)
                                                                                      tACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCt                                                 >  2:686261/1‑141 (MQ=255)
                                                                                       aCGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCtt                                                >  1:521857/1‑141 (MQ=255)
                                                                                          ttcttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgt                                             <  4:152463/141‑1 (MQ=255)
                                                                                           tcttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgtt                                            >  3:425042/1‑141 (MQ=255)
                                                                                            cttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgttg                                           >  3:312633/1‑141 (MQ=255)
                                                                                            cttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgttg                                           >  2:485408/1‑141 (MQ=255)
                                                                                            cttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgttg                                           >  3:350697/1‑141 (MQ=255)
                                                                                            cttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTCATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTGATGAGTGCTTCgttg                                           >  3:155455/1‑141 (MQ=255)
                                                                                             ttcACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCgttgt                                          >  1:571300/1‑141 (MQ=255)
                                                                                               cACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGttt                                        >  4:326129/1‑141 (MQ=255)
                                                                                                 ccAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCa                                      >  3:125557/1‑141 (MQ=255)
                                                                                                  cAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCaa                                     >  4:134068/1‑141 (MQ=255)
                                                                                                   aaTTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAAt                                    <  1:383056/141‑1 (MQ=255)
                                                                                                   aaTTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAAt                                    <  1:1079192/141‑1 (MQ=255)
                                                                                                    aTTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATg                                   <  1:548896/141‑1 (MQ=255)
                                                                                                       gTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCGAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGGTTCGTTGTTTCAATGtat                                >  4:343000/1‑141 (MQ=255)
                                                                                                       gTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATCATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCTACTTGAAGTAAACATAATGCGTGGTTCGCTGTTTCAATGtat                                >  1:1059422/1‑141 (MQ=255)
                                                                                                         tACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTatat                              <  1:828283/141‑1 (MQ=255)
                                                                                                          aCGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATAtt                             >  1:113440/1‑141 (MQ=255)
                                                                                                            ggCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTg                           <  1:903809/141‑1 (MQ=255)
                                                                                                             gCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGt                          <  1:243690/141‑1 (MQ=255)
                                                                                                              cAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGtt                         <  2:1063996/141‑1 (MQ=255)
                                                                                                              cAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGtt                         <  1:849441/141‑1 (MQ=255)
                                                                                                              cAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGtt                         >  2:455810/1‑141 (MQ=255)
                                                                                                                  cGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATc                     <  1:437114/141‑1 (MQ=255)
                                                                                                                     aTCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttc                  >  1:742228/1‑141 (MQ=255)
                                                                                                                     aTCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttc                  >  2:857982/1‑141 (MQ=255)
                                                                                                                     aTCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttc                  >  3:20054/1‑141 (MQ=255)
                                                                                                                      tCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttca                 <  4:14357/141‑1 (MQ=255)
                                                                                                                       cTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttcat                <  1:378052/141‑1 (MQ=255)
                                                                                                                         aaaCCAACAGGTCCGCCAGTAGACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAGACCTAATGCGTGCTTCGTTGTGTCAATGTATATTTGTTCATCttcattc              <  2:373309/141‑1 (MQ=255)
                                                                                                                          aaCCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCttcattca             >  4:134727/1‑141 (MQ=255)
                                                                                                                               aCAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTTTTCATCTTCATTCACTTGc        >  1:287843/1‑141 (MQ=255)
                                                                                                                                 aGGTCCACCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTg      >  2:636368/1‑141 (MQ=255)
                                                                                                                                  ggTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTgg     >  3:448663/1‑141 (MQ=255)
                                                                                                                                  ggTCCACCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTgg     >  2:123356/1‑141 (MQ=255)
                                                                                                                                    tCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTGGaa   <  4:314910/141‑1 (MQ=255)
                                                                                                                                     ccGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTGGAAg  >  3:315168/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                          
TAATGGTGTTGCTTTTCTGCCACGTGGCGTACGTTCTAAAAAGCCTTTCTGAATAAGAAATGGCTCATAAACGTCCTCAATTGTAATACGTTCTTCACCAATTGTTACGGCAATCGTATCTAAACCAACAGGTCCGCCATTATACTGCTTAATAATACAGTTCATCATTTTATGATCAATGTAATCTAGTCCGTGTTGATCAACTTGAAGTAAACCTAATGCGTGCTTCGTTGTTTCAATGTATATTTGTTCATCTTCATTCACTTGCTGGAAG  >  USA300TCH1516_ALE/1752233‑1752506

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