Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,470,089 G→A 6.5% G582G (GGC→GGT)  emrB ← Multidrug export protein EmrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,470,0890GA 6.5% 171.9 / 2.6 64G582G (GGC→GGT) emrBMultidrug export protein EmrB
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CTTTTGCTTTTTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGATTGCTGTA  >  USA300TCH1516_ALE/2469950‑2470227
                                                                                                                                            |                                                                                                                                          
cttttgcttttTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTCACACTTTCATTgcgc                                                                                                                                           >  1:213193/1‑141 (MQ=255)
   tttcttttTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCctt                                                                                                                                        <  1:470737/138‑1 (MQ=255)
    tgcttttTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCcttc                                                                                                                                       <  1:169569/141‑1 (MQ=255)
     gcttttTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCcttct                                                                                                                                      <  1:351100/141‑1 (MQ=255)
            tatcaCTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCCAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGc                                                                                                                               <  2:588128/141‑1 (MQ=255)
             cttccTTTCTAAAATTAAAGATAAAATTAAGTCGGTGGTTGTGATTGTTGTTGCATGTATAAATGGTTCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTCGTAACACTTTCATTGCGCCTTCTTGTCCGcc                                                                                                                              >  2:283741/5‑141 (MQ=255)
               taCTTTGTAAATATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCat                                                                                                                            <  2:292713/141‑1 (MQ=255)
                 cTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTACGAGTATATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCatat                                                                                                                          >  1:275412/1‑141 (MQ=255)
                      tAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGa                                                                                                                     >  1:496424/1‑141 (MQ=255)
                      tAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGa                                                                                                                     >  1:260417/1‑141 (MQ=255)
                          aaTAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCt                                                                                                                 <  1:221871/141‑1 (MQ=255)
                            aaaaCATAAAATTAAGGAGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTaa                                                                                                               <  1:126357/140‑1 (MQ=255)
                             aaaCAGAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAAt                                                                                                              <  2:608198/141‑1 (MQ=255)
                               aCATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTc                                                                                                            >  1:519540/1‑141 (MQ=255)
                                    aaaTTAAGGCTATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCa                                                                                                       <  1:95427/141‑1 (MQ=255)
                                     aaTTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGTCTCATTCATGATTATGAAAACCT‑CTAGTTTATTTACAAATTGTAGCAGCACTTTCATGGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCat                                                                                                      <  2:590379/141‑1 (MQ=255)
                                      aTTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCatt                                                                                                     >  2:687221/1‑140 (MQ=255)
                                       ttAAGGGGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACTTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTTTCCGCCATATTGTGATGCTAATTCATGCatat                                                                                                    <  2:259928/141‑1 (MQ=255)
                                            tcGAGGATGGGAAATAATGTTGAAACTATAACTACATCATTAAT‑ACCTTCAACAGTTGCTGGGTTATGTACAAATTGTAGTAACACTTTTATTGCTCCCTCTTGCCCGACATATTGTGATGTTATTTCATGCATATGGTCt                                                                                               <  2:59008/140‑1 (MQ=255)
                                                  atgCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTaca                                                                                         <  1:555537/141‑1 (MQ=255)
                                                                    aCTATTAATTCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCt                                                                       <  1:309404/141‑1 (MQ=255)
                                                                        taAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAAct                                                                   >  1:264620/1‑141 (MQ=255)
                                                                           aTGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCtt                                                                >  1:475577/1‑141 (MQ=255)
                                                                           aTGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAATAACGCTTTCATTGCGCCTTCTTGTCCGCAATATTGTGATGCTACTTCATGCATATGATATTGTACCAACGCATTAGTTTTATCTaagcatt                                                                >  2:13412/1‑136 (MQ=255)
                                                                                      aT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCt                                                     >  2:783841/1‑141 (MQ=255)
                                                                                       t‑ACCTTCAACCGTTGCTAGTTTATTTACACATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCACAAAAACTg                                                    >  1:269729/1‑141 (MQ=255)
                                                                                       t‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTg                                                    <  1:656448/141‑1 (MQ=255)
                                                                                          ccTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCACCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGAt                                                  >  2:654698/1‑141 (MQ=255)
                                                                                          ccTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCACCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGAt                                                  >  2:238345/1‑141 (MQ=255)
                                                                                          ccTTCAACCGTGGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGATTATGACATTGTACTGCATGTTTCGTTTTATCTAACTCTTCCCCAAAAGCTgtt                                                  >  2:379182/1‑139 (MQ=255)
                                                                                            ttCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGAGCCTTCTTGTCGGCCATATTGTGATGCTAATTCATGCATATGATATTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCGGATac                                                >  2:154816/1‑140 (MQ=255)
                                                                                            ttCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACAATTTCATTGCGCCTTCTTGTCCTCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTCACTCTTCCCCAAAAGCTGATaa                                                >  2:256775/1‑141 (MQ=255)
                                                                                             tCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTACTTCACCTACCTTGTCAGGTTCAACAGGGTTCGTTTTATCTAACTCTTCCCCAAAAGCTGGTAAg                                               >  1:233298/1‑141 (MQ=255)
                                                                                                tccGTTCGTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCTCTTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAgtgt                                            <  2:605399/140‑1 (MQ=255)
                                                                                                aCCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGTTAAgttt                                            >  1:53627/1‑139 (MQ=255)
                                                                                                  cGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGtgttg                                          >  1:52817/1‑141 (MQ=255)
                                                                                                     tGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTgtgt                                       <  2:20690/141‑1 (MQ=255)
                                                                                                      gCTAGTTTATTTACAAATTGTAGGAACACTTTCAGTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTa                                      <  1:343420/141‑1 (MQ=255)
                                                                                                       cTAGTTTATTTACAAAGTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAg                                     <  1:182073/141‑1 (MQ=255)
                                                                                                         aGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTGGtt                                   >  2:816225/1‑141 (MQ=255)
                                                                                                         aGTTTATTTACAAATTGTAGTAACACTTTCATGGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGtt                                   <  2:412207/141‑1 (MQ=255)
                                                                                                          gTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGttt                                  <  1:796453/141‑1 (MQ=255)
                                                                                                            ttatttaCCCATTGTAGTAACACTTTCATTGCGTCTTCTTGTCCTTTATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTtgt                                <  2:442586/141‑1 (MQ=255)
                                                                                                             tatttaCAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATCTTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGCATTATCTAACTCTTCGCCACAAGCTGATAAGTGTTGTGTAGTTtgtg                               >  1:482920/1‑141 (MQ=255)
                                                                                                              ttttaCTAATTGTAGTAACACTTTCATTGCGCCTTCTTGTACGCCATATTGTGATGCTAATTCATGTATATGATCTTGTACAACTGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTgtgt                              <  1:269754/140‑1 (MQ=255)
                                                                                                              atttaCAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTgtgt                              <  1:672523/141‑1 (MQ=255)
                                                                                                              atttaCAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTgtgt                              <  1:573454/141‑1 (MQ=255)
                                                                                                                ttaCAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGCTAAGTGTTGATTATTTtgaaaca                            <  2:598095/141‑6 (MQ=255)
                                                                                                                 taCGAATTGTAGGAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGtgttgt                           <  2:325146/141‑1 (MQ=255)
                                                                                                                 taCAAATTGTAGTAACACTTTCATTGCGCCTTCTTTTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACCGGATTCGTTTTCTCTGACTCTTCCCCAAAAGCTGCTACGCGTTGTGTAGtttgttttgt                           >  2:202530/1‑141 (MQ=255)
                                                                                                                   cAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCa                         >  1:712636/1‑141 (MQ=255)
                                                                                                                   cAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCa                         >  1:586484/1‑141 (MQ=255)
                                                                                                                         gtagtaACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAg                   >  2:589475/1‑141 (MQ=255)
                                                                                                                         gtagtaACACTTTCATTGCACCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAg                   <  1:654698/141‑1 (MQ=255)
                                                                                                                         gtagtaACACTTTCATTGCACCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAg                   <  1:238345/141‑1 (MQ=255)
                                                                                                                           agtaACACAATCATTGCGCCTGCTTGTCCGCCATAGTATGATGCCAATTCATGCATATGATCTTGTACCACGTGTTTCGCGTTATCGAACTCTACCCCAAAACCTGTTAACTTATTGGTAGTTTGTGTTGTCATTACAGTa                 >  2:329315/1‑141 (MQ=255)
                                                                                                                            gtaACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTaa                >  2:22345/1‑141 (MQ=255)
                                                                                                                               acacTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTa             >  2:416026/1‑141 (MQ=255)
                                                                                                                                caatttCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATATTGTGTAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTaa            <  1:194358/138‑1 (MQ=255)
                                                                                                                                  cTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGa          >  1:46886/1‑141 (MQ=255)
                                                                                                                                    ttCATTGCGCCTTCTTGTCCGTCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGAtt        <  1:325868/141‑1 (MQ=255)
                                                                                                                                    ttCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGAtt        <  1:456906/141‑1 (MQ=255)
                                                                                                                                         tGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTCACTAAGATTGCTGt   >  2:679515/1‑141 (MQ=255)
                                                                                                                                         tGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGATTGCTGt   >  1:207927/1‑141 (MQ=255)
                                                                                                                                         tGCCCCTTTTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAATAGCTGATAAGTGTTGTGTAGTTTGTGTTTTCATTACAGTAACTAAGATTGCTGt   <  1:55599/141‑1 (MQ=255)
                                                                                                                                          gcgcCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGATTGCTGTa  <  1:643860/141‑1 (MQ=255)
                                                                                                                                          gcgcCCTCTTGTCCGCCAAATTGGGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGATTGCTGTa  <  1:150942/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CTTTTGCTTTTTTATTACTTTGTAAAAATAAACATAAAATTAAGGCGATGATGCTAAATATCGTTGCAACTATAAATGCATCATTAAT‑ACCTTCAACCGTTGCTAGTTTATTTACAAATTGTAGTAACACTTTCATTGCGCCTTCTTGTCCGCCATATTGTGATGCTAATTCATGCATATGATCTTGTACAACAGGATTCGTTTTATCTAACTCTTCCCCAAAAGCTGATAAGTGTTGTGTAGTTTGTGTTGTCATTACAGTAACTAAGATTGCTGTA  >  USA300TCH1516_ALE/2469950‑2470227

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: ‑