Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,373,667 G→A 19.3% Q37* (CAA→TAA)  pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,6670GA19.3% 168.1 / 30.9 86Q37* (CAA→TAA) pbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base G (29/40);  new base A (9/8);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

AACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAA  >  USA300TCH1516_ALE/2373528‑2373807
                                                                                                                                           |                                                                                                                                             
aaCGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACCCCTGCTAAACTTAAAACTTGc                                                                                                                                              <  1:702730/141‑1 (MQ=255)
  gtcAAGTGGATATTTAGCGATTAGCCCCATGAATAGCGAGCCTACAGATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCgg                                                                                                                                            <  2:218961/139‑1 (MQ=255)
  cGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAATATAGCACCTTGGTCCATTTTCATTTTTTTTGATACGCCTTTACCACCTGCTAAACTTAAAACTTACgg                                                                                                                                            <  2:130658/141‑1 (MQ=255)
      aTTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTa                                                                                                                                        >  1:539767/1‑141 (MQ=255)
        tGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAAc                                                                                                                                      <  2:534902/141‑1 (MQ=255)
        tGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAAc                                                                                                                                      <  1:583932/141‑1 (MQ=255)
        tGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAAc                                                                                                                                      <  2:714634/141‑1 (MQ=255)
         ggATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAACTGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTATCACCTCCTAAACTTAAAACTTACGGGTTCACt                                                                                                                                     >  2:592262/1‑141 (MQ=255)
         ggATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACt                                                                                                                                     <  1:713847/141‑1 (MQ=255)
            tatTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCt                                                                                                                                  <  2:650861/141‑1 (MQ=255)
              tttAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTaa                                                                                                                                <  1:107532/141‑1 (MQ=255)
              tttAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTaa                                                                                                                                <  2:697736/141‑1 (MQ=255)
                tAGCTATCAGTCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTaaaa                                                                                                                              <  2:76808/141‑1 (MQ=255)
                 aGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAat                                                                                                                             <  1:672838/141‑1 (MQ=255)
                       tAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGc                                                                                                                       <  1:754263/141‑1 (MQ=255)
                          ccccATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCAt                                                                                                                    >  1:416315/1‑141 (MQ=255)
                           ctcATGAATAGCGAGCCGACAAATGCTGATAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATa                                                                                                                   <  2:590385/139‑1 (MQ=255)
                           cccATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTAGGATGAGGGATTGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATa                                                                                                                   <  2:351490/141‑1 (MQ=255)
                              gtGAATAGCGAGCCTACAAAGGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTGTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGAt                                                                                                               <  2:748522/140‑1 (MQ=255)
                               tGAAAAGCGAGCCTACAGATGCTGCTAATGCAATCGCTACAAAAATGGCAACTTGGTCCATTTTCATATCTTCTGATACGCCGTTAACAACTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATa                                                                                                              <  2:25877/141‑1 (MQ=255)
                                gAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATaa                                                                                                             >  2:500566/1‑141 (MQ=255)
                                  aTAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATCTAGGCCATAGAAAAGa                                                                                                           >  2:414998/1‑141 (MQ=255)
                                  aTAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGa                                                                                                           <  2:733029/141‑1 (MQ=255)
                                  aTAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGa                                                                                                           >  1:727166/1‑141 (MQ=255)
                                        aGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCCCCTTTGTCCATTGCCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTtg                                                                                                     <  1:650652/141‑1 (MQ=255)
                                         gCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTtgt                                                                                                    >  1:77679/1‑141 (MQ=255)
                                         gCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTtgt                                                                                                    >  1:662496/1‑141 (MQ=255)
                                          ccTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGACTTTAACACTTGCTAAACTTAAAACTTGCTGGTTAACTGCTAAAATACAGTCCATAGATAAGAAAGTtgtg                                                                                                   <  2:239954/141‑1 (MQ=255)
                                                aaTGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCg                                                                                             >  1:483014/1‑141 (MQ=255)
                                                aaTGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCg                                                                                             <  1:449376/141‑1 (MQ=255)
                                                     tgctAATGCAGTCGCTACAAAAATCGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTaa                                                                                        <  2:495133/141‑1 (MQ=255)
                                                      gctAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAg                                                                                       >  1:519619/1‑141 (MQ=255)
                                                          aTGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGAttt                                                                                   >  1:317756/1‑141 (MQ=255)
                                                           tGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGAAAACAAAGTTGTGATACCGCCTAAAATTtc                                                                                  >  2:722982/1‑141 (MQ=255)
                                                           tGCAGGCGCTACAAAAATGGCACCTTGGTCCATGTTCATATCCTCTGATACGCCTTTAACACCCGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGGGATACCGCCTAAGATTtc                                                                                  <  2:248416/141‑1 (MQ=255)
                                                            gCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTtct                                                                                 <  1:628833/141‑1 (MQ=255)
                                                               gTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCtt                                                                              <  1:314903/141‑1 (MQ=255)
                                                                   cTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTaa                                                                          <  1:585672/141‑1 (MQ=255)
                                                                   cTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTaa                                                                          <  2:539767/141‑1 (MQ=255)
                                                                     aCAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAg                                                                        >  1:541693/1‑141 (MQ=255)
                                                                     aCAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAg                                                                        <  2:778762/141‑1 (MQ=255)
                                                                         aaaTGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTg                                                                    >  1:617343/1‑141 (MQ=255)
                                                                              gCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCa                                                               >  2:701180/1‑141 (MQ=255)
                                                                                aCCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCata                                                             <  1:652092/141‑1 (MQ=255)
                                                                                      gTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAACACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCCACTAAGATTTCTCTTTTAACGTTTCTTCCCCAGTTATc                                                       >  1:289789/1‑141 (MQ=255)
                                                                                          aTTTTCATATCTTCTGATACGCCTTTCACACCTGCTAAACTTAAAACTTGTGGGGTAACTGCGAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAAc                                                   <  1:569626/141‑1 (MQ=255)
                                                                                              tCATATCTTCTGATACGCCTTTAACACCTTCTAAACTTAAAACTTGCTGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGa                                               <  2:63210/141‑1 (MQ=255)
                                                                                               cATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAAATTTCTCTTTTAAAGTTTGTTACATATTTATCGAACTGGaa                                              >  2:573566/1‑141 (MQ=255)
                                                                                                 tatCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAAt                                            >  1:678143/1‑141 (MQ=255)
                                                                                                  atCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGACTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTCGAAATa                                           >  2:731022/1‑141 (MQ=255)
                                                                                                   tcttctGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTACAGTTTGTTCCATATTTATCGAACTGGCAATAg                                          >  1:504712/1‑141 (MQ=255)
                                                                                                           tACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCAc                                  >  1:405810/1‑141 (MQ=255)
                                                                                                           tACGCCTTTAACACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCAc                                  >  1:742637/1‑141 (MQ=255)
                                                                                                             cGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGa                                >  2:248519/1‑141 (MQ=255)
                                                                                                              gccgtgcacacCTGCGAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGaa                               <  2:52571/134‑1 (MQ=255)
                                                                                                              gCCTGTAACTCCTGATAAATTTAATACTTGCGGGTTAACTGCTAAAAAATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGaa                               <  1:513392/141‑1 (MQ=255)
                                                                                                               ccTTTACCACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATTCCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAAt                              <  2:139873/141‑1 (MQ=255)
                                                                                                                    aaCACCTGCTAAACTTAAAACTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATg                         >  2:577448/1‑141 (MQ=255)
                                                                                                                     acacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAATTAGTTTTTCACGAATGGATGc                        >  1:629125/1‑141 (MQ=255)
                                                                                                                      cacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGTATGCt                       >  2:15535/1‑141 (MQ=255)
                                                                                                                      cacCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATTCt                       >  2:573082/1‑141 (MQ=255)
                                                                                                                        ccTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGTTGCTcc                     >  1:426164/1‑141 (MQ=255)
                                                                                                                         cTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCt                    >  1:400240/1‑141 (MQ=255)
                                                                                                                         cTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCt                    >  1:400300/1‑141 (MQ=255)
                                                                                                                            cTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCtttt                 <  2:99598/141‑1 (MQ=255)
                                                                                                                            cTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCtttt                 >  1:401257/1‑141 (MQ=255)
                                                                                                                             tAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttt                >  2:398751/1‑141 (MQ=255)
                                                                                                                             tAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttt                >  1:251082/1‑141 (MQ=255)
                                                                                                                               aaCTTAAAGCTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTTATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttttt              <  1:251436/141‑1 (MQ=255)
                                                                                                                               aaCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttttt              >  1:48592/1‑141 (MQ=255)
                                                                                                                               aaCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttttt              <  1:27529/141‑1 (MQ=255)
                                                                                                                               aaCTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTGATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCttttttt              <  1:237810/141‑1 (MQ=255)
                                                                                                                               aaCTTAAAACTGGCGGGTTAACTGCTAAAATATAGGCCATAGATTAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCtttttt               <  1:222681/141‑1 (MQ=255)
                                                                                                                                aCTTAAAACGTGCAGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTTATACCGTCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTa             <  2:237281/141‑1 (MQ=255)
                                                                                                                                 cgtAAAACTTGCTGGTTAACTGCTAATATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTACAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaa            <  2:476207/139‑1 (MQ=255)
                                                                                                                                  ttAAAACTTGCGGTTTAACTGCTAAAATTTCGATCATACATAAGTAAGTTGTTATACCGCCGAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaaa           <  1:108541/141‑1 (MQ=255)
                                                                                                                                  ttAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGGAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaaa           <  1:743361/141‑1 (MQ=255)
                                                                                                                                  ttAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTaaa           <  1:121260/141‑1 (MQ=255)
                                                                                                                                    ataaCTTGCTTGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAAt         <  1:102033/139‑1 (MQ=255)
                                                                                                                                    aaaaCTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAAt         <  1:303603/141‑1 (MQ=255)
                                                                                                                                    aaaaCTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTAGAAATAGTTTTTCACGACTAGATGCTCCTTTTTTTAAAAt         >  2:258242/1‑141 (MQ=255)
                                                                                                                                    aaaaCTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCACATTTATCGAACTAGAAATAGTTGTTCACGCATGGATGCTCCTTTTTTTTAAAt         >  1:106714/1‑141 (MQ=255)
                                                                                                                                    aaaaCTTGCGGGTTAACTGCTAAAAAAAAGCCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTCTAAATAGTTTTTCACGAATGGATGCTCCTCTTTCTAAAAt         >  2:246033/1‑141 (MQ=255)
                                                                                                                                      aaCTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATAc       >  2:243875/1‑141 (MQ=255)
                                                                                                                                       aCTTACGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTAACGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATAcc      <  2:519619/141‑1 (MQ=255)
                                                                                                                                        cTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCt     <  1:508632/141‑1 (MQ=255)
                                                                                                                                           gCGGGTTAACTGATAAAATATAGGTCATAGATAAGATAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATAGATGCTAAGTTTTTTAAAATACCTTaa  <  2:684543/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
AACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAA  >  USA300TCH1516_ALE/2373528‑2373807

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: