Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,646,189 G→C A155G (GCA→GGA)  ypdF ← Aminopeptidase YpdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,646,1890GC97.7% 300.5 / NA 88A155G (GCA→GGA) ypdFAminopeptidase YpdF
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (37/48);  total (38/50)
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.99e-01

TTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACGG  >  USA300TCH1516_ALE/1646053‑1646326
                                                                                                                                        |                                                                                                                                          
tgtctCACGTGCAACACCATGTGGTAATGCACCTTTAGGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                                        <  2:51737/137‑1 (MQ=255)
 ttATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCtttt                                                                                                                                       <  1:623984/141‑1 (MQ=255)
   aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                                     >  2:326454/1‑141 (MQ=255)
   aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                                     >  2:874194/1‑141 (MQ=255)
   aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                                     <  1:947372/141‑1 (MQ=255)
     cACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGCTGCTACAATCGTATCGAATGATGGTCCATCTGCTCATAATTCTAGCATTTTGCTTTCTAATATTGACTTCAATTCTTTTTCGGTCAGGCCTCCTTTTacca                                                                                                                                   >  2:526551/1‑139 (MQ=255)
      aCTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                                  <  2:616023/141‑1 (MQ=255)
      aCTAGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                                  <  2:845648/141‑1 (MQ=255)
         ggCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTAGCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTTTTTTTCAGTCATGCCGCCTTTTACATCAGt                                                                                                                               <  2:604682/140‑1 (MQ=255)
          gCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGtt                                                                                                                              <  1:934202/141‑1 (MQ=255)
          gCAACACCATCTGGTAATGCCCCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGtt                                                                                                                              <  2:899018/141‑1 (MQ=255)
           cAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTa                                                                                                                             <  1:608536/141‑1 (MQ=255)
              cacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                          <  1:700454/141‑1 (MQ=255)
              cacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                          <  2:363877/141‑1 (MQ=255)
              cacCATGTGGTAACGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGTCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                          <  1:193344/141‑1 (MQ=255)
              agccATGTGGTAATGCACCCGTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCACTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                          <  2:382322/139‑1 (MQ=255)
               acCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACACCAGTTAAAat                                                                                                                         >  1:806832/1‑141 (MQ=255)
                ccgtgtgGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                        <  1:866510/138‑1 (MQ=255)
                      ggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCa                                                                                                                  <  1:857587/141‑1 (MQ=255)
                      ggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCa                                                                                                                  <  2:273241/141‑1 (MQ=255)
                       gTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                                                                                 <  1:77609/141‑1 (MQ=255)
                       gTAATGCACCTCCATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                                                                                 <  1:303412/141‑1 (MQ=255)
                           tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                             <  2:845483/141‑1 (MQ=255)
                            gCACCTCTAGGACCACATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGtt                                                                                                            <  2:33532/141‑1 (MQ=255)
                              aCCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtc                                                                                                          >  1:788968/1‑141 (MQ=255)
                               ccTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtca                                                                                                         <  2:138856/141‑1 (MQ=255)
                               ccTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtca                                                                                                         <  1:984770/141‑1 (MQ=255)
                               ccTCTATGACCAGATGCTACAATCGTATCGAAAGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtct                                                                                                         >  1:713146/1‑140 (MQ=255)
                                  ctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATCTGTTtcatca                                                                                                      >  2:708632/1‑141 (MQ=255)
                                        ccAGATGCTACATTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAata                                                                                                <  2:26494/141‑1 (MQ=255)
                                        ccAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAata                                                                                                <  1:487675/141‑1 (MQ=255)
                                                       gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGa                                                                                >  2:182850/1‑141 (MQ=255)
                                                         aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAt                                                                              >  2:321078/1‑141 (MQ=255)
                                                           cGAATGATGGTCCATCTGCTCATAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCGTCAACAATGTTAGCGGCTGGTTGAtttt                                                                            >  2:682661/1‑140 (MQ=255)
                                                              atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCACCAATATTAGCTGCTTTTTGAATTagag                                                                         >  1:948676/1‑141 (MQ=255)
                                                              atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGATTTaaac                                                                         >  2:928776/1‑140 (MQ=255)
                                                               tgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGc                                                                        >  1:370408/1‑141 (MQ=255)
                                                                     tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTc                                                                  >  1:818549/1‑141 (MQ=255)
                                                                     tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTACTTTTTGAATTAAAGCAATTTc                                                                  >  1:545902/1‑141 (MQ=255)
                                                                      ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTct                                                                 >  2:444817/1‑140 (MQ=255)
                                                                      ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCACTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                                 >  2:116784/1‑141 (MQ=255)
                                                                      ccATCTGCTCCTAAGTCTAGCATTTTGCTTTCTATAAGTGCCTTTACTTCTTCGTCAGGCATCCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                                 >  2:453025/1‑141 (MQ=255)
                                                                          cTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAg                                                             >  1:612792/1‑141 (MQ=255)
                                                                           tGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCCGc                                                            >  1:288751/1‑141 (MQ=255)
                                                                            gccccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCCTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                                           <  1:87194/138‑1 (MQ=255)
                                                                               ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt                                                        <  1:874194/141‑1 (MQ=255)
                                                                               ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt                                                        <  1:326454/141‑1 (MQ=255)
                                                                                 tcaTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAATAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCttt                                                      <  1:830841/139‑1 (MQ=255)
                                                                                 tAATTCTATCCTTTTGCTTTCTCATCTTGCCTTCCATTATTTTTGAGTCGTGTCTCTTTTTACGGCCCCTAAAATTTTTTCGGAGGTTGCATCAACAATATGATCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCttt                                                      <  2:477094/141‑1 (MQ=255)
                                                                                   attgtcgCAGTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACATTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGa                                                    <  2:823075/135‑1 (MQ=255)
                                                                                   aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGa                                                    >  2:13647/1‑141 (MQ=255)
                                                                                   aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTAAGCATCTTTGa                                                    >  2:61485/1‑141 (MQ=255)
                                                                                        aGCATTTTGCTTTCTAATCTTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtct                                               <  1:275416/141‑1 (MQ=255)
                                                                                        aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtct                                               >  2:527279/1‑141 (MQ=255)
                                                                                         gCATTTTGCTTTCTACTATTGCCTTTAATTCTATTTCAGAAGTTGCTGCTTTTACAACAGTTAAAATATATTCATATGTTTTATCCACAATATTAGCTGCTTTTTGAATTAAAGTAATTTCGTCTGGATTTATGACGtctc                                              <  2:449601/141‑1 (MQ=255)
                                                                                         ccATTTTGTTTTCTAATATTGCCTTTAATTCTTTTTCATTCATTCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTATCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                              <  2:765555/140‑1 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTTCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                             <  1:420035/141‑1 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                             <  2:119290/141‑1 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                             <  2:64048/141‑1 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTCCGTCAGCATCTTTGACGTctct                                             <  2:727282/141‑1 (MQ=255)
                                                                                          cATTTTGCTTTCTAATATTGCCTTTAACTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                             <  2:277933/141‑1 (MQ=255)
                                                                                           aTTTTGCTTTCTAATATTGCCTTTACTTCTTTTTCCGTCATGCCTCCTTTTACAACAGTTAAAATTTATTCATATGTTTCATCAACAATATTAGCTTCTTTTGAAATTAAAGCACTTTAGTAAGCATCTTTGACGtctttt                                            >  2:277290/1‑138 (MQ=255)
                                                                                           aTTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTa                                            >  2:405516/1‑141 (MQ=255)
                                                                                           aTTTTGCTTTCTAATATTGCCTTTAATTATTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGAAATTTGGTCAGCATCTTTGACGTCTCTa                                            >  2:500898/1‑141 (MQ=255)
                                                                                                 cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTa                                      <  1:319099/141‑1 (MQ=255)
                                                                                                 cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTa                                      <  1:372651/141‑1 (MQ=255)
                                                                                                  tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCCACAGTATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTAt                                     <  2:248961/141‑1 (MQ=255)
                                                                                                  tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTAt                                     <  1:981464/141‑1 (MQ=255)
                                                                                                  tttCTAATATTCCCCTTACTGCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCCTTGACGTCTCTAATTTTAt                                     <  2:89543/141‑1 (MQ=255)
                                                                                                    tCTATTATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCt                                   <  2:868842/141‑1 (MQ=255)
                                                                                                     cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGATTTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTa                                  >  1:673902/1‑141 (MQ=255)
                                                                                                     cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTTTCTa                                  >  2:48717/1‑141 (MQ=255)
                                                                                                     cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTa                                  >  1:192364/1‑141 (MQ=255)
                                                                                                     cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTa                                  >  1:630890/1‑141 (MQ=255)
                                                                                                     cTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAACTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTa                                  >  2:511774/1‑141 (MQ=255)
                                                                                                            tGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAtt                           <  2:788968/141‑1 (MQ=255)
                                                                                                             gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTTACGTCTCTAATTTTATTTACAGTAtt                           >  2:74828/1‑141 (MQ=255)
                                                                                                             gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTCCAGTAtt                           >  2:977792/1‑141 (MQ=255)
                                                                                                             gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAtt                           >  2:58792/1‑141 (MQ=255)
                                                                                                                    aTTGTTTTTCAGGCATGCTTTCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGc                   <  2:491119/141‑1 (MQ=255)
                                                                                                                    aTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGc                   <  1:522879/141‑1 (MQ=255)
                                                                                                                         ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatt              >  1:689606/1‑141 (MQ=255)
                                                                                                                         ttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatt              >  1:198427/1‑141 (MQ=255)
                                                                                                                            tCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAAt           <  2:872885/141‑1 (MQ=255)
                                                                                                                             cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATg          >  2:757521/1‑141 (MQ=255)
                                                                                                                                  aTGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGatat      >  1:685025/1‑140 (MQ=255)
                                                                                                                                  aTGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGatat      <  2:685025/140‑1 (MQ=255)
                                                                                                                                    gCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACg   <  1:807200/141‑1 (MQ=255)
                                                                                                                                     cctcctTTGACAACAGTTATAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACgg  <  2:437148/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATATACGG  >  USA300TCH1516_ALE/1646053‑1646326

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: