Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,231,407 G→A E18K (GAA→AAA)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,4070GA77.9% 269.3 / 79.8 137E18K (GAA→AAA) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base A (57/49);  total (67/70)
Fisher's exact test for biased strand distribution p-value = 6.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 138 total aligned reads displayed.

TACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAAT  >  USA300TCH1516_ALE/1231270‑1231547
                                                                                                                                         |                                                                                                                                            
tACGTAGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGGAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAg                                                                                                                                           >  4:22384/1‑141 (MQ=255)
   gtcgtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATg                                                                                                                                        >  1:25504/1‑141 (MQ=255)
     cgtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCg                                                                                                                                      >  1:519149/1‑141 (MQ=255)
            aaGGTTCGGATATTGAAGCAGCTTTTGGTAAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTa                                                                                                                               <  1:389437/141‑1 (MQ=255)
             aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAt                                                                                                                              >  2:395318/1‑141 (MQ=255)
               gTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAttt                                                                                                                            <  1:879322/141‑1 (MQ=255)
                  cGGATATTGACGCAGCTTGTGGTCCCATAAGAGCAAAGGAACGACACGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGACGGAAGATGCGGGTGGCAtcattt                                                                                                                         >  2:962745/1‑136 (MQ=255)
                  cGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGTAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAtttttt                                                                                                                         >  2:689412/1‑141 (MQ=255)
                  cGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTAtttttt                                                                                                                         <  2:429261/141‑1 (MQ=255)
                  cGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTAtttttt                                                                                                                         >  1:985014/1‑141 (MQ=255)
                  cGGATATTGACGCAGATTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTAtttttt                                                                                                                         <  2:429278/141‑1 (MQ=255)
                   ggATATTGACGCATCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGCATAAAGATGCGGGTGGTATTTTTta                                                                                                                        <  2:466643/141‑1 (MQ=255)
                    gATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTtat                                                                                                                       >  1:242266/1‑141 (MQ=255)
                     atatTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTtata                                                                                                                      >  2:812595/1‑141 (MQ=255)
                       atTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAAt                                                                                                                    >  1:723493/1‑141 (MQ=255)
                       atTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAAt                                                                                                                    <  1:375303/141‑1 (MQ=255)
                        ttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATc                                                                                                                   >  4:200914/1‑141 (MQ=255)
                        ttGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATc                                                                                                                   >  4:343675/1‑141 (MQ=255)
                         tGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCa                                                                                                                  >  3:151087/1‑141 (MQ=255)
                            cGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAAc                                                                                                               >  3:162755/1‑141 (MQ=255)
                              cAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTa                                                                                                             <  3:235946/141‑1 (MQ=255)
                                 cTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATc                                                                                                          <  1:356486/141‑1 (MQ=255)
                                      ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                     >  3:309750/1‑141 (MQ=255)
                                      ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATcaacaa                                                                                                     <  3:226496/141‑1 (MQ=255)
                                         cAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACtt                                                                                                  >  1:473187/1‑141 (MQ=255)
                                          aaTTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACttt                                                                                                 <  1:660078/141‑1 (MQ=255)
                                            ttAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTa                                                                                               <  1:401869/141‑1 (MQ=255)
                                              aaGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGt                                                                                             <  4:445505/141‑1 (MQ=255)
                                                gagCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTc                                                                                           >  2:19149/1‑141 (MQ=255)
                                                  gCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCtg                                                                                         >  1:443510/1‑141 (MQ=255)
                                                     aaGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACttt                                                                                                 <  3:171495/130‑1 (MQ=255)
                                                       ggAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGa                                                                                    <  3:266757/141‑1 (MQ=255)
                                                        gAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGAt                                                                                   <  2:357311/141‑1 (MQ=255)
                                                         aaCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATg                                                                                  >  2:168133/1‑141 (MQ=255)
                                                          aCGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATgg                                                                                 <  3:408859/141‑1 (MQ=255)
                                                           cGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGt                                                                                >  1:118439/1‑141 (MQ=255)
                                                               aaGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgg                                                                            <  4:389769/141‑1 (MQ=255)
                                                               aaGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgg                                                                            >  1:539770/1‑141 (MQ=255)
                                                                aGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATggg                                                                           <  4:223878/141‑1 (MQ=255)
                                                                  tAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACAAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTg                                                                         <  1:758595/141‑1 (MQ=255)
                                                                   aGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTgg                                                                        <  2:967140/141‑1 (MQ=255)
                                                                    gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                       <  2:743956/141‑1 (MQ=255)
                                                                    gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                       >  2:706450/1‑141 (MQ=255)
                                                                    gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                       <  2:760514/141‑1 (MQ=255)
                                                                     aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                      <  1:600226/141‑1 (MQ=255)
                                                                     aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                      <  3:27868/141‑1 (MQ=255)
                                                                     aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                      >  3:128451/1‑141 (MQ=255)
                                                                     aaaCGAGGTAAAGAAAAATGATAGAGGAAAAATTTTTTACTGATAATGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTACTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                      <  4:36669/141‑1 (MQ=255)
                                                                      aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCa                                                                     <  4:323488/141‑1 (MQ=255)
                                                                       aagAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCAt                                                                    <  1:763997/139‑1 (MQ=255)
                                                                       aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGAGAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCAt                                                                    <  1:1013537/141‑1 (MQ=255)
                                                                         gAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATaa                                                                  <  2:78264/141‑1 (MQ=255)
                                                                          aGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATaaa                                                                 >  2:256840/1‑141 (MQ=255)
                                                                              aaaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAg                                                             >  4:337119/1‑141 (MQ=255)
                                                                               aaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATATAGATAAGAATGACGCTGCGGGTGGTATTTTTTATAATCACACTACTCAACAAATTTTAGTCCTGTGTTATGGTATGGGTGGCCATAAAGCAgg                                                            <  1:153088/141‑1 (MQ=255)
                                                                                aGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGga                                                           >  3:157600/1‑141 (MQ=255)
                                                                                 gACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                          <  2:816300/141‑1 (MQ=255)
                                                                                 gACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                          <  1:1046668/141‑1 (MQ=255)
                                                                                  aCAAATGCTAGAGGCACAATTTTTTACTTATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgaga                                                         >  2:598072/1‑141 (MQ=255)
                                                                                  aCAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgaga                                                         >  1:4974/1‑141 (MQ=255)
                                                                                      aTGCTAGAGGCACAATTTTTTACTTTTACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGtt                                                     >  2:854832/1‑141 (MQ=255)
                                                                                       tGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTg                                                    >  4:27486/1‑141 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACGAGATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                  <  4:169573/141‑1 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                  >  1:420798/1‑141 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                  >  2:282976/1‑141 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                  >  2:283098/1‑141 (MQ=255)
                                                                                         cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                  <  1:149447/141‑1 (MQ=255)
                                                                                           agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAg                                                <  4:326188/141‑1 (MQ=255)
                                                                                            gagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGt                                               <  2:990321/141‑1 (MQ=255)
                                                                                             agGCACAATTTTTTACTGATACTGGACAGCATAGAGATAAGAATAAAGAGGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                                              <  2:90698/141‑1 (MQ=255)
                                                                                             agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                                              >  1:960696/1‑141 (MQ=255)
                                                                                                   aaTTTTTTACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTg                                        >  1:587622/1‑141 (MQ=255)
                                                                                                   aaTTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCCTAACGCAGGAGAAGTTGCAAGTAAACTTg                                        >  1:35057/1‑141 (MQ=255)
                                                                                                   aaTTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTg                                        >  2:804526/1‑141 (MQ=255)
                                                                                                    aTTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGt                                       <  1:146835/141‑1 (MQ=255)
                                                                                                     ttttttACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGtt                                      <  1:945111/141‑1 (MQ=255)
                                                                                                      tttttACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                     >  3:136501/1‑141 (MQ=255)
                                                                                                      tttttACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTa                                     >  3:195748/1‑141 (MQ=255)
                                                                                                        tttACTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACa                                   <  4:405666/141‑1 (MQ=255)
                                                                                                           aCTGATACTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGAt                                <  3:344736/141‑1 (MQ=255)
                                                                                                            cTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATg                               >  1:272795/1‑141 (MQ=255)
                                                                                                                  ctagACAACATAGAGATAAGAAAAAAGATGCGGGTGGTATTTTTTATAATAAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                         <  4:434636/138‑1 (MQ=255)
                                                                                                                  cTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                         <  1:68018/141‑1 (MQ=255)
                                                                                                                  cTGGACAACATATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                         <  1:650337/141‑1 (MQ=255)
                                                                                                                   tGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGaa                        <  1:305088/141‑1 (MQ=255)
                                                                                                                       caacaTATAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTCTTAGTTCTGTGTGATGGTATGGGTGGACATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATCAGTTGAAATcc                    >  2:790817/1‑141 (MQ=255)
                                                                                                                         acaTAGAGATAAGAATAAAGATGCGGGTGGTCTTTTTTATACTCAAACTAATCAACAACTTTAAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGCAGTGGCAAGTAAAATTGTTACAGATGATTTGAAATCCCg                  >  2:618314/1‑141 (MQ=255)
                                                                                                                         acaTAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                  <  1:442696/141‑1 (MQ=255)
                                                                                                                         acaTAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                  <  2:921274/141‑1 (MQ=255)
                                                                                                                         acaTAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                  <  2:723493/141‑1 (MQ=255)
                                                                                                                            tAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt               <  2:519149/141‑1 (MQ=255)
                                                                                                                            tAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt               <  2:538110/141‑1 (MQ=255)
                                                                                                                            tAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt               <  3:23665/141‑1 (MQ=255)
                                                                                                                              gagaTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTAGTTACAGATGAGTTGAAATCCCGtttgg             >  4:141512/1‑139 (MQ=255)
                                                                                                                                 aTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg          <  1:706450/141‑1 (MQ=255)
                                                                                                                                   aaGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg        >  2:435732/1‑141 (MQ=255)
                                                                                                                                   aaGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg        >  2:460456/1‑141 (MQ=255)
                                                                                                                                   aaGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg        >  2:34457/1‑141 (MQ=255)
                                                                                                                                    aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAAGTAATCAACAACTTTTAGTTCACTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATCTGTTAGCGATGAGTTAACATCCCGTTTTGAAGCgg       >  2:570220/1‑141 (MQ=255)
                                                                                                                                    aGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCgg       >  3:89679/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
TACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAAT  >  USA300TCH1516_ALE/1231270‑1231547

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