Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,807,163 C→T 6.8% G487R (GGG→AGG)  asp2 ← Accessory Sec system protein Asp2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,807,1630CT 6.8% 200.0 / 3.3 73G487R (GGG→AGG) asp2Accessory Sec system protein Asp2
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base T (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

GTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGA  >  USA300TCH1516_ALE/2807023‑2807299
                                                                                                                                            |                                                                                                                                         
gTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCTTCTTTATTTGGccc                                                                                                                                           >  1:355164/1‑141 (MQ=255)
   tgtcAAATATAATTCCAAGGAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAgg                                                                                                                                        <  1:406778/141‑1 (MQ=255)
   tactaaaTAAGATTCCAAGTAACTTGAAATGTTTTGTATTTTAGCAAACTGAACTCTCCCAAACGTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAgg                                                                                                                                        <  1:567997/137‑1 (MQ=255)
    gtgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTCTGTCGCCCAGGa                                                                                                                                       >  1:813091/1‑141 (MQ=255)
    gtgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTAAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACAAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGa                                                                                                                                       >  2:620152/1‑141 (MQ=255)
     tgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCGACGaa                                                                                                                                      >  1:61517/1‑141 (MQ=255)
     tgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCAACATTACGTCGCCCaagca                                                                                                                                      >  2:118645/1‑137 (MQ=255)
      gAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCACCATTATGTCGCCCAGGAAt                                                                                                                                     <  1:545677/141‑1 (MQ=255)
      gAAATATAATCCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTAAATACTTTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAAt                                                                                                                                     >  1:607963/1‑141 (MQ=255)
          tataATTCCAAGTAACTTTAAATGTTTTGTTTTTCTGCATGTTGAACTCTCCCAAACTTGTCTTCCAAATTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTg                                                                                                                                 >  2:234033/1‑141 (MQ=255)
             aaTTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATATTGAACGCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAAGGGGAGGTGAATCATCATTATGTCGCCACGAAATACTGCGa                                                                                                                              >  2:457934/1‑141 (MQ=255)
              tttCTAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGAt                                                                                                                             <  2:247592/140‑1 (MQ=255)
                tCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCGTTATGTCGCCCAGGACTACTGCGATtc                                                                                                                           >  1:392619/1‑141 (MQ=255)
                  cAAGTACCTTTAATTGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAATATAATGTTGTAAAAATTAACAAACCAACTTGCTATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATtcat                                                                                                                         <  2:27677/141‑1 (MQ=255)
                   aaGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATtcatc                                                                                                                        <  2:443037/141‑1 (MQ=255)
                     gtgaccttAAATGTTTTCTATTTCAGCATGTTGAACTCTCCCACACTTGTCTTCCAAAAAAAAAAAATAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  2:744907/135‑1 (MQ=255)
                     gTAACTTTAATTGTTTTGTTTTTCAGCATGATGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  2:668072/141‑1 (MQ=255)
                     gTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCACAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      >  2:139636/1‑141 (MQ=255)
                     gTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  1:693890/141‑1 (MQ=255)
                          tttAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTcttg                                                                                                                 >  1:580358/1‑141 (MQ=255)
                              aaTGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAACACTGCGATTCATCACTcttgcttg                                                                                                             <  1:533739/141‑1 (MQ=255)
                                tgttttgtttTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGgt                                                                                                           >  1:545428/1‑141 (MQ=255)
                                   tttgttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGtgtg                                                                                                        <  2:179232/141‑1 (MQ=255)
                                      gttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTg                                                                                                     <  2:418316/141‑1 (MQ=255)
                                      gttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTg                                                                                                     <  1:624318/141‑1 (MQ=255)
                                       tttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGt                                                                                                    >  2:177457/1‑141 (MQ=255)
                                       tttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATCCATCACTCTTGCTTGGTGTGCTGt                                                                                                    >  2:320870/1‑141 (MQ=255)
                                        ttttCAGCATTTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAGAATTAACAAACCAACTTGCAATGGTAGGTGAATAATCATTATGTCGCCCAGATATACTTCGATTTATCACTTTTGCTTGGTTTGCTGTc                                                                                                   >  2:91826/1‑141 (MQ=255)
                                        ttttCAGCATGTTTAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTc                                                                                                   <  2:265290/141‑1 (MQ=255)
                                        ttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTACTGTc                                                                                                   <  2:731612/141‑1 (MQ=255)
                                        tttgCATCATGTTGATCTCCCCCAACCTCGTCTGACAACATAAGGTTGTACAAATTAACAAACCAACTTGCAATGGTATGGGACTCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTc                                                                                                   <  2:91629/141‑1 (MQ=255)
                                          ttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACTAACTTGCAATGGTTGGTGGATCATTATTATGTTGAGCAGGGATTGTGGGATTCATCACTTTTGCTTGGTTTGCTGTCat                                                                                                 >  2:281556/1‑140 (MQ=255)
                                           tCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAAt                                                                                                >  2:220216/1‑141 (MQ=255)
                                           tCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGCCCTGTCAAt                                                                                                >  2:571213/1‑141 (MQ=255)
                                                  gTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTa                                                                                         <  1:821332/141‑1 (MQ=255)
                                                      aaCTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAg                                                                                     <  1:95429/141‑1 (MQ=255)
                                                        ctctCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCCGGAATACTGCGAGTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTaaaaact                                                                                   >  1:784046/1‑136 (MQ=255)
                                                        ctctCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGct                                                                                  >  2:130079/1‑141 (MQ=255)
                                                         tctcCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAAAAAACCAACGTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGctc                                                                                 <  1:782104/141‑1 (MQ=255)
                                                                   ttGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCa                                                                       <  2:437225/141‑1 (MQ=255)
                                                                    tgttatCCAAAATAATGTGGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCTATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCAt                                                                      <  2:463598/136‑1 (MQ=255)
                                                                        ttCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTgg                                                                  >  1:357852/1‑141 (MQ=255)
                                                                        ttCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTgg                                                                  <  2:740293/141‑1 (MQ=255)
                                                                         tCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTAGGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAAGACAGGTAATAGCTCTTGAAATGCATGTGGa                                                                 >  1:379303/1‑141 (MQ=255)
                                                                            aaaaTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCa                                                              >  2:373376/1‑141 (MQ=255)
                                                                            aaaaTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCa                                                              <  1:837952/141‑1 (MQ=255)
                                                                                 aatGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATAATTATGCAGCTCACGAATACTGCGATTCATCACCCTTGCTTGGTGTGCTGTCAATACAGGTAATAGACCTTGCACTGCAGGTGGATCATaacc                                                         >  2:346339/1‑139 (MQ=255)
                                                                                  atGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAatca                                                        >  2:54068/1‑141 (MQ=255)
                                                                                  atGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATCTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCCTAatca                                                        <  2:526287/141‑1 (MQ=255)
                                                                                   tgttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAATTGCATGTGGATCATAatcgt                                                       >  2:638333/1‑139 (MQ=255)
                                                                                   tgttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAatcat                                                       <  1:599849/141‑1 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  1:380586/1‑141 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAAtcatca                                                     <  1:93331/141‑1 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  2:192151/1‑141 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGAAATAGCTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  2:65104/1‑141 (MQ=255)
                                                                                          aaaaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTg                                                <  2:399648/141‑1 (MQ=255)
                                                                                           aaaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTCAAATTACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGc                                               >  1:649372/1‑141 (MQ=255)
                                                                                            aaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATATAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCa                                              <  2:256373/141‑1 (MQ=255)
                                                                                             aTTAACAAACCAACATGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTCGTGTGCTGTCAAGACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCat                                             >  2:572695/1‑141 (MQ=255)
                                                                                                    aaCCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCt                                      >  1:572277/1‑141 (MQ=255)
                                                                                                           ttGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGGTCATAATCGACATGTTGCATATATGCTTTTGCaa                               >  1:189433/1‑141 (MQ=255)
                                                                                                            tGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCaaa                              <  1:177457/141‑1 (MQ=255)
                                                                                                            tGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGGCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCaaa                              <  1:320870/141‑1 (MQ=255)
                                                                                                               aaTGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGATTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAAc                           <  2:297948/141‑1 (MQ=255)
                                                                                                                atcgTAGGTGAATAATTATTCTGTCGCCCAGGACTATTGCGATTTCTCGCTTTTCCTTGGTGGTCGGGCAATAAAGGTCTCAGACTGTTTCTTGCATTTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACa                          <  2:804179/138‑1 (MQ=255)
                                                                                                                 tGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGCCACGCATGCGGATCATAATAATCATGTTGCATCTATGCTACGGCAAAAACAg                         >  1:798979/1‑141 (MQ=255)
                                                                                                                 tGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAg                         >  1:268742/1‑141 (MQ=255)
                                                                                                                     aGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGATtg                     >  1:558620/1‑141 (MQ=255)
                                                                                                                          aaTCATCATTATGTCGCCCAGGAACACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAATACAGATTGTGACa                >  1:611474/1‑141 (MQ=255)
                                                                                                                           atcatcATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACaa               <  1:14303/141‑1 (MQ=255)
                                                                                                                              atcatTATGTGGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGa            <  2:255320/141‑1 (MQ=255)
                                                                                                                                 atTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGAttt         >  2:194817/1‑141 (MQ=255)
                                                                                                                                        cGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGa  <  1:384780/141‑1 (MQ=255)
                                                                                                                                        cGCACAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGa  <  1:255370/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAGCTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGA  >  USA300TCH1516_ALE/2807023‑2807299

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: