Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,231,470 A→C 6.4% M39L (ATG→CTG)  stp → Serine/threonine phosphatase stp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,231,4700AC 6.4% 227.1 / 5.0 86M39L (ATG→CTG) stpSerine/threonine phosphatase stp
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base C (3/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTCAG  >  USA300TCH1516_ALE/1231330‑1231610
                                                                                                                                            |                                                                                                                                            
gACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTCTTAGTTCTGTGAGATGGTa                                                                                                                                              >  1:406182/1‑141 (MQ=255)
aaCAAGTACAACAGAGGTAAACAAAAAAGCTAGAGGCAGAATTCATAAAAGATAATGAACAACATAGAGAACAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTa                                                                                                                                              <  1:839630/140‑1 (MQ=255)
 aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                                             <  2:530715/141‑1 (MQ=255)
 aCAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTAt                                                                                                                                             <  1:1175103/141‑1 (MQ=255)
   aaGTAGAAACGAGGAAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATgg                                                                                                                                           >  1:622542/1‑141 (MQ=255)
     aaaaaaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGt                                                                                                                                         <  1:647453/137‑1 (MQ=255)
      tAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACGTTTAGTTCTGTGTGATGGTATGGGTg                                                                                                                                        >  1:857282/1‑141 (MQ=255)
         aaaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGcc                                                                                                                                     >  2:1193478/1‑141 (MQ=255)
          aaCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCa                                                                                                                                    <  1:190037/141‑1 (MQ=255)
           aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAGCTTTTAGTTCTGTGTGATGGTCTGGGTGGCCAt                                                                                                                                   <  2:385886/141‑1 (MQ=255)
           aCGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCAt                                                                                                                                   <  1:1168664/141‑1 (MQ=255)
            tgAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACATCATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGGCATa                                                                                                                                  <  2:726259/140‑1 (MQ=255)
             gAGGTAAAGACAAATGCTAGATGCAGAGTTTTTTACTGATACTGGACAACATAGATATAATAATGAACATGCCGGTGGTATTTTTTATAATCAAACTAATCAACAACATTTAGTTCTGTGTGATGGTATGGGTGGCCATaa                                                                                                                                 <  2:1001150/141‑1 (MQ=255)
             gAGGGAAAGGCAATGGCTAGAGGCCCAATTTTTTCCTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATaa                                                                                                                                 <  2:270936/141‑1 (MQ=255)
              cagaaaaGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATaaa                                                                                                                                <  1:265849/137‑1 (MQ=255)
                ggaaaGACAAATTCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGTGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGc                                                                                                                              <  2:480284/139‑1 (MQ=255)
                 tAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCGTAAAGCa                                                                                                                             >  2:351347/1‑141 (MQ=255)
                    aGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGCATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATAAACCACTTTTAGTTATGTGTGATGGTAAGGGGGACCATAGAGCAggg                                                                                                                          >  2:630842/1‑140 (MQ=255)
                         aaTGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGt                                                                                                                     >  1:796053/1‑141 (MQ=255)
                          aTGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGtt                                                                                                                    >  1:479219/1‑141 (MQ=255)
                          aTGCTAGAGGCACAATTTTTTACTGATACTGGACAAAATAGAGATAGGAATGAAGATGCGGGTGGTATTTTTTATAATCAAAGTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGTGAAGtt                                                                                                                    >  2:787128/1‑141 (MQ=255)
                             cTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCa                                                                                                                 <  2:568651/141‑1 (MQ=255)
                                 agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGGTGGTATGGGTGGCCATCAAGCAGGAGAAGTTGCAAGTa                                                                                                             >  1:135126/1‑141 (MQ=255)
                                 agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTCCGGGTGCCGTAAAAGCAGGAGAAGTTGCAAGTa                                                                                                             >  2:537060/1‑141 (MQ=255)
                                 agGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTa                                                                                                             <  2:452442/141‑1 (MQ=255)
                                   gCACGAATTTTTAATGATACTGGACAAAATAGTGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTaaa                                                                                                           <  2:82621/141‑1 (MQ=255)
                                   gCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTaaa                                                                                                           >  1:1079192/1‑141 (MQ=255)
                                      caATTTTTTACTGATACTGGTCAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATAAACAACTTTTAGTTCTGTGTAATGGTATGGGTGGCCTTAAAGCAGGAGAGGTTGTGAGTAAAttg                                                                                                        >  2:942539/1‑140 (MQ=255)
                                         ttttttACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGtt                                                                                                     >  1:1080096/1‑141 (MQ=255)
                                          tttttACTGATACTGGACAACATAGAGATAAGAATGAAAATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTGTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTACATTTGTTa                                                                                                    >  2:803928/1‑141 (MQ=255)
                                              tACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGa                                                                                                >  1:946533/1‑141 (MQ=255)
                                                cTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAGGCAGGAGAAGTTGCAAGTAAATTTGTTACAGagg                                                                                              >  2:374795/1‑139 (MQ=255)
                                                cTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATg                                                                                              >  1:505935/1‑141 (MQ=255)
                                                 tGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAAGCAACAACTTTTAGTTCTGTGTGATGTTAGTGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGa                                                                                             <  1:169731/141‑1 (MQ=255)
                                                  gATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTGGTTCTGTTTGATGGTCTGGGTTGACATAAAGAAGGGGCATTTGCAAGTTAAGTTGTTACAGATGAg                                                                                            >  2:748019/1‑141 (MQ=255)
                                                  gATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTGCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAg                                                                                            <  2:505585/141‑1 (MQ=255)
                                                     aCTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTg                                                                                         <  1:926923/141‑1 (MQ=255)
                                                      cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                                                                        >  1:160853/1‑141 (MQ=255)
                                                      cTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGa                                                                                        >  2:470099/1‑141 (MQ=255)
                                                            aacaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTGTTAGTTCTGTGTGAAGGTATGGGTGGCCTTAAAGCAGGGGAAGTTGCAAGTAAATTTGTTACAGATGAGTGGAGAtacg                                                                                  >  1:390753/1‑138 (MQ=255)
                                                             acaTAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCg                                                                                 <  2:713974/141‑1 (MQ=255)
                                                                tAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGttt                                                                              >  1:479820/1‑141 (MQ=255)
                                                                    gctAAGAATGAAGATGCGGGTGGTATTTTTTAGAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGaa                                                                          <  1:327809/139‑1 (MQ=255)
                                                                        aGAATGAAGAGGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACATTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATTAATTGAAATCCCGTTTTGAAGCgg                                                                      >  2:525842/1‑141 (MQ=255)
                                                                          aaTGAAGATGGGGGTGGTATTTTTTATAATACAACTATGAAACAACGTTTAGTTCTGTGTGATGGTATGGGTGCCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATTCCGTTTTGAAGCGGaa                                                                    <  2:538460/141‑1 (MQ=255)
                                                                            tGAAGATGCGGGGGGTATTTTTTATAATCAAACTAATCAACGACTTTTAGTTATGTGTGATGGTATGGGGGGCCCTCCAGCAGGAGACGTTGCTGGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGaaaa                                                                  >  2:348184/1‑141 (MQ=255)
                                                                             gAAGATGCGGGTGGTATTTTTTATAATCAAACTAAGCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAAt                                                                 <  2:861782/141‑1 (MQ=255)
                                                                                   gCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtata                                                           >  2:75873/1‑141 (MQ=255)
                                                                                   gCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTtata                                                           >  2:952847/1‑141 (MQ=255)
                                                                                    cGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAg                                                          <  1:482306/141‑1 (MQ=255)
                                                                                      ggtggtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGGCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAAATCTATAGaa                                                        >  2:1000783/1‑141 (MQ=255)
                                                                                       gtggtATTTTTTATAAGCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaac                                                       <  2:1043592/141‑1 (MQ=255)
                                                                                          gtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaac                                                    >  2:569713/1‑141 (MQ=255)
                                                                                          gtATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGaacaac                                                    >  1:281941/1‑141 (MQ=255)
                                                                                              tttttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGGGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCa                                                <  1:1103664/141‑1 (MQ=255)
                                                                                                tttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAg                                              <  1:1017243/141‑1 (MQ=255)
                                                                                                tttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAg                                              <  1:706136/141‑1 (MQ=255)
                                                                                                tttATAATCAAACTAAGCAATAACTTTTAGTTCTGTGTGATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAg                                              <  1:133841/141‑1 (MQ=255)
                                                                                                 ttATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGc                                             >  2:170422/1‑141 (MQ=255)
                                                                                                    aaaTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGa                                          <  1:859463/140‑1 (MQ=255)
                                                                                                      aTCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGCTGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGaaa                                        >  2:1063986/1‑141 (MQ=255)
                                                                                                       tCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGaaaa                                       <  1:355061/141‑1 (MQ=255)
                                                                                                       tCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGaaaa                                       <  1:1039530/141‑1 (MQ=255)
                                                                                                       acAAACTAATCAACGACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGaaaa                                       <  1:190119/140‑1 (MQ=255)
                                                                                                         aaaCTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAAtt                                     >  2:598033/1‑141 (MQ=255)
                                                                                                            cTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTgtt                                  >  2:315420/1‑139 (MQ=255)
                                                                                                            cTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGt                                  <  2:707602/141‑1 (MQ=255)
                                                                                                              aaTCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTg                                <  2:622542/141‑1 (MQ=255)
                                                                                                                tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCg                              >  2:945213/1‑141 (MQ=255)
                                                                                                                tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCg                              >  1:662580/1‑141 (MQ=255)
                                                                                                                tCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCg                              <  2:716853/141‑1 (MQ=255)
                                                                                                                tCAACAACTTTTAGGTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCg                              <  2:423484/141‑1 (MQ=255)
                                                                                                                tCAACAACTATTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGGAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGATGAAAATTGGTTGCg                              >  2:534347/1‑141 (MQ=255)
                                                                                                                       cTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTaataat                       <  2:210050/141‑1 (MQ=255)
                                                                                                                         tttAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAAGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAatat                     <  2:751797/141‑1 (MQ=255)
                                                                                                                         tttAGTTCTGTGTGATGGGATGGGTGGCCATAGAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAatat                     <  2:751825/141‑1 (MQ=255)
                                                                                                                         tttAGGTCTGTGTGATGGTATGGGGTGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAatat                     <  2:476117/141‑1 (MQ=255)
                                                                                                                           tAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGTTTGCGgaataaaaaaa                   >  1:447518/1‑136 (MQ=255)
                                                                                                                           tAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATaaaa                   >  2:3694/1‑141 (MQ=255)
                                                                                                                           tAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATaa                   >  1:338966/1‑141 (MQ=255)
                                                                                                                            aGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATaaa                  >  1:368872/1‑141 (MQ=255)
                                                                                                                            aGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATaaa                  >  1:469735/1‑141 (MQ=255)
                                                                                                                             gTTCTGTGTGATGGTCTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATaaaa                 >  1:394070/1‑141 (MQ=255)
                                                                                                                             gTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATaaaa                 >  1:1087394/1‑141 (MQ=255)
                                                                                                                                 tgtgtgATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGata             >  1:335747/1‑141 (MQ=255)
                                                                                                                                       aTGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAAtt       >  1:1001652/1‑141 (MQ=255)
                                                                                                                                       aTGGTATGGGGGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAAtt       <  2:44512/141‑1 (MQ=255)
                                                                                                                                          gTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAGATTTGTTGCGGATCAGTCGAACTCCCGTTTTGAAGCGGATGAGCTAATAGAACGACATCAAGCTGGAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTc    <  2:876410/141‑1 (MQ=255)
                                                                                                                                            aTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTCAg  >  2:1075535/1‑141 (MQ=255)
                                                                                                                                            aTGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAAATTTCAg  >  2:161655/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GACAAGTAGAAACGAGGTAAAGACAAATGCTAGAGGCACAATTTTTTACTGATACTGGACAACATAGAGATAAGAATGAAGATGCGGGTGGTATTTTTTATAATCAAACTAATCAACAACTTTTAGTTCTGTGTGATGGTATGGGTGGCCATAAAGCAGGAGAAGTTGCAAGTAAATTTGTTACAGATGAGTTGAAATCCCGTTTTGAAGCGGAAAATCTTATAGAACAACATCAAGCTGAAAATTGGTTGCGTAATAATATAAAAGATATAAATTTTCAG  >  USA300TCH1516_ALE/1231330‑1231610

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: