Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,926,533 A→G 5.6% K114E (AAA→GAA)  USA300HOU_RS09555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,926,5330AG 5.6% 211.6 / 3.6 71K114E (AAA→GAA) USA300HOU_RS09555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/40);  new base G (2/2);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

TAGAAGACTTTGCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCTAAATATAATGAA  >  USA300TCH1516_ALE/1926393‑1926671
                                                                                                                                            |                                                                                                                                          
tAGAAGACTTTGCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAgaga                                                                                                                                            <  1:189881/141‑1 (MQ=255)
tAGAAGACTTGGCGGCATTTAATTTATTAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAgaga                                                                                                                                            <  2:632185/141‑1 (MQ=255)
 agaagcGTTTGCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGaa                                                                                                                                           <  2:659777/141‑1 (MQ=255)
 agaagaCTTTGCGGCATTTAATTGATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGaa                                                                                                                                           <  1:304994/141‑1 (MQ=255)
       cTTTGCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGAc                                                                                                                                     <  1:314220/141‑1 (MQ=255)
           gCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACatat                                                                                                                                 >  2:492484/1‑141 (MQ=255)
               tatttaatttaTCCATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGTCATATaaaa                                                                                                                             <  2:140803/140‑1 (MQ=255)
               caattaatttaTCAATTGAAGACAAGTTTACAGTCTGTCTACAAAGAACTTAAGAACTAATAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACAtatatta                                                                                                                             <  2:130165/138‑4 (MQ=255)
               cATTTAATTTATAAATGGAAGAAACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAAGTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATaaaa                                                                                                                             <  2:63874/141‑1 (MQ=255)
                   taatttaTCAATTGAAGACACAGTAAGATGCTGTTTACAAAGAAATCAAGCACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCg                                                                                                                         <  2:448277/141‑1 (MQ=255)
                         aTCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTa                                                                                                                   <  1:414481/141‑1 (MQ=255)
                                    acacAGTTACAGTCTGTTTACAAAGAACTGAGGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGCCATTTAAAAGTCGAATGTAACTACGTGAGt                                                                                                        >  1:451671/1‑141 (MQ=255)
                                    acacAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGt                                                                                                        >  1:967285/1‑141 (MQ=255)
                                           tccAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGAtt                                                                                                 <  1:536936/139‑1 (MQ=255)
                                           tACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGAtt                                                                                                 <  2:258856/141‑1 (MQ=255)
                                            aCAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGAAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGAttt                                                                                                >  1:505460/1‑141 (MQ=255)
                                                gcTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTAtt                                                                                            <  1:281569/140‑1 (MQ=255)
                                                 taagttACAAGGAAGCGAAGAACAATTAAAACAATACAAACATTTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTc                                                                                           <  2:156256/137‑1 (MQ=255)
                                                 ccgTTTAACAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTc                                                                                           <  1:971323/139‑1 (MQ=255)
                                                 cTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTc                                                                                           <  1:325275/141‑1 (MQ=255)
                                                   gattAAAAAGAAGTGAAGATCCATTAAAACAATACAAAAATGAGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTGTTGATTTTATTCtt                                                                                         <  1:722683/139‑1 (MQ=255)
                                                   gTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCtt                                                                                         >  1:744124/1‑141 (MQ=255)
                                                    tttACCAAGAACTGAAGAACTATGAAAACGATACAAAAATGTGGAATTCAATGACAAAATATTGATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCttt                                                                                        <  1:823709/141‑1 (MQ=255)
                                                       aCAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGGAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTa                                                                                     >  1:210975/1‑141 (MQ=255)
                                                       aCAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGGAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTa                                                                                     >  2:396488/1‑141 (MQ=255)
                                                          aaGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAAt                                                                                  <  2:264319/141‑1 (MQ=255)
                                                          aaGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAAt                                                                                  <  2:77242/141‑1 (MQ=255)
                                                          aaGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAAt                                                                                  <  2:264338/141‑1 (MQ=255)
                                                            gAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTa                                                                                >  1:278341/1‑141 (MQ=255)
                                                              aCTGAAGAACTATTAAAACCAAACAAAAATGTGGAAATCAATGACAAAATAGTAATCAATTGGCGTATTATACAAGATAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACa                                                                              <  2:60730/141‑1 (MQ=255)
                                                              aCTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACa                                                                              <  1:30229/141‑1 (MQ=255)
                                                                tGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAg                                                                            >  1:336445/1‑141 (MQ=255)
                                                                tGAAGAACTATTAAAACAAGACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAg                                                                            >  1:608575/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCa                                                                          >  1:252249/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCa                                                                          >  1:252161/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCa                                                                          >  1:252120/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCa                                                                          >  2:1020903/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCa                                                                          >  1:136111/1‑141 (MQ=255)
                                                                  aagaaCTATTAAAACAATACAAAAATGTGAAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGGCGAATGTAACTACGTGAGTATTGCTTTTAGTCTTTTTAAATTAAAAGCa                                                                          >  1:794837/1‑141 (MQ=255)
                                                                    gaaCTATTAAAACAATACAAAAATGTGAAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAAAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGCTTTTATTCTTTGTAAATTACAAGCAtt                                                                        >  2:107199/1‑141 (MQ=255)
                                                                         aTTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCata                                                                   >  1:772885/1‑141 (MQ=255)
                                                                         aTTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCata                                                                   >  1:788701/1‑141 (MQ=255)
                                                                         aTTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCata                                                                   >  1:152404/1‑141 (MQ=255)
                                                                                cAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGCTTTTATTCTTTGTAAATTACCACCATTTCATATTATaaa                                                            >  1:1068969/1‑141 (MQ=255)
                                                                                   tACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGtt                                                         <  2:744124/141‑1 (MQ=255)
                                                                                    aCAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATCGAGTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGttt                                                        <  1:485443/141‑1 (MQ=255)
                                                                                     aaaaaaTGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTg                                                       <  2:552040/140‑1 (MQ=255)
                                                                                      aaaaaTGTGGAAGGCANTGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGa                                                      <  1:937866/141‑1 (MQ=255)
                                                                                              ggAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGACATCGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGt                                              >  1:699996/1‑141 (MQ=255)
                                                                                              ggAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGt                                              >  1:654123/1‑141 (MQ=255)
                                                                                                  tttaaTGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTatat                                          <  2:84257/138‑1 (MQ=255)
                                                                                                  ttCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTatat                                          <  2:621277/141‑1 (MQ=255)
                                                                                                        gACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAAt                                    <  2:959231/141‑1 (MQ=255)
                                                                                                         aCAAAATATTAATCAATTTGCGTATTATACAAGAAAAATAGACTTATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATg                                   <  2:1003839/141‑1 (MQ=255)
                                                                                                         aCAAAATATTAATCAATTGGCGTATTATACAAGAGGAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATg                                   <  1:396488/141‑1 (MQ=255)
                                                                                                         aCAAAATATTAATCAATTGGCGTATTATACAAGAGGAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATg                                   <  2:210975/141‑1 (MQ=255)
                                                                                                         aCAAAATATTAATCAATTGGCGGATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATg                                   <  1:764979/141‑1 (MQ=255)
                                                                                                          cAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATgg                                  <  2:981684/141‑1 (MQ=255)
                                                                                                          aaaaaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATgg                                  <  2:604880/140‑1 (MQ=255)
                                                                                                            aaaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGgag                                <  2:1090196/141‑1 (MQ=255)
                                                                                                            aaaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGgag                                <  2:325664/141‑1 (MQ=255)
                                                                                                             aaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGgaga                               >  2:635694/1‑141 (MQ=255)
                                                                                                             aaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGgaga                               >  2:12700/1‑141 (MQ=255)
                                                                                                             aaTATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATAGTATAAAGTTTGAAAAGAGGTATATTGAAATGaaga                               >  2:868814/1‑141 (MQ=255)
                                                                                                                 ttAATCAATTGGGGTATTAAACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGaaaaa                           <  1:431102/141‑1 (MQ=255)
                                                                                                                 ttAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATGTTGAAATGGAGaaaaa                           <  1:547117/141‑1 (MQ=255)
                                                                                                                   aaTCACTTCGCGTATTATCAAAAAGAAATAGACATATAAAACTCGAATGTAACTACGTGAGTATTGGTTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTTAAAAGAGGTATATTGAAATGGAGAAAAATg                         <  2:554673/141‑1 (MQ=255)
                                                                                                                   aaTCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATg                         <  2:1066400/141‑1 (MQ=255)
                                                                                                                     acAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGaa                       <  1:143343/140‑1 (MQ=255)
                                                                                                                       aaTTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAata                     >  2:854209/1‑141 (MQ=255)
                                                                                                                        aTTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAatat                    >  2:1049621/1‑141 (MQ=255)
                                                                                                                        aTTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTAAGTGAGTATTAATGGTATTCTTTCTAAGCTTGTCGCTGTTCAGAGTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAatat                    >  2:502465/1‑141 (MQ=255)
                                                                                                                              gtgttatATAAGAGAAATAGACATATAAAAGTCGAATATAACTACGTGAGTATTGATTTCATTCTTTGTAAATTACAAGCATGTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCt              <  2:59014/138‑1 (MQ=255)
                                                                                                                                 ttatACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCTaaa           <  2:278341/141‑1 (MQ=255)
                                                                                                                                        agagAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCTAAATATAATg    <  1:515843/141‑1 (MQ=255)
                                                                                                                                          agaAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCTAAATATAATGaa  <  1:1020903/141‑1 (MQ=255)
                                                                                                                                          agaAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTAAAATGGAGAAAAATGAATATATAGCTAAATATAATGaa  <  2:136111/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TAGAAGACTTTGCGGCATTTAATTTATCAATTGAAGACACAGTTACAGTCTGTTTACAAAGAACTGAAGAACTATTAAAACAATACAAAAATGTGGAATTCAATGACAAAATATTAATCAATTGGCGTATTATACAAGAGAAATAGACATATAAAAGTCGAATGTAACTACGTGAGTATTGATTTTATTCTTTGTAAATTACAAGCATTTCATATTATAAAGTTTGAAAAGAGGTATATTGAAATGGAGAAAAATGAATATATAGCTAAATATAATGAA  >  USA300TCH1516_ALE/1926393‑1926671

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