Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,407 G→A 6.5% 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,4070GA 6.5% 228.1 / 4.7 77E18K (GAA→AAA) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base A (2/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAA  >  USA300TCH1516_ALE/1231267‑1231545
                                                                                                                                            |                                                                                                                                          
aaaTACGTCGTGAAAAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATg                                                                                                                                            <  2:258589/140‑1 (MQ=255)
   tACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAAAGATAAGAATgcaa                                                                                                                                         >  2:772909/1‑138 (MQ=255)
    aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATATTGGACAACATAGCGATAAGAATGAAGa                                                                                                                                        >  1:401254/1‑141 (MQ=255)
    aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGa                                                                                                                                        >  1:344588/1‑141 (MQ=255)
    aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGa                                                                                                                                        >  1:2955/1‑141 (MQ=255)
    aCGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGATAGGGGCTCAATTTTTTATTGATACTGCAAGACATAAAGCTAAGAATGAAGa                                                                                                                                        >  1:830361/1‑141 (MQ=255)
        cgtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCg                                                                                                                                    <  1:587657/141‑1 (MQ=255)
         gtGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCgg                                                                                                                                   >  1:335669/1‑141 (MQ=255)
           gAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCgggg                                                                                                                                 >  2:446042/1‑140 (MQ=255)
            aaCAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGggtg                                                                                                                                <  2:301609/141‑1 (MQ=255)
               aaGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTa                                                                                                                             <  2:219413/141‑1 (MQ=255)
                aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGCGGCACAATTTTTTACTGATAGTGGACGAAATAGAGGTAAGAATGAAGATGCGGGTGGTAt                                                                                                                            >  1:716990/1‑141 (MQ=255)
                aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTAt                                                                                                                            >  1:821526/1‑141 (MQ=255)
                aGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTAGACAACATAGAGATAAGAAAGAAGATGGGGGTGGAAt                                                                                                                            >  2:541044/1‑141 (MQ=255)
                   ttCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTAtttt                                                                                                                         <  1:176275/141‑1 (MQ=255)
                     cGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTAtttttt                                                                                                                       <  2:404657/141‑1 (MQ=255)
                          atTGAGGCAGCTTGAGGTCAATGAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTAGTATTTTTTATAAt                                                                                                                  <  2:640613/141‑1 (MQ=255)
                                    cTTGGGGTCAATTAAGATCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATATTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATc                                                                                                        <  2:632275/141‑1 (MQ=255)
                                         ggTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTGTTTATAATCAAACTAATcaacca                                                                                                   >  1:118181/1‑139 (MQ=255)
                                            cAATTAAGAGCAAAGGAACGACAAGTAGGAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACtt                                                                                                <  2:658791/141‑1 (MQ=255)
                                                         aGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCtgtgtg                                                                                   <  1:863496/141‑1 (MQ=255)
                                                             aCGACAAGTAGATACGAGGTAAAGACAAATGCTAGAGGCACAAGTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATACTCAAACTAATCAACAGCCTTCAGTTCTGCGTGAtag                                                                               <  2:790344/141‑3 (MQ=255)
                                                             aCGACAAGTAGAAACGAGGTAAAGACAAATACTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAACCACCCTAATCATCAACCTTTAGTCCTGTGTGATgg                                                                               <  2:85848/141‑1 (MQ=255)
                                                               aaCAAGAAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTa                                                                             <  2:232004/140‑1 (MQ=255)
                                                                aggaGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTAt                                                                            <  2:231647/138‑1 (MQ=255)
                                                                aaaaGTAGAAACCGGGTAAAGACAAATGCTAGAGGCACAATTTTTTACGGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTTGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                            <  2:31981/139‑1 (MQ=255)
                                                                aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                            >  1:208187/1‑141 (MQ=255)
                                                                aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATAAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                            >  2:116114/1‑141 (MQ=255)
                                                                aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGAAACTGTACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                            <  1:26727/141‑1 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATg                                                                           <  2:870968/141‑1 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                           >  1:392580/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                           >  2:282117/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTATg                                                                           >  1:766846/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAAATTTTAGTTCTGTGTGATGGTATg                                                                           >  2:425447/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAAAAGAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATg                                                                           >  2:287129/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATGTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGAAGCGGGTGGTATTTTTTATAATCAAACTAATAANCAANTNTTAGTCTTGTGTGATGGTATg                                                                           >  1:486228/1‑141 (MQ=255)
                                                                 cAAGTAGAAACGAGCTAAAGACAAATGCTAGCGGCCCAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATg                                                                           <  2:867881/141‑1 (MQ=255)
                                                                    gTAGAAACGGGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                        <  2:143985/141‑1 (MQ=255)
                                                                    gTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                        <  2:567734/141‑1 (MQ=255)
                                                                      aGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTgg                                                                      <  1:781462/141‑1 (MQ=255)
                                                                       gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                     >  1:372698/1‑141 (MQ=255)
                                                                       gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGc                                                                     >  1:552862/1‑141 (MQ=255)
                                                                       gAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTGTTAGTTCTGTGTGATGGTATGGGTGGc                                                                     >  1:403299/1‑141 (MQ=255)
                                                                              ggTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGATATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAg                                                              >  1:202900/1‑141 (MQ=255)
                                                                              ggTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGAATGGGTGGCCATaaac                                                              >  2:284948/1‑140 (MQ=255)
                                                                               gTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGc                                                             <  1:322265/141‑1 (MQ=255)
                                                                                 aaaGACAAATGCTATAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTAGCCATAAAGCAg                                                           <  1:606125/141‑1 (MQ=255)
                                                                                 aaaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAg                                                           >  1:43822/1‑141 (MQ=255)
                                                                                    gACAAATGCTAGCGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                        <  2:479314/141‑1 (MQ=255)
                                                                                    gACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGgag                                                        <  2:641979/141‑1 (MQ=255)
                                                                                            cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                <  2:678433/141‑1 (MQ=255)
                                                                                            cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                >  2:470547/1‑141 (MQ=255)
                                                                                            cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAGGGGTGGCCATAAATCAGGAGAAGTTGcg                                                >  1:195954/1‑140 (MQ=255)
                                                                                              agagGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAg                                              <  2:43822/141‑1 (MQ=255)
                                                                                                   cacaATTTTTTACTGATACTGGACAACATAGAGATCAGACTGACGAAAAGGGCGGTATTCTTTCTACTCAAAATAATAGAAGATAATCACGTTTGCGCAGTGGCACGGGTGCCCATAAAGCATGATAAGGTGCAAGTAAAt                                         >  2:315396/1‑141 (MQ=255)
                                                                                                   cacaATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAt                                         >  1:79735/1‑141 (MQ=255)
                                                                                                     caATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAAttt                                       >  1:177733/1‑141 (MQ=255)
                                                                                                          ttttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCGTAAAGCAGGAGAAGTTGCAAGTAAATTTTTTAc                                  >  1:436399/1‑141 (MQ=255)
                                                                                                            ttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAAGCAAACTAATCAACAACTTGTAGTTCTGTGTGATGGTAGGGGTGGCCATAAAGCGGGAGAAGTTGCCAGTAAATTTGTTACAg                                >  1:902276/1‑141 (MQ=255)
                                                                                                              aCTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGATATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGAt                              <  2:675182/141‑1 (MQ=255)
                                                                                                               cTGATACTGGACAACATAGAGATAAGCATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTACATAAACGTAAAGGTATGGGAGGCCAGAAAGCAGGAGAAGTTACAAGTAAATTTGTTACAGATg                             >  1:926276/1‑141 (MQ=255)
                                                                                                                  aTACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAAATAATCAACAACTGTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGt                          <  2:335669/141‑1 (MQ=255)
                                                                                                                   tgcTGTACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTATAATCAAACTAATGAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTg                        <  2:231848/139‑1 (MQ=255)
                                                                                                                         acaacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCCGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGACATc                   >  1:206823/1‑141 (MQ=255)
                                                                                                                         acaacaTAGAGATAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATc                   >  1:41336/1‑141 (MQ=255)
                                                                                                                           aacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGTGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATccc                 >  1:264043/1‑141 (MQ=255)
                                                                                                                              aTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTTATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGtt              >  1:151463/1‑141 (MQ=255)
                                                                                                                              aTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTATTAGTTCTGTGTTATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTCAATTTGTTACAGATGAGTTGAAATCCCgct              >  1:724845/1‑139 (MQ=255)
                                                                                                                                 gagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg           >  1:74564/1‑141 (MQ=255)
                                                                                                                                 gagaTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAATAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTg           >  1:300826/1‑141 (MQ=255)
                                                                                                                                 gagaTAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGCGAGGTATTAAAGCAGGATAAGTTGCAAGAAAACTTGTTACAGATGCGTTGAAATCCCGTTTTg           <  2:429967/141‑1 (MQ=255)
                                                                                                                                    agaaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAACCACCAACTTTTAGTTGTGTGTGATGGTATGGGTGGCCAGAAAGCAGGAGAAGTTGCAAGTAAATTTGGTACAGAGGAGTTGAGATCCCTTGTTTAAg        <  2:309154/139‑1 (MQ=255)
                                                                                                                                    aTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg        <  2:448847/141‑1 (MQ=255)
                                                                                                                                    aTAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAg        <  1:882744/141‑1 (MQ=255)
                                                                                                                                     tAAGAATAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGc       <  1:331907/141‑1 (MQ=255)
                                                                                                                                      aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg      <  1:776831/141‑1 (MQ=255)
                                                                                                                                      aaGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCg      <  2:177733/141‑1 (MQ=255)
                                                                                                                                         aaTAAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGACGCGGaa   >  2:124701/1‑141 (MQ=255)
                                                                                                                                          aTGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACGTTTAGTTACGTGTGATGGTATGGGTGACCATACAGCAGGAGAAGTTGCAAGTAAAGTTATTACAGATGAGTGGAAATCGCGTATTGAGGCGGaaa  >  2:411099/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CAATACGTCGTGAACAAGGTTCGGATATTGACGCAGCTTGTGGTCAATTAAGAGCAAAGGAACGACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAA  >  USA300TCH1516_ALE/1231267‑1231545

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: