Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,373,506 A→T 15.7% N90K (AAT→AAA pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,5060AT15.7% 143.1 / 24.6 65N90K (AAT→AAApbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base A (36/18);  new base T (3/7);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-02

TAAGGAATTGCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAAC  >  USA300TCH1516_ALE/2373366‑2373646
                                                                                                                                             |                                                                                                                                             
tAAGGAATTGCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCa                                                                                                                                               <  1:531348/141‑1 (MQ=255)
    gAATTGCATTAATAATAACTTCCCTAAGCCCGGGGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTa                                                                                                                                           >  1:455479/1‑141 (MQ=255)
           aTTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTACTCCCAt                                                                                                                                   >  2:739871/1‑141 (MQ=255)
               ataataaCTTCCCTAAACCCTATGACGGTTAATATCGCAAAGAATATACCTGAAAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCt                                                                                                                               <  1:511443/141‑1 (MQ=255)
                 aataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTgg                                                                                                                             <  1:816272/141‑1 (MQ=255)
                 aataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTgg                                                                                                                             <  2:10533/141‑1 (MQ=255)
                  ataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAACACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGt                                                                                                                            >  1:835051/1‑141 (MQ=255)
                   taaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTg                                                                                                                           >  1:98723/1‑141 (MQ=255)
                   taaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTg                                                                                                                           <  1:461210/141‑1 (MQ=255)
                    aaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAAAATGCATTTAATCCAATACCTGGTGc                                                                                                                          >  2:277026/1‑141 (MQ=255)
                     aCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCt                                                                                                                         >  2:926441/1‑141 (MQ=255)
                      cTTCCCTAAACCCTGTGACGGTTAATATCGCATAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCACGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTa                                                                                                                        <  1:630810/141‑1 (MQ=255)
                            tAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGCCTCCCTCATGTGTACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCaa                                                                                                                  >  1:448083/1‑141 (MQ=255)
                            tAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGTTTAACACAACAGTAAATGCAAAGAATGCTTTTAAGCCCATACCTGGTGCTAACGCaa                                                                                                                  <  2:269877/141‑1 (MQ=255)
                             aaaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAAt                                                                                                                 <  1:104679/141‑1 (MQ=255)
                                cccTGCGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTAACGCCACAACCTGGTGCTAACGCAATTgg                                                                                                              <  1:456888/141‑1 (MQ=255)
                                 ccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACCTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGa                                                                                                             <  2:7943/141‑1 (MQ=255)
                                    gtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATAGGatat                                                                                                          >  2:937816/1‑141 (MQ=255)
                                    gtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACATGGTGCTAACGCAATTGTatat                                                                                                          >  1:186414/1‑141 (MQ=255)
                                    gtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACCCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTACTAACGCAATTGGatat                                                                                                          >  1:490161/1‑141 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAGTACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGAAATTTa                                                                                                       >  2:447048/1‑141 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACCACAGTAAATGCACAGAATGCATTTAATAGAATAACTGGGGCTAACGCAATTGGATAtttc                                                                                                       >  2:150402/1‑140 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGGTATTTa                                                                                                       >  2:161501/1‑141 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATAttca                                                                                                       >  1:331576/1‑139 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTa                                                                                                       >  1:865261/1‑141 (MQ=255)
                                       aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGGGCTAACGCAATTGGATAtttg                                                                                                       >  1:419867/1‑140 (MQ=255)
                                        cGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAg                                                                                                      >  1:920449/1‑141 (MQ=255)
                                              atatCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCGATTGGATATTTAGCTATTa                                                                                                <  1:562849/141‑1 (MQ=255)
                                                 tCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGcc                                                                                             <  2:98723/141‑1 (MQ=255)
                                                   gCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGcccc                                                                                           <  1:716926/141‑1 (MQ=255)
                                                   ccAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGcccc                                                                                           <  1:23736/140‑1 (MQ=255)
                                                    cAAAGAATCTACCTGAGCATAAAACGCCTGCTAAACCAACTTGCCACGCAATACCCATGGTTAACACAACAGTCAACGCCAAGAATGCACTTACTCCCATACCTGGTGCTAACGCACTTGCCCACTTAGCTATCCGccccc                                                                                          >  2:864575/1‑140 (MQ=255)
                                                    cAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCa                                                                                          <  2:455479/141‑1 (MQ=255)
                                                       aGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTGGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGa                                                                                       <  1:779826/141‑1 (MQ=255)
                                                           tataCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAg                                                                                   <  2:187452/141‑1 (MQ=255)
                                                              aCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCAATGGTTAACACAACAGTAAATGCAAAGAATGCAGTTAATCCAATACCTGGTGCTAACGCAAATGGATAGTTAGCTATTAGCCCCATGAATAGCgg                                                                                >  2:844930/1‑140 (MQ=255)
                                                                   agaATAAACCGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTa                                                                           <  2:72705/141‑1 (MQ=255)
                                                                     aaTAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACa                                                                         >  1:577079/1‑141 (MQ=255)
                                                                                 gTTAAACCAACTTGCCAAGGAATACCCATGGTTCACACCACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAAGTCCATCGGCGCGTGATGTAGGAGCCCCATGAATAGCGAGCATACAAATGCTGCTaag                                                             >  1:306826/1‑140 (MQ=255)
                                                                                    aaaCCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGATGCTAATGCa                                                          >  1:633476/1‑141 (MQ=255)
                                                                                        cAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGAGATTTAGCTATTAGCCCCATGAATAGCGAGCGTACAAATGCCGCTAATGCGGTcc                                                      >  2:330861/1‑140 (MQ=255)
                                                                                         aaCTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGc                                                     <  2:362337/141‑1 (MQ=255)
                                                                                           cctGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTa                                                   <  1:855822/139‑1 (MQ=255)
                                                                                             ggCCAAGGAATACCCAGGGTTAACACAACGGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACa                                                 <  1:864575/140‑1 (MQ=255)
                                                                                              gCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaa                                                >  1:492129/1‑141 (MQ=255)
                                                                                                 aaGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaaa                                             >  2:819116/1‑141 (MQ=255)
                                                                                                  atgAATACCCATGTTTAACACAACATTAAATGCAAAGAATGCATTTAACCCCATACCTGGTGCTAACGCAATTGGATATTTAGCGATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAAt                                            <  2:539654/139‑1 (MQ=255)
                                                                                                  aGGAATACCCATGGCTAACACGACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATGTTTAGCTATTAGACCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAAt                                            >  2:646544/1‑141 (MQ=255)
                                                                                                        aCCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCAcc                                      >  1:752701/1‑141 (MQ=255)
                                                                                                         cccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTGTTTATTATCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCt                                     >  1:686564/1‑141 (MQ=255)
                                                                                                              tGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGt                                 <  2:308482/141‑1 (MQ=255)
                                                                                                              tGGTTAAAACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCGAACTTAATTGGATATTTAGCTATTCGCCCGATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGt                                 <  2:705206/141‑1 (MQ=255)
                                                                                                                    agacaacaGTCAATACAAGGAATGCATTTAATGTCATAACTAATGGAACTGCATTTATTGAGTAAGTTGTTGGCGCCAGGAATAGCGAGACTACAAATGCTGCTAATGCAGTGGCTAAAAAAGGGGCCACGTGGTCAAttt                           >  2:429974/3‑141 (MQ=255)
                                                                                                                    acacAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                           >  1:578818/1‑141 (MQ=255)
                                                                                                                    acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCCGTCGCTACAAAAATGGCACCTTGGTGCCttt                           >  1:843/1‑141 (MQ=255)
                                                                                                                    acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                           >  1:144763/1‑141 (MQ=255)
                                                                                                                    acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                           >  1:493949/1‑141 (MQ=255)
                                                                                                                     tacaacaGTAAATGCAAAGTATGCATTTAATCCCATTCCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAtttt                          <  1:157232/140‑1 (MQ=255)
                                                                                                                     cacaACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAtttt                          <  2:114710/141‑1 (MQ=255)
                                                                                                                              aaaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttc                 <  2:733594/141‑1 (MQ=255)
                                                                                                                               aaTGGAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACCAAAATGGCACCATGGTCAATGTTCGAAtcttct                >  2:571315/1‑141 (MQ=255)
                                                                                                                               aaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTGCATTTTCATAtcttct                >  2:412281/1‑141 (MQ=255)
                                                                                                                               aaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct                >  2:828842/1‑141 (MQ=255)
                                                                                                                               aaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAccttct                >  2:224222/1‑141 (MQ=255)
                                                                                                                                      aGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc         >  1:146781/1‑141 (MQ=255)
                                                                                                                                           gCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTa    >  2:876030/1‑141 (MQ=255)
                                                                                                                                             aTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACACAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTacc  >  2:259595/1‑139 (MQ=255)
                                                                                                                                             aTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCCCATTGGTCCCTTTTCCTATCTTCTGAAAAACCTTTAac  >  1:377954/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TAAGGAATTGCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAAC  >  USA300TCH1516_ALE/2373366‑2373646

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: