Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 19,264 G→A G307D (GGC→GAC)  gdpP → Cyclic‑di‑AMP phosphodiesterase GdpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE19,2640GA96.2% 301.9 / NA 107G307D (GGC→GAC) gdpPCyclic‑di‑AMP phosphodiesterase GdpP
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (59/41);  minor base C (2/2);  total (62/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

AAGTATCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATT  >  USA300TCH1516_ALE/19127‑19397
                                                                                                                                         |                                                                                                                                     
aaGTATCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACgg                                                                                                                                    >  1:159803/1‑141 (MQ=255)
    aTCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGTTAAg                                                                                                                                >  2:285914/1‑141 (MQ=255)
    aTCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAg                                                                                                                                >  1:888479/1‑141 (MQ=255)
     tCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGa                                                                                                                               >  1:647857/1‑141 (MQ=255)
      cGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAAGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTGTCTATGGCGGTGAgcc                                                                                                                              >  2:343700/1‑139 (MQ=255)
       ggTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACt                                                                                                                             >  2:717918/1‑141 (MQ=255)
       ggTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATCAATGGTAATGTGGGTTTCTATGACGGTAAGACt                                                                                                                             >  2:551386/1‑141 (MQ=255)
           ttGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGAcc                                                                                                                         <  1:895095/141‑1 (MQ=255)
           ttGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGAcc                                                                                                                         <  2:854199/141‑1 (MQ=255)
              gTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTCAGACTGACCCGa                                                                                                                      >  1:1667/1‑141 (MQ=255)
                  aGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATggg                                                                                                                  >  1:380825/1‑140 (MQ=255)
                    gTACTGAAAATTTAATCGACTTAGGTGAATTATCATAATCAGGCCTAGACTTAGCATTAGGTCGCGGTGGCGACCAGGTTGCAATTAAAAGTATTAATGGTAGTGTGCGTTTCTATGACGGTTAGACTGCGCCGGTGgagg                                                                                                                >  2:687304/1‑140 (MQ=255)
                    gTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGggga                                                                                                                >  1:913932/1‑141 (MQ=255)
                    gTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGgaga                                                                                                                >  2:761586/1‑141 (MQ=255)
                    gTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGgaga                                                                                                                <  2:808870/141‑1 (MQ=255)
                      aCTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGaaa                                                                                                              <  2:179220/141‑1 (MQ=255)
                      aCTGAAAATTTAATCGACTTAGGTGAATTATCAAAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGaaa                                                                                                              <  2:832765/141‑1 (MQ=255)
                       ctcaaaaTTAAATCGATTTATGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAAc                                                                                                             <  1:169695/138‑1 (MQ=255)
                        aaacaaTTAACTCGACTTACGTGAATTATCACAATCAGGCCTAGACTTATCATTAGGACGCGGTGGCTACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACg                                                                                                            <  2:943129/137‑1 (MQ=255)
                            aaTTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACt                                                                                                        <  2:758983/141‑1 (MQ=255)
                            aaTTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACt                                                                                                        <  2:833078/141‑1 (MQ=255)
                              tttAATCGACTTAGGTGAATTAGCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCg                                                                                                      <  1:82016/141‑1 (MQ=255)
                               ttAATCGGGTTAGGTGAATTATCACGATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCgt                                                                                                     <  2:308364/141‑1 (MQ=255)
                                  aTCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTa                                                                                                  >  1:187515/1‑141 (MQ=255)
                                  aTCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTa                                                                                                  >  2:648330/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGCCTAAGTTGCAATTAAAAGTATTAATGGTATTGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGGGTaa                                                                                                 >  1:867953/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTaa                                                                                                 >  2:975876/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTaa                                                                                                 >  2:830753/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTaa                                                                                                 >  1:211606/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGAATGACCCGATGGAGAAACGTACTCGTGTaa                                                                                                 >  1:894125/1‑141 (MQ=255)
                                   tCGACTTAGGTGAATAATCGCAATCAGGCCTAGATTTAGCATTAGGACACGGTGGCGACAACGTAGCTATTATATGTATTAATAGTAATGAGCGTTTCTATGACGGTAAGGCCGACCCGGTGGAGAAACGTACTCGTGTaa                                                                                                 >  1:738242/1‑141 (MQ=255)
                                       cTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTCAGAGc                                                                                             >  1:801746/1‑141 (MQ=255)
                                       cTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGc                                                                                             >  1:458845/1‑141 (MQ=255)
                                              gAATTATCACGATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGtgtg                                                                                      >  1:716776/1‑141 (MQ=255)
                                                 ttATCACAATCAGGCATAGACTTAGCATTAGGACGCGATGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAGGACTGACCCGATGGAGAAACGTACTAGTGTAAGAGCACGTGTGAtc                                                                                   <  1:695981/141‑1 (MQ=255)
                                                     cacaTTCAGGCCTTGCATTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGTCTGTCCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTcac                                                                               <  2:901737/141‑1 (MQ=255)
                                                     cacaATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGGAACGTACTCGTGTAAGAGCACGTGTGATCTcac                                                                               >  2:457067/1‑141 (MQ=255)
                                                     cacaATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTcac                                                                               >  2:721984/1‑141 (MQ=255)
                                                     cacaATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTcac                                                                               >  1:654602/1‑141 (MQ=255)
                                                       caATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGGGTGGTCTcccag                                                                             >  1:877243/1‑140 (MQ=255)
                                                       caATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGGGTAAGAGCACGTGTGATCTcacag                                                                             >  1:279747/1‑140 (MQ=255)
                                                       caATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGCTTCTATCCCGTTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGGGCACGTGTGATCGCACAt                                                                             >  1:718802/1‑141 (MQ=255)
                                                        aaTCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGATCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTCGGACGGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATg                                                                            <  1:384023/141‑1 (MQ=255)
                                                        aaTCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTTCAAGTATAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTATTGGTGGAGGGGCACGTGTGATCTCACATg                                                                            >  1:433929/1‑141 (MQ=255)
                                                        aaTCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGGGCACGTGTGATCCCACATg                                                                            >  2:957526/1‑141 (MQ=255)
                                                        aaTCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTAACATg                                                                            >  2:480617/1‑141 (MQ=255)
                                                        aaTCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACGCGTGTAAAAGCACGTGTGATCTcacaag                                                                            >  2:378611/1‑139 (MQ=255)
                                                          tCATGCCTAGACGTAGCATAAGAATTGTGTTGCTACCAAGTTGCGATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCg                                                                          <  2:690845/141‑1 (MQ=255)
                                                                cTAGACGTGGCATTAGGGCGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTaaaa                                                                    <  2:25329/141‑1 (MQ=255)
                                                                 tagaaaaaGCATTAGGACGCGGTGGCGAACCAGTTGAATTTAACAGCGTTAATGGCAATGTGCGTTTCTATGACGGTAAGATTGACTCGAAGTAGAATCGTACTCGTGTAAGAGCACGCGTGATCTCACATGCGTTAAAAg                                                                   <  1:573261/134‑1 (MQ=255)
                                                                   gACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGat                                                                 >  1:291693/1‑141 (MQ=255)
                                                                       tAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATcc                                                             <  1:496940/141‑1 (MQ=255)
                                                                       tAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGGTAAAAAATATcc                                                             >  1:375936/1‑141 (MQ=255)
                                                                                gACGGGGTGGCGACTAATTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAgg                                                    <  1:784108/141‑1 (MQ=255)
                                                                                  cgcgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAgggg                                                  >  2:260501/1‑140 (MQ=255)
                                                                                  cgcgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGt                                                  <  1:834663/141‑1 (MQ=255)
                                                                                    cgGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGCAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGa                                                >  2:890018/1‑141 (MQ=255)
                                                                                     ggTGGCGACCAAGTTGCAATGAAAAGTATGAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGAc                                               <  1:528650/141‑1 (MQ=255)
                                                                                        ggCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGGCaaa                                            >  1:907594/1‑141 (MQ=255)
                                                                                        ggCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACaaa                                            >  2:399362/1‑141 (MQ=255)
                                                                                         gCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAg                                           >  1:92115/1‑141 (MQ=255)
                                                                                         gCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAg                                           >  1:387594/1‑141 (MQ=255)
                                                                                          cGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCGCGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGt                                          <  1:830753/141‑1 (MQ=255)
                                                                                             ccAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACACGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGGCAAAGTCAt                                       >  2:630714/1‑141 (MQ=255)
                                                                                              cAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCAtt                                      <  2:321982/141‑1 (MQ=255)
                                                                                               aaGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGACGGTAAGTCTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTa                                     <  2:647857/141‑1 (MQ=255)
                                                                                                            aaatATTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATaa                        <  2:320286/138‑1 (MQ=255)
                                                                                                            aaGTATTAAGGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATaa                        <  2:320295/141‑1 (MQ=255)
                                                                                                                aTTAATGGTAATGTGCGTTTCTATGCCGGTAAGACTGACCCGATGGAGAAACGTACACGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATACACga                    >  1:406768/1‑140 (MQ=255)
                                                                                                                aTTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTTCTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGt                    >  1:658216/1‑141 (MQ=255)
                                                                                                                aTTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGGTCTCACATGCGTTAAACGATATCCTTGCAGAGGGTGACAAAGTCATTAGCATGGGACATAAACGt                    >  1:692925/1‑141 (MQ=255)
                                                                                                                aTTAATGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGt                    >  2:689045/1‑141 (MQ=255)
                                                                                                                  tactGGCAATGTGCCTTCCTACGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTCAAAAGATATCCTTGCCGAGGGTGACAAAGTCATTATCATGGGACATAAACGTcc                  <  1:442768/138‑1 (MQ=255)
                                                                                                                   aaTGGTAGTGTGCGTTTCTATGACGGTAAGACTGACCCGAAGGAGAAACGTACTCGTGTCAGAGCACGTGTGATCTCACATGCGTTAAGAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCt                 <  2:880268/141‑1 (MQ=255)
                                                                                                                   aaTGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGAAGAGGGTGACAAAGTCATAATCATGGACCATAAACGTcag                 >  1:695727/1‑139 (MQ=255)
                                                                                                                     tGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGa               <  2:187515/141‑1 (MQ=255)
                                                                                                                     tGGTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGa               <  1:42633/141‑1 (MQ=255)
                                                                                                                      ggTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGAc              <  2:663141/141‑1 (MQ=255)
                                                                                                                       ggaaTGCGCGTTTATATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt             <  1:914653/139‑1 (MQ=255)
                                                                                                                       gTAATGTGCGTTTCTATGTCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt             <  1:91330/141‑1 (MQ=255)
                                                                                                                       gTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt             <  1:568368/141‑1 (MQ=255)
                                                                                                                       gTAATGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt             <  1:79597/141‑1 (MQ=255)
                                                                                                                       gTAATGCGCCTTTCTATGCCCGTCAGACTGACCCGATGGAGAAACGTACGCGGTTGAGAGAACTTGTGATCTCACATGCGTTAAAAGAGAGCTTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACt             <  1:375288/141‑1 (MQ=255)
                                                                                                                          aTGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAg          <  1:218393/141‑1 (MQ=255)
                                                                                                                          aTGTGCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAg          <  2:888479/141‑1 (MQ=255)
                                                                                                                           tgtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGa         >  1:950169/1‑141 (MQ=255)
                                                                                                                            gtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGt                                                  >  2:896201/1‑99 (MQ=255)
                                                                                                                            gtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGt                                                  <  1:896201/99‑1 (MQ=255)
                                                                                                                            gtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGAGAAAGTCGTTATCATGGGACATAAAGGTCCTGAGTTAgtt        >  1:975662/1‑139 (MQ=255)
                                                                                                                            gtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTGGAt        >  1:340896/1‑141 (MQ=255)
                                                                                                                            gtgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGAt        >  1:590897/1‑141 (MQ=255)
                                                                                                                             tgCGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGGCCTGACTTAGATg       >  1:73507/1‑141 (MQ=255)
                                                                                                                               cGTTTCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTAATTATCATGGGACCTAAAAGTCCTGACTTAGATGCa     >  2:358726/1‑141 (MQ=255)
                                                                                                                               cGTTGCTATGCTGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCa     <  2:83505/141‑1 (MQ=255)
                                                                                                                               agtctCTATGACGCTAAGACTGACGTGATGGAGAAAGGTACTCGTGTAAGAGAACGTGTGACCTCACATGCGTTAAAAGATCTCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATTCa     <  2:309709/137‑1 (MQ=255)
                                                                                                                                gTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCaa    <  1:286883/141‑1 (MQ=255)
                                                                                                                                gTTTCTATGACGGTAAGACTGACCCGATAGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCaa    <  1:687752/141‑1 (MQ=255)
                                                                                                                                 tttCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGGGACAAAGTCATTATCATGGGACGTAAACGACCTGACGTAGATGCAAt   >  2:371749/1‑141 (MQ=255)
                                                                                                                                  ttCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAAtt  >  1:47345/1‑141 (MQ=255)
                                                                                                                                  ttCTATGACGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAAtt  >  2:32839/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                     
AAGTATCGGTGTTGGTGAAGGTACTGAAAATTTAATCGACTTAGGTGAATTATCACAATCAGGCCTAGACTTAGCATTAGGACGCGGTGGCGACCAAGTTGCAATTAAAAGTATTAATGGTAATGTGCGTTTCTATGGCGGTAAGACTGACCCGATGGAGAAACGTACTCGTGTAAGAGCACGTGTGATCTCACATGCGTTAAAAGATATCCTTGCAGAGGGTGACAAAGTCATTATCATGGGACATAAACGTCCTGACTTAGATGCAATT  >  USA300TCH1516_ALE/19127‑19397

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: