Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 2,671,907 G→A H46Y (CAC→TAC)  USA300HOU_RS13760 ← Putative esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,671,9070GA98.7% 243.3 / NA 75H46Y (CAC→TAC) USA300HOU_RS13760Putative esterase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (43/30);  minor base C (0/1);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTGCTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTGTTGTCA  >  USA300TCH1516_ALE/2671775‑2672041
                                                                                                                                    |                                                                                                                                      
tCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATtaa                                                                                                                                >  1:555347/1‑141 (MQ=255)
        gcttTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGGACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCt                                                                                                                        <  2:283358/141‑1 (MQ=255)
             aaCTAAATAAATACGACTTACTGGTTTGTCATGTATTGTGACAGGTACTTTATTCTCTGCTTCCGTTCTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTtct                                                                                                                   <  1:278976/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTTAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                               <  1:868866/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                               <  1:705749/141‑1 (MQ=255)
                 aaGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCaaa                                                                                                               <  1:648825/141‑1 (MQ=255)
                  aGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAAt                                                                                                              <  1:414842/141‑1 (MQ=255)
                    tAAATACGTCCTATTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCt                                                                                                            <  2:810566/141‑1 (MQ=255)
                    tAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCt                                                                                                            <  1:13672/141‑1 (MQ=255)
                    tAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCt                                                                                                            <  1:815526/141‑1 (MQ=255)
                     aaaTACGACCTACTGGTTCGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCtt                                                                                                           <  1:317121/141‑1 (MQ=255)
                      aaTACGACCTAGTGGTTTGTCATGTATTGTGACAGTTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGTTTAATAAACGCTGTTCTCAAATCTTc                                                                                                          <  2:107396/141‑1 (MQ=255)
                      aaTACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTc                                                                                                          >  1:368076/1‑141 (MQ=255)
                      aaTACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCAGGGACTAGTTAATAAAGGCTGTTCTCAAATCTTc                                                                                                          >  2:251162/1‑141 (MQ=255)
                       aTACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCg                                                                                                         >  2:692543/1‑141 (MQ=255)
                       aTACGACCTACTGGTTTGTCATGTATTGTGACAGGTACATTATACGCTGCTTCCGATTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAGTAAACGCTGTTCTCAAATCTTca                                                                                                         >  2:693679/1‑140 (MQ=255)
                         aCGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCGCAAATCTTCGaa                                                                                                       >  2:165516/1‑141 (MQ=255)
                            aCCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGCTTAATAACCGCTGTTCTCAAATCTTCAAAccc                                                                                                    >  2:552158/1‑140 (MQ=255)
                            aCCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCTAACCa                                                                                                    >  1:306523/1‑141 (MQ=255)
                            aCCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCa                                                                                                    >  2:150115/1‑141 (MQ=255)
                             ccTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCac                                                                                                   >  2:457831/1‑140 (MQ=255)
                             ccTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCAt                                                                                                   >  2:128428/1‑141 (MQ=255)
                               tACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCGAATCTTCGAACCATCt                                                                                                 <  1:472408/141‑1 (MQ=255)
                               tACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCt                                                                                                 >  2:141432/1‑141 (MQ=255)
                                 cTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTaa                                                                                               >  1:229478/1‑141 (MQ=255)
                                 cTGGTTTGTCATGTAGTGTGACAGGAACTTTATTCTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTaa                                                                                               <  2:415467/141‑1 (MQ=255)
                                  tGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAat                                                                                              >  1:919830/1‑141 (MQ=255)
                                      ttGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATaa                                                                                          <  2:858994/141‑1 (MQ=255)
                                                    gACAGGTACTTTATACTCTACTTCAACTGGCAGAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAAc                                                                            >  1:678154/1‑141 (MQ=255)
                                                        ggTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGaca                                                                        >  2:389270/1‑141 (MQ=255)
                                                            cTTTATACTCTGCTTCCGTATTCATAAGTATAGGCGAAATGCGTTGATAGATCATTGTTGAGTAATTGATGTGCAGATTAATAGACGCTGGTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTGACGACACCCa                                                                    >  1:688942/1‑141 (MQ=255)
                                                               tatcCTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGACCATTGTTGAGTAATTAATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATAAAAACGATATTACTAAAGACACCCATCg                                                                 >  2:112401/1‑141 (MQ=255)
                                                                ataCTCTGCCTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGc                                                                <  1:501315/141‑1 (MQ=255)
                                                                      tGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAAt                                                          <  1:356241/141‑1 (MQ=255)
                                                                       gCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATg                                                         >  1:359798/1‑141 (MQ=255)
                                                                            cGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTa                                                    >  1:514513/1‑141 (MQ=255)
                                                                             gTTTTCATAAGGATAGGTGAAATGCCTTGATTGATCGTTGTTGAGTAAATAACACACTGTGTAAGAAAAGCTGTTATCACATCTTCTAACCATCTAAAATAACCGATATTACTAACAACACACACCGCACCTATGCCGTaa                                                   >  2:166806/1‑141 (MQ=255)
                                                                               tttCATAAGTATATGTGAAATGCCTTGTTTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAAcc                                                 <  1:392397/141‑1 (MQ=255)
                                                                                 tCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGGCGTAACCAt                                               >  1:92135/1‑141 (MQ=255)
                                                                                 tCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAAGAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAAccgt                                               >  1:156513/1‑139 (MQ=255)
                                                                                       gTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACtt                                         >  2:458774/1‑141 (MQ=255)
                                                                                       gTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACtt                                         >  2:192744/1‑141 (MQ=255)
                                                                                             gTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAttt                                   <  1:33918/141‑1 (MQ=255)
                                                                                              tGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttc                                  <  2:355541/141‑1 (MQ=255)
                                                                                              tGAAACGCCGTGATTGATCATTGTTGAGTAATTCATGTACAGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttc                                  <  2:356336/141‑1 (MQ=255)
                                                                                                  aTGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttctttc                              >  1:753740/1‑141 (MQ=255)
                                                                                                  aTGCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGAtttctttc                              <  1:899508/141‑1 (MQ=255)
                                                                                                    gCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATGTTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTtctc                            >  1:25793/1‑141 (MQ=255)
                                                                                                    gCCTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACAAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTtctc                            >  1:714758/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCGTTGTTGAGCAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                           >  1:476607/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAGCGACAACCATAGCATCAATCTCGTAACCATTCACTTCGATTTCTTTctct                           >  2:72551/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGCTTTCTTTctct                           >  1:868409/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                           >  1:843016/1‑141 (MQ=255)
                                                                                                     ccTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctct                           >  1:362113/1‑141 (MQ=255)
                                                                                                      tttgtttgaTCATTGTTTTGTCATTCACGTACTGATTAATAACAGCTGTTTTCAGAGCGTCGATCCATCTAATATATACGATATTATGTACGACTCCCATCGCTTCAAGGTCGTAATCATTCACTTCGATTTCTTTCTCTg                          <  2:130557/140‑1 (MQ=255)
                                                                                                      cTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTctctc                          >  2:603581/1‑140 (MQ=255)
                                                                                                      cTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTg                          >  1:185928/1‑141 (MQ=255)
                                                                                                      cTTGATTGATCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTg                          >  1:138850/1‑141 (MQ=255)
                                                                                                       ttgattgaTCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGc                         <  2:833324/141‑1 (MQ=255)
                                                                                                       ttgattgaTCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGc                         >  1:113239/1‑141 (MQ=255)
                                                                                                       ttgattgaTCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGc                         >  1:188722/1‑141 (MQ=255)
                                                                                                              atCATTGTTGAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTTCAACTAAt                  >  2:71499/1‑141 (MQ=255)
                                                                                                                cATTGTTGAGTAATTCATGTACTGATAAATAAACGCTGATCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCACTGTCGTAAGCATCCACTTCGATTTCTTTCTCTGCAACTAAtga                >  2:323606/1‑139 (MQ=255)
                                                                                                                   tgttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGAcc             <  1:697887/141‑1 (MQ=255)
                                                                                                                    gttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCt            <  1:329696/141‑1 (MQ=255)
                                                                                                                     ttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCtt           <  1:141432/141‑1 (MQ=255)
                                                                                                                     ttgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCtt           <  1:185973/141‑1 (MQ=255)
                                                                                                                      tgAGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACACATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTc          >  2:689893/1‑141 (MQ=255)
                                                                                                                       gAGTAAGTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCt         >  2:508534/1‑141 (MQ=255)
                                                                                                                        aGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCtg        <  1:251162/141‑1 (MQ=255)
                                                                                                                        aGTAATTCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCtg        <  2:368076/141‑1 (MQ=255)
                                                                                                                           aaTTCATGTCCTGATTAATAAACGCTGTTACCAAAGCTTCGAACCATCTAATAGAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCtgttg     <  2:548273/141‑1 (MQ=255)
                                                                                                                             ttCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTTTTGTc   >  2:378415/1‑141 (MQ=255)
                                                                                                                             ttCATGTACTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTGTTGTc   >  2:533278/1‑141 (MQ=255)
                                                                                                                              tCATGTACTAATTAATACCCGCTGCCCGCAAATCTTCGTACCATCTAATATAAACGATCTTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCTATTTCTTTCTCTGCAACTAATACACCTTCTGTTGTCa  <  2:688942/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TCTTGCTTGCTTTAACTAAGTAAATACGACCTACTGGTTTGTCATGTATTGTGACAGGTACTTTATACTCTGCTTCCGTTTTCATAAGTATAGGTGAAATGCCTTGATTGATCATTGTTGAGTAATTCATGTGCTGATTAATAAACGCTGTTCTCAAATCTTCGAACCATCTAATATAAACGATATTACTAACGACACCCATCGCATCAATGTCGTAACCATTCACTTCGATTTCTTTCTCTGCAACTAATAGACCTTCTGTTGTCA  >  USA300TCH1516_ALE/2671775‑2672041

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