Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,744,544 T→C 58.1% Y43C (TAC→TGC)  apt ← Adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,744,5440TC58.1% 38.0 / 66.3 69Y43C (TAC→TGC) aptAdenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (22/18);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCT  >  USA300TCH1516_ALE/1744406‑1744682
                                                                                                                                          |                                                                                                                                          
cGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTAt                                                                                                                                          >  1:285698/1‑141 (MQ=255)
   aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAGTCCCCATCGAATAAGCTACAGGACATCCAATGATAAAGCCACGCGCTTCAGCTCCTACAACGATCTCTAATGCCCTGTCTTTTGCGTATTCt                                                                                                                                       <  1:814707/141‑1 (MQ=255)
   aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCt                                                                                                                                       <  2:29181/141‑1 (MQ=255)
   aTGACTGCACGAGGTAATTTCCCTTCGTTTCTAACAGTTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                                                                                                                       <  2:264397/141‑1 (MQ=255)
     gACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCAGGGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc                                                                                                                                     >  1:767257/1‑141 (MQ=255)
         tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAAtt                                                                                                                                 <  2:46622/141‑1 (MQ=255)
         tCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGACTTTTGCGCATTCTACAAtt                                                                                                                                 >  1:556177/1‑141 (MQ=255)
                    tttCCCCACATTTCCAACAAGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCCACAACGATATCAACAACTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGc                                                                                                                      <  1:397118/141‑1 (MQ=255)
                            cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCat                                                                                                              >  1:892562/1‑141 (MQ=255)
                            cttTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCat                                                                                                              >  1:38271/1‑141 (MQ=255)
                             ttttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCata                                                                                                             >  2:1013264/1‑141 (MQ=255)
                               ttCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATg                                                                                                           >  2:856364/1‑141 (MQ=255)
                                 cTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCt                                                                                                         >  1:710823/1‑141 (MQ=255)
                                  tAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTGTTTTATAATTTTTTTTGTTGGATGGCCATATGCtt                                                                                                        >  1:99996/1‑141 (MQ=255)
                                  tAACAGGTGCAAAGCCAATCCCCAGTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTATTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCtt                                                                                                        >  2:835417/1‑141 (MQ=255)
                                    aCAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCa                                                                                                      <  1:286288/141‑1 (MQ=255)
                                     cAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCAc                                                                                                     <  1:490854/141‑1 (MQ=255)
                                     cAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAGGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCCTAGCCATATGCTTCAc                                                                                                     >  2:930818/1‑141 (MQ=255)
                                      aGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCAcc                                                                                                    <  1:425393/141‑1 (MQ=255)
                                      aGGTGCAAAGCCAATCCACATTGAAGAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCAcc                                                                                                    <  2:517821/141‑1 (MQ=255)
                                        gTGCAAAGCCAATCCCCATTGAATAAGGTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATGCACATCTCTGTCTTTTGCTTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCAt                                                                                                  >  2:827840/1‑141 (MQ=255)
                                             aacccAATCCCCATTGAATAAGCAACAGGACAGCCAATGATAAAGCCACGCGCCTCAGGTCCTACAACGATATCAACATCTCTTTCTTTTGCCCATCCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                                                                                             <  2:156345/138‑1 (MQ=255)
                                             aaGCCCATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                                                                                             <  1:684485/141‑1 (MQ=255)
                                             aaGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACACCTCTGTCTTTTGCGCATTCTACAATTTTAGCTGTTGCATAGCCATATGCTTCACCATTATcc                                                                                             <  1:296734/141‑1 (MQ=255)
                                               cccAATCCCCGTGGAATAACCTACCGGACATCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCAt                                                                                           <  2:40698/140‑1 (MQ=255)
                                                ccAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATa                                                                                          >  2:949422/1‑141 (MQ=255)
                                                       ccATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTgtag                                                                                   <  2:851990/141‑1 (MQ=255)
                                                            gAATAAGCTCCAGGCCAGCCAATGATAAAGCCACGCGCTTCAGGGCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAata                                                                              <  2:759984/141‑1 (MQ=255)
                                                             aaTAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAatat                                                                             >  2:788095/1‑141 (MQ=255)
                                                               tAAGCTACAGGACAGCCAATGATAAAGCCACGCGCGTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTCTCTGTTGCATAGCCATATGCTTCCCCAGTATCCATACTTgtagtcagctcc                                                                           >  2:150584/1‑134 (MQ=255)
                                                                aaGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCCCCATTATCCATATTTGTAGTAATATCCt                                                                          >  1:1023901/1‑141 (MQ=255)
                                                                  gCTACAGGACAGCCAATGATCAAGCCACGCGCTCCAGGTCCTACAACGATATCAACCTCTCTGTCTTTTGCGTATTCGACCATTTTTTCGGTTGCATAGCCTTATACTTCGCCATTTTCCATAATTGTAGTAATCACCTTg                                                                        >  2:756519/1‑141 (MQ=255)
                                                                          aCAGCCAATGATAAAGCCAGGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTCTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAac                                                                >  2:617008/1‑141 (MQ=255)
                                                                           cAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAcca                                                               >  2:828737/1‑141 (MQ=255)
                                                                           cAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAaca                                                               >  2:579505/1‑141 (MQ=255)
                                                                               cAATGATAAAGCCAGGCGCTTCAGGTCCGACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTg                                                           >  1:721464/1‑141 (MQ=255)
                                                                               cAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTg                                                           >  1:772550/1‑141 (MQ=255)
                                                                                 aTGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTTTAGTAATATCCTTGAAACTAACACCTGGt                                                         >  2:416970/1‑141 (MQ=255)
                                                                                  tGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGtt                                                        >  2:131632/1‑141 (MQ=255)
                                                                                  tGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCAATAATTTTAGTATTATCCTTGAGACTAACACCTGGtt                                                        >  1:952141/1‑141 (MQ=255)
                                                                                       aaGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGc                                                   >  2:918303/1‑141 (MQ=255)
                                                                                          ccACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCaa                                                <  2:1000887/141‑1 (MQ=255)
                                                                                          ccACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCaa                                                <  1:619287/141‑1 (MQ=255)
                                                                                               gccTTCAGGTCCTACAACGATCTCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                           <  2:122716/139‑1 (MQ=255)
                                                                                               cgcTTCAGGTCCTACCACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                           <  1:1110513/141‑1 (MQ=255)
                                                                                               cgcTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                           <  1:763323/141‑1 (MQ=255)
                                                                                               cgcTTCAGGCTCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCCGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTACTTTCCTTGAAACTAACACCTGGTTTCGGCCAATCTTg                                           <  1:41264/141‑1 (MQ=255)
                                                                                                 cTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGaa                                         >  1:16606/1‑141 (MQ=255)
                                                                                                          cTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTATCACCTTGGTTCGGCCAATCTTGAACTTCTGATa                                <  2:285698/141‑1 (MQ=255)
                                                                                                           tACAACGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATAc                               >  2:88089/1‑141 (MQ=255)
                                                                                                           tACAACGATATCAACAGCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATag                               >  1:1110175/1‑140 (MQ=255)
                                                                                                              aaCGATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTa                            >  2:1014410/1‑141 (MQ=255)
                                                                                                                 gATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTg                         >  1:364787/1‑141 (MQ=255)
                                                                                                                 gATATCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTg                         >  1:557718/1‑141 (MQ=255)
                                                                                                                    atCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCtt                      <  2:899848/141‑1 (MQ=255)
                                                                                                                    atCAACATCTCTGTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCGGGTTTCGGCCAATCTTGACCTGCTGATACGTATTGCtt                      <  1:617008/141‑1 (MQ=255)
                                                                                                                     tCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCACTCTTGAACTTCTGATTCTTGTTGCttt                     >  2:594485/1‑141 (MQ=255)
                                                                                                                     tCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCttt                     >  1:342088/1‑141 (MQ=255)
                                                                                                                     tCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCttt                     >  2:44597/1‑141 (MQ=255)
                                                                                                                     tCAACAGCTCTGGCTGTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAAATCACAGATTGGGTTAGCCATTCTTGAACTTCTGATACGTATTGCttt                     >  2:228797/1‑141 (MQ=255)
                                                                                                                             cctGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCAt             <  1:599266/140‑1 (MQ=255)
                                                                                                                                ttCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATCTGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGACCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTaa          <  2:262485/140‑1 (MQ=255)
                                                                                                                                gTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCCTATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGACCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTaa          <  1:1054318/141‑1 (MQ=255)
                                                                                                                                gTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTaa          <  1:854842/141‑1 (MQ=255)
                                                                                                                                gTCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTaa          >  1:899368/1‑141 (MQ=255)
                                                                                                                                 tCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAat         <  1:1046529/141‑1 (MQ=255)
                                                                                                                                 tCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAat         <  1:616257/141‑1 (MQ=255)
                                                                                                                                 tCTTTTGCGCATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTATCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAATTCCATTAat         <  2:710823/141‑1 (MQ=255)
                                                                                                                                        cGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATACATACTTGTAGTAATATCCTTGAAACTACAATCTGGTTTATCTTAAGCTGTAATTTCTGATACGTATTGCTTTAAATCCATTAATATTtcct  >  2:1062217/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                          
CGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAACTAACACCTGGTTTCGGCCAATCTTGAACTTCTGATACGTATTGCTTTAAATCCATTAATATTTCCT  >  USA300TCH1516_ALE/1744406‑1744682

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: