Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,622,775 A→T 9.5% K43* (AAA→TAA)  marA → Multiple antibiotic resistance protein MarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,622,7750AT 9.5% 498.5 / 35.1 147K43* (AAA→TAA) marAMultiple antibiotic resistance protein MarA
Reads supporting (aligned to +/- strand):  ref base A (63/70);  new base T (9/5);  total (72/75)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 155 total aligned reads displayed.

TATTGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTACAC  >  USA300TCH1516_ALE/1622687‑1622862
                                                                                        |                                                                                       
tatTGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTaaa                                                                                       <  1:1929312/91‑1 (MQ=255)
tatTGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTaaa                                                                                       <  1:4117887/91‑1 (MQ=255)
 atTGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTtaaa                                                                                      <  1:2408028/91‑1 (MQ=255)
  ttGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTTAAAt                                                                                     <  1:4026878/91‑1 (MQ=255)
  ttGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTTAAAt                                                                                     <  1:793018/91‑1 (MQ=255)
   tGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTCAAAATg                                                                                    <  1:3226317/91‑1 (MQ=255)
      aGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTTAAATGAtt                                                                                 <  1:13070/91‑1 (MQ=255)
      aGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGAtt                                                                                 >  1:2529390/1‑91 (MQ=255)
           gTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCgg                                                                            >  1:2059061/1‑91 (MQ=255)
           gTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCgg                                                                            >  1:3293049/1‑91 (MQ=255)
           gTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCgg                                                                            >  1:1519305/1‑91 (MQ=255)
           gTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCgg                                                                            >  1:837352/1‑91 (MQ=255)
            tttATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATCACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGc                                                                           >  1:3483425/1‑91 (MQ=255)
            tttATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATCACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGc                                                                           >  1:3483411/1‑91 (MQ=255)
            tttATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATAACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGc                                                                           >  1:3483313/1‑91 (MQ=255)
             ttattaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCt                                                                          <  1:451565/91‑1 (MQ=255)
             ttattaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCt                                                                          >  1:1895689/1‑91 (MQ=255)
             ttattaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCt                                                                          >  1:1895989/1‑91 (MQ=255)
              tattaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCtt                                                                         >  1:3909419/1‑91 (MQ=255)
               attaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCGTTCGCCATCCCAGCTTGATCAATCATTTAAAATGATTGTCGGCTTa                                                                        <  1:3319241/91‑1 (MQ=255)
                ttaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTAt                                                                       >  1:489308/1‑91 (MQ=255)
                ttaGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTAt                                                                       >  1:3646932/1‑91 (MQ=255)
                 taGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTAtc                                                                      >  1:2141003/1‑91 (MQ=255)
                 taGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTAtc                                                                      <  1:2006406/91‑1 (MQ=255)
                     gCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCa                                                                  >  1:970827/1‑91 (MQ=255)
                     gCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCa                                                                  <  1:3365530/91‑1 (MQ=255)
                     gCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCa                                                                  >  1:3359353/1‑91 (MQ=255)
                     gCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCa                                                                  >  1:970814/1‑91 (MQ=255)
                       cTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGa                                                                <  1:2086199/91‑1 (MQ=255)
                        tttCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGaa                                                               <  1:3945825/91‑1 (MQ=255)
                        tttCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGaa                                                               >  1:567710/1‑91 (MQ=255)
                         ttCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAg                                                              <  1:3429961/91‑1 (MQ=255)
                           cAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCt                                                            <  1:3700064/91‑1 (MQ=255)
                           cAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCt                                                            >  1:114416/1‑91 (MQ=255)
                            aaTTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCtt                                                           <  1:348246/91‑1 (MQ=255)
                             aTTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTa                                                          <  1:2892052/91‑1 (MQ=255)
                              tttGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTAt                                                         >  1:679526/1‑91 (MQ=255)
                                tGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGc                                                       >  1:3558518/1‑91 (MQ=255)
                                 gCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCa                                                      >  1:4116917/1‑91 (MQ=255)
                                 gCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCa                                                      <  1:904375/91‑1 (MQ=255)
                                  cAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAc                                                     >  1:2901257/1‑91 (MQ=255)
                                     gAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcg                                                  <  1:3032722/91‑1 (MQ=255)
                                     gAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcg                                                  <  1:2520170/91‑1 (MQ=255)
                                      aaTTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgc                                                 >  1:685127/1‑91 (MQ=255)
                                      aaTTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgc                                                 <  1:227550/91‑1 (MQ=255)
                                      aaTTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgc                                                 >  1:2640529/1‑91 (MQ=255)
                                         tAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgcgcg                                              <  1:2872109/91‑1 (MQ=255)
                                         tAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgcgcg                                              <  1:2876672/91‑1 (MQ=255)
                                         tAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCAcgcgcgcg                                              <  1:3431529/91‑1 (MQ=255)
                                          aaGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGTGCGCGt                                             >  1:1605135/1‑91 (MQ=255)
                                            gTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTaa                                           >  1:602715/1‑91 (MQ=255)
                                             tGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTaaa                                          <  1:1286594/91‑1 (MQ=255)
                                              gATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTaaaa                                         <  1:2517001/91‑1 (MQ=255)
                                              gATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTaaaa                                         >  1:3902385/1‑91 (MQ=255)
                                              gATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTaaaa                                         >  1:2978992/1‑91 (MQ=255)
                                               aTTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAAt                                        <  1:83268/91‑1 (MQ=255)
                                                ttACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATg                                       <  1:969419/91‑1 (MQ=255)
                                                 tACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGa                                      <  1:2370868/91‑1 (MQ=255)
                                                 tACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGa                                      <  1:2674185/91‑1 (MQ=255)
                                                  ctgTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGac                                     <  1:1166363/89‑1 (MQ=255)
                                                  aCGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGac                                     <  1:1263364/91‑1 (MQ=255)
                                                    gTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGacac                                   <  1:1630402/91‑1 (MQ=255)
                                                      tGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTc                                 >  1:662702/1‑91 (MQ=255)
                                                      tGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTc                                 >  1:662733/1‑91 (MQ=255)
                                                         tCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCt                              >  1:3440006/1‑91 (MQ=255)
                                                         tCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCt                              <  1:3645278/91‑1 (MQ=255)
                                                          cTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTg                             <  1:1014554/91‑1 (MQ=255)
                                                          cTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTg                             <  1:3689978/91‑1 (MQ=255)
                                                            ttCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCa                           >  1:1302874/1‑91 (MQ=255)
                                                              cGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCaaa                         >  1:579978/1‑91 (MQ=255)
                                                              cGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCaaa                         <  1:1601061/91‑1 (MQ=255)
                                                               gCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAat                        <  1:408952/91‑1 (MQ=255)
                                                                  aTACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAatgat                     >  1:2784471/1‑91 (MQ=255)
                                                                  aTACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAatgat                     <  1:4131279/91‑1 (MQ=255)
                                                                   tACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAatgatg                    <  1:1485654/91‑1 (MQ=255)
                                                                     ccATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGAtgtg                  >  1:2696025/1‑91 (MQ=255)
                                                                     ccATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGAtgtg                  >  1:2467534/1‑91 (MQ=255)
                                                                      cATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGa                 <  1:4070237/91‑1 (MQ=255)
                                                                        tCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGAtt               <  1:670537/91‑1 (MQ=255)
                                                                         cTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTa              >  1:2455325/1‑91 (MQ=255)
                                                                         cTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTa              <  1:487896/91‑1 (MQ=255)
                                                                          ttGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTaa             >  1:1806058/1‑91 (MQ=255)
                                                                           tGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAAt            <  1:169670/91‑1 (MQ=255)
                                                                           tGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAAt            <  1:2647495/91‑1 (MQ=255)
                                                                             atcaatcaTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATgg          >  1:3864009/1‑91 (MQ=255)
                                                                                aatcaTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGc       >  1:3795732/1‑91 (MQ=255)
                                                                                 atcaTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCt      >  1:1029528/1‑91 (MQ=255)
                                                                                   caTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTac    <  1:3561843/91‑1 (MQ=255)
                                                                                   caTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTac    <  1:2121727/91‑1 (MQ=255)
                                                                                   caTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTac    >  1:2577814/1‑91 (MQ=255)
                                                                                   caTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTac    >  1:1076847/1‑91 (MQ=255)
                                                                                    aTTTTAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   >  1:4015765/1‑91 (MQ=255)
                                                                                    aTTTTAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   >  1:1674642/1‑91 (MQ=255)
                                                                                    aTTTTAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   >  1:2192142/1‑91 (MQ=255)
                                                                                    aTTTTAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   >  1:2192158/1‑91 (MQ=255)
                                                                                    aTTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   <  1:2975237/91‑1 (MQ=255)
                                                                                    aTTTAAAATGATAGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTaca   >  1:3198024/1‑91 (MQ=255)
                                                                                     tttAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTacac  <  1:1447920/91‑1 (MQ=255)
                                                                                     tttAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTacac  >  1:1867190/1‑91 (MQ=255)
                                                                                     tttAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTacac  <  1:1180986/91‑1 (MQ=255)
                                                                                        |                                                                                       
TATTGAAGATCGTTTATTAGAGCCTTTCAATTTGCAAGAATTAAGTGATTACGTTGGTCTTTCGCCATACCATCTTGATCAATCATTTAAAATGATTGTCGGCTTATCTCCAGAAGCTTATGCACGCGCGCGTAAAATGACACTCGCTGCAAATGATGTGATTAATGGTGCTACAC  >  USA300TCH1516_ALE/1622687‑1622862

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