Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,373,525 G→T 50.0% A84E (GCA→GAA)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,5250GT50.0% 5.9 / 86.7 67A84E (GCA→GAA) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base G (17/14);  new base T (18/14);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

AAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAAC  >  USA300TCH1516_ALE/2373395‑2373664
                                                                                                                                  |                                                                                                                                           
aaaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAACCCAACTTGCCAAGGAATACCCATGGTTCACACAACAGTAAATGCAAAGAATACACTTAATCCCATACCTGGTTCTAACGCAAt                                                                                                                                   >  2:569626/1‑141 (MQ=255)
aaaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATGCCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTCGTTCTAACGCAAt                                                                                                                                   >  1:476207/1‑141 (MQ=255)
   cccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTgg                                                                                                                                >  2:121260/1‑141 (MQ=255)
   cccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTgg                                                                                                                                >  1:2541/1‑141 (MQ=255)
     cTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGat                                                                                                                              <  1:516704/141‑1 (MQ=255)
      tgtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGata                                                                                                                             <  1:246478/141‑1 (MQ=255)
       gtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGatat                                                                                                                            >  1:316536/1‑141 (MQ=255)
       gtgACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTCGctat                                                                                                                            >  2:232991/1‑141 (MQ=255)
             gTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCt                                                                                                                      <  1:318463/141‑1 (MQ=255)
             gTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCt                                                                                                                      <  1:39223/141‑1 (MQ=255)
             gTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCt                                                                                                                      <  2:764171/141‑1 (MQ=255)
                   atCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGc                                                                                                                <  2:410479/141‑1 (MQ=255)
                      gCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGcccc                                                                                                             >  1:480819/1‑141 (MQ=255)
                                taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCg                                                                                                   >  1:690300/1‑141 (MQ=255)
                                    tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTATCCCCATGAATCGCGAGcc                                                                                               >  2:222681/1‑141 (MQ=255)
                                    tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGcc                                                                                               <  2:346940/141‑1 (MQ=255)
                                          tAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACaaa                                                                                         <  1:384901/141‑1 (MQ=255)
                                          tAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACaaa                                                                                         >  1:714634/1‑141 (MQ=255)
                                           aaaaCGACTGTTAAACCAACTTGCCAAGGAATGCCTATGTGTAACACAACAGTAATTTCTAAGAATGCATTTAATCCTATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAAt                                                                                        <  2:316536/141‑1 (MQ=255)
                                           aaaaCGACTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAAt                                                                                        <  2:598888/141‑1 (MQ=255)
                                              aCGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAAtgct                                                                                     >  1:99598/1‑141 (MQ=255)
                                                      tAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGc                                                                             >  1:121870/1‑141 (MQ=255)
                                                      tAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCACCATGAATAGTGATTCTACAAATGCTGCTAATGc                                                                             >  2:672838/1‑141 (MQ=255)
                                                      tAAACCAACTTGCCAAAGAATACCCATGGTTAACACAAGAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGc                                                                             <  2:361477/141‑1 (MQ=255)
                                                             aCTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCCTGTATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCt                                                                      >  2:163746/1‑141 (MQ=255)
                                                              cTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTACTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTa                                                                     >  1:580424/1‑141 (MQ=255)
                                                                tGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACa                                                                   <  1:68338/141‑1 (MQ=255)
                                                                tGCCAACGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACa                                                                   <  2:46125/141‑1 (MQ=255)
                                                                  ccAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATTCTGCTAATGCAGTCGCTACaaa                                                                 >  2:406288/1‑141 (MQ=255)
                                                                  ccAAGAAATACCCATGGTTAAAACAACAGTAAATGCAAAGAATGCATTTAATTCAATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCAATGAATAGCGACCCTAAAAATGATACTAATACGGTCGCTAaaaa                                                                 >  2:108541/1‑141 (MQ=255)
                                                                   cAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTCCTAACGCAATTGGATATTTAGCTATTATCCCCATACATAGCGAGCCTACAAATCCTGCTAATGCAGTCGCTACaaaa                                                                >  2:237810/1‑141 (MQ=255)
                                                                   cAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATGTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaa                                                                <  1:33263/141‑1 (MQ=255)
                                                                        aaTACTCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGc                                                           <  2:702787/141‑1 (MQ=255)
                                                                            cccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCt                                                       >  2:449376/1‑141 (MQ=255)
                                                                            cccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTA‑CCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCtt                                                      >  1:237281/1‑141 (MQ=255)
                                                                               aTGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTgg                                                    <  1:272015/141‑1 (MQ=255)
                                                                                tGGTTAACACAACAGTAAATTCAAAGAATGCATTTAATCCCATACCTGGTACTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGt                                                   <  2:459464/141‑1 (MQ=255)
                                                                                      acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                                             >  2:102033/1‑141 (MQ=255)
                                                                                      acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTAGCTACAAAAATGGACCCTTGGTCCAttt                                             >  1:344483/1‑141 (MQ=255)
                                                                                      acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCACTCGCTACAACAATGGCACCTTGGTCCAttt                                             >  2:303603/1‑141 (MQ=255)
                                                                                      acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                                             >  2:743361/1‑141 (MQ=255)
                                                                                         caacaGTAAATGCAAAGAATGCATTTAATCCCATACCTCGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAGATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCa                                          <  2:738737/141‑1 (MQ=255)
                                                                                           acaGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGCCACCTTGGTCCATTTTCatt                                        >  2:314903/1‑140 (MQ=255)
                                                                                            caGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCatat                                       >  1:602101/1‑141 (MQ=255)
                                                                                             tgTAAATGCAAATAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtc                                      <  1:135024/140‑1 (MQ=255)
                                                                                              gTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTA‑CCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtctt                                    <  1:681387/141‑1 (MQ=255)
                                                                                               tAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAACATGGCCCCTTGGTCCCTCTTCATAtctt                                    >  2:132346/1‑141 (MQ=255)
                                                                                                  aTGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTg                                 <  1:605724/141‑1 (MQ=255)
                                                                                                  aTGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTGGGTCCATTTTCATGTCTTCTg                                 >  1:666087/1‑141 (MQ=255)
                                                                                                   tGCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGa                                >  2:652092/1‑141 (MQ=255)
                                                                                                    gCAAAGCATGCATTTAAGCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGAt                               <  1:649989/141‑1 (MQ=255)
                                                                                                    gCAAAGAATGCATTTAATCCCATACCTGGTTCTAACGCCATTGGATATTTAGCTCTTAGCCACCTAAATAGAGAGCCTTTATATGCTGCTATTGCAGTCGCTACTCAAATGGGTTTTTGCTTCTTTTTCATATCTTCTGAt                               <  2:681952/141‑1 (MQ=255)
                                                                                                        ggAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc                           <  2:700393/140‑1 (MQ=255)
                                                                                                          aaTGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCt                         >  2:754263/1‑141 (MQ=255)
                                                                                                          aaTGCATTTAATCCCATACCTGGTACTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCt                         <  2:73541/141‑1 (MQ=255)
                                                                                                           aTGCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATTCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCtt                        <  1:406288/141‑1 (MQ=255)
                                                                                                           aTGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCtt                        <  1:112017/141‑1 (MQ=255)
                                                                                                            tGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATACTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCttt                       >  2:238159/1‑141 (MQ=255)
                                                                                                             gCATTTAATCCCATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTa                      >  1:351490/1‑141 (MQ=255)
                                                                                                                  tAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGACACCTTGGTCCATTTTCATATCTTCTGATACTCCTTTAACAcc                 >  2:508632/1‑141 (MQ=255)
                                                                                                                        cATACCTGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTaa           >  2:585672/1‑141 (MQ=255)
                                                                                                                         aTACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTaaa          >  2:628833/1‑141 (MQ=255)
                                                                                                                          tACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCCGCTAAAc         >  2:582040/1‑141 (MQ=255)
                                                                                                                              tGGTTCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaa     <  2:74245/141‑1 (MQ=255)
                                                                                                                                gTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaaaa   >  1:110854/1‑141 (MQ=255)
                                                                                                                                 tGCTAACGCAATTTGATATTTAGCTATTAACCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAAc  <  2:655750/141‑1 (MQ=255)
                                                                                                                                 tGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACTCTTGCTTCTATTGCAGTCGCTAGATTGATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAAc  <  2:167584/141‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                           
AAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAAC  >  USA300TCH1516_ALE/2373395‑2373664

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