Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,504,989 T→C 9.0% intergenic (+293/+111) USA300HOU_RS12890 → / ← gltT hypothetical protein/Proton/sodium‑glutamate symport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,504,9890TC 9.0% 300.4 / 7.6 140intergenic (+293/+111)USA300HOU_RS12890/gltThypothetical protein/Proton/sodium‑glutamate symport protein
Reads supporting (aligned to +/- strand):  ref base T (58/62);  new base C (3/15);  total (62/78)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

TTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGACCTTT  >  USA300TCH1516_ALE/2504850‑2505129
                                                                                                                                           |                                                                                                                                             
ttGCGCATGCATAAACGCCCCTTCCAACCTGAGGTCACTGCGCTTCCGTTCGGCTGGGCTTTAGCGCATGCATAAAAGCCCCTAACACCCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTc                                                                                                                                             <  2:36124/141‑1 (MQ=14)
  gcgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTccc                                                                                                                                           >  2:755482/1‑141 (MQ=255)
   cgctcgcatTAACGCCCCTATCATCCTGAGGCCACCCCGCGCCGGTTCGCCTGCGCTTTAGCGCATGCATAAACGCCCCTAACACCCGGAGGTCACTCCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCt                                                                                                                                          <  1:104779/136‑1 (MQ=14)
   cgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCt                                                                                                                                          <  1:562407/141‑1 (MQ=255)
    gcatgcatCAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTa                                                                                                                                         <  1:141176/141‑1 (MQ=255)
     catgcatAAGAGCCCCTATCCACCTGGGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaa                                                                                                                                        <  2:571479/141‑1 (MQ=255)
     catgcatAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCGCTACGCGCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaa                                                                                                                                        >  1:770940/1‑141 (MQ=255)
     catgcatAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaa                                                                                                                                        >  1:276673/1‑141 (MQ=255)
     catgcatAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCATTAGCACATGCATAAAAGCCCCTAACAACCGGAGGTCACTACGCGTCGGGTCGCCTCCGCATTACCGCACACATAAAAACCCCTaa                                                                                                                                        >  2:242204/1‑141 (MQ=11)
      aagcatTAAAGCCCCTAACAACCTGAGTGCACGTCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaac                                                                                                                                       >  2:935948/3‑141 (MQ=14)
        gcatAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaacaa                                                                                                                                     >  1:887676/1‑141 (MQ=255)
             aaaGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGa                                                                                                                                >  2:672379/1‑141 (MQ=255)
               aGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAgg                                                                                                                              >  1:284060/1‑141 (MQ=255)
               aGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAgg                                                                                                                              >  1:249673/1‑141 (MQ=255)
                      aacaacCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACACCCTAAGGTCACTAc                                                                                                                       >  2:503147/1‑141 (MQ=255)
                       acaacCTGAGGTCACTGCGCCCCGGGTCGCCCGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACg                                                                                                                      <  2:380367/141‑1 (MQ=255)
                        caacCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACgg                                                                                                                     >  1:868598/1‑140 (MQ=255)
                             tGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTTCGCTTcct                                                                                                                >  2:295491/1‑140 (MQ=255)
                              gAGGTCACTGCGCTCCGCGTCGCCTTGGCTTTAGCGCATGCATAAAGGCCCCTAACGACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGt                                                                                                               <  2:173199/141‑1 (MQ=17)
                               gggTCCCTGCGCCCGCGTTCGCCTGCGCTTCAGCGCATGCATAAAAGCCCCTAACAACCGTAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAACAGCCCCAATCAACCTGAGGTCACTACGCTTCCGtt                                                                                                              <  1:224855/140‑1 (MQ=14)
                               aGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGCGGTCACTAAGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGTTCACTGCGCTTCCGtt                                                                                                              >  1:102218/1‑141 (MQ=255)
                                ggTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACACCCTGGGGTTACTAcggttcggttt                                                                                                             >  2:548854/1‑140 (MQ=17)
                                ggTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTc                                                                                                             >  2:708690/1‑141 (MQ=255)
                                ggTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTAACTACGCTTCCGTTc                                                                                                             >  1:804462/1‑141 (MQ=255)
                                 gTCACTGCGCTCCGGTTCGACTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAATCTGAGGTCACTACGCTTCCGTTCg                                                                                                            <  1:601059/141‑1 (MQ=255)
                                   cACTGCCCTCCGGGTCGCCCGGGCTTTCGCGCATGCATACAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGcc                                                                                                          <  2:282630/141‑1 (MQ=17)
                                      tGCGCCCCGGTTCGCCTGCGCTTTAGCCCACGCATAAAAGCCCCTAACCTCCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgc                                                                                                       <  2:465593/141‑1 (MQ=17)
                                       ccgcTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgcg                                                                                                      <  1:672279/140‑1 (MQ=255)
                                          cTCCGGTTCGCCTGCGCTTTCGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGATTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCtt                                                                                                   <  2:702620/141‑1 (MQ=255)
                                          cTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCGCTAACAACCTGAGGTAACTACGCTTCCGTTCGCCTGCGCtt                                                                                                   >  1:530784/1‑141 (MQ=255)
                                          cTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACGACGCTTCCATTCGCCCgcgcat                                                                                                   >  2:27778/1‑139 (MQ=255)
                                               gTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAgcg                                                                                              >  1:461025/1‑141 (MQ=255)
                                                 tCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCCCCTGCGCTTTAGCgca                                                                                            <  1:485205/141‑1 (MQ=255)
                                                  cGCCTTCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACAACGCTTCCGTTCGCCTGCGCTTTAGCgcat                                                                                           <  2:697362/141‑1 (MQ=255)
                                                    ccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACGACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAATGTCCCTACCACCCTGTGGTGCCCACGCTTTCGTTCGCCTGCGCTTTAGCgcttgc                                                                                         >  1:513661/1‑141 (MQ=14)
                                                    ccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgg                                                                                         >  1:376040/1‑140 (MQ=255)
                                                    ccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTCCCGTTCGCCTGCGCTTTAGGgcatga                                                                                         >  2:736909/1‑140 (MQ=255)
                                                     gtgggcTTTAGCGTATGTATAGGGGCCACCAGCAATCTGAGGTCCCTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgca                                                                                        <  2:138482/137‑1 (MQ=14)
                                                      tGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgcat                                                                                       <  1:451068/141‑1 (MQ=255)
                                                      tGCCCTTTAGCGCATTCATAAATGCCCCTAACAACCGGGGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgcat                                                                                       <  2:900730/141‑1 (MQ=14)
                                                        cgcTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTTCGCATTAGCGCTTGCATaa                                                                                     >  2:260159/1‑141 (MQ=255)
                                                          cTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaaa                                                                                   <  1:395973/141‑1 (MQ=255)
                                                             tAGCGCATGCATAACAGCGCCTAACCACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcc                                                                                <  1:908969/141‑1 (MQ=255)
                                                               gcgcATTCATAAAAGCCCCTTACAACCTGAGGTCACTCCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGcccc                                                                              <  2:873555/141‑1 (MQ=255)
                                                                cgcATGCATAAAAGCCCCTAACAACCTGGGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCt                                                                             <  1:163076/141‑1 (MQ=255)
                                                                cgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTGTAGCGCATGCATAAAAGCCCCt                                                                             >  1:652299/1‑141 (MQ=255)
                                                                     gcatAAAAGCCCCTAACAACCTGAGGGCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaa                                                                        <  1:97084/141‑1 (MQ=255)
                                                                     gcatAAAAGCCCCTAACAACCGGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaa                                                                        <  1:359880/141‑1 (MQ=255)
                                                                       atAACAGCCCCTTAGAATCTGAGGGCACTACGCTTCGGCTCCACTGCGCTTTAGCGCATGCATAAATGTCCCTAACAGCCTTAGTTCACGACGCTTCCGTTCGCCTGCGCTTTATCGCATGCATAAAAGCCCCTAACAAcc                                                                      <  2:564971/141‑1 (MQ=14)
                                                                          aaaGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCaaa                                                                   >  1:356132/1‑141 (MQ=255)
                                                                           aaGCCCCTAACACCCGGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAg                                                                  <  1:277165/141‑1 (MQ=255)
                                                                           aaGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAACCCCCTAACAACCAAAg                                                                  <  1:469751/141‑1 (MQ=255)
                                                                                ccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGt                                                             >  2:131457/1‑141 (MQ=255)
                                                                                ccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTCCTCTTTATCGCATGCATAAAAGCCCATACACACCAAAGGCTGt                                                             >  2:27541/1‑141 (MQ=14)
                                                                                acTAACAACCGGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTACCGTTCGCCTGCGCTTTACCGCTTACATAAAAGCCCCTAACTACCAACGATTgg                                                             >  2:817402/2‑140 (MQ=14)
                                                                                 cTAACCACCTGTGGTCACTCCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTa                                                            <  2:774622/141‑1 (MQ=255)
                                                                                 cTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTa                                                            <  2:164849/141‑1 (MQ=255)
                                                                                  caactacCTGGGGTCACTACGCTTCGGTTCCCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTaa                                                           <  2:57038/140‑1 (MQ=255)
                                                                                  caacaccCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTaa                                                           <  2:502361/140‑1 (MQ=255)
                                                                                  caacaacCTGAGGTCCCTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTaa                                                           <  2:818320/140‑1 (MQ=255)
                                                                                   aacaacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTAGGCATTAGCGCATGCATCAAAGCCCCTAACAACCCAAGGTTGTAAg                                                          >  2:326801/1‑141 (MQ=17)
                                                                                    acaacCTGACGTCCATACGCTTCGGTTCGCCTGAGCTCTAGCGCATGCATAAAAGACCCTAACAACCTGAGGTCACTACGCGTCCGAACGCCTACCCATTAGCGCATGCATCAAACCCACAGACAACTAAACGTTGtcgtt                                                         >  2:672671/1‑137 (MQ=9)
                                                                                     caacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAggg                                                        <  1:727929/141‑1 (MQ=255)
                                                                                      aacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAgggg                                                       >  1:478657/1‑141 (MQ=255)
                                                                                       acCTTAGGTCACCACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc                                                      <  2:233951/141‑1 (MQ=255)
                                                                                       acCTGGGGTCACTCCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGGCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc                                                      <  2:704852/141‑1 (MQ=255)
                                                                                        ccTGAGGTCACTGCGCTTCGGTTCGCCTGCCCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGAATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCt                                                     <  2:312634/141‑1 (MQ=255)
                                                                                             ggTCCCTACGCTTCGGTTCGGCCGCGCTTTAGCGAATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCGTCCGTTCGCCTGCGCTGTAGCGCAGGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTaaa                                                <  1:246779/141‑1 (MQ=17)
                                                                                               tCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGTGCATGCATAAAAGCCCCTTACAACCAAAGGTTGTATGGGGCTTTAAATg                                              >  1:449706/1‑141 (MQ=255)
                                                                                               tCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAACTg                                              >  1:325005/1‑141 (MQ=255)
                                                                                                cACTACGCTCCGGTCGCCCTTCGCTTTAGCGCATGCATAAAAGTCCCTAACACCCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGt                                             <  1:876938/141‑1 (MQ=18)
                                                                                                 aCTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTc                                            >  1:793291/1‑141 (MQ=255)
                                                                                                  cTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCa                                           >  1:104063/1‑141 (MQ=255)
                                                                                                         ggtgTTCGCCTGCGCTGTAGCGCATGCATAAAAGTCCCTTACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACCACGAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTa                                    <  2:439511/138‑1 (MQ=18)
                                                                                                          cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTTCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAt                                   <  1:603345/141‑1 (MQ=255)
                                                                                                          cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAt                                   <  2:607372/141‑1 (MQ=255)
                                                                                                            gTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAttt                                 >  1:107095/1‑141 (MQ=255)
                                                                                                            gTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAttt                                 >  1:741287/1‑141 (MQ=255)
                                                                                                               ccccTGCGCTTTAGCTCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAtgtttt                              <  2:575945/139‑1 (MQ=255)
                                                                                                               cGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAtttttt                              <  1:755482/141‑1 (MQ=255)
                                                                                                                  cTGCGGTTTATCGCATGCAGAAGAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTaaa                           <  2:744223/141‑1 (MQ=21)
                                                                                                                  cTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTaaa                           <  2:790929/141‑1 (MQ=255)
                                                                                                                  cTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTaaa                           <  1:423061/141‑1 (MQ=255)
                                                                                                                  cTGCGCTTGAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTaaa                           <  1:107614/141‑1 (MQ=255)
                                                                                                                     cgcTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTg                        >  1:718815/1‑141 (MQ=255)
                                                                                                                         ttAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTCAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGAttt                    >  2:515186/1‑141 (MQ=255)
                                                                                                                           tgcgcATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTa                  <  2:770940/140‑1 (MQ=255)
                                                                                                                                 tgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATAtt            >  2:774552/1‑141 (MQ=255)
                                                                                                                                    atAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCtt         >  2:933948/1‑141 (MQ=255)
                                                                                                                                    atAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGCTGTGTTTTATTTTAAAAGACGGGTAGTGTTTTTTAAATTGATTTTAAATATTCtt         >  1:443676/1‑141 (MQ=17)
                                                                                                                                    atAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCtt         <  2:215745/141‑1 (MQ=255)
                                                                                                                                     tAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTg        <  1:6255/141‑1 (MQ=255)
                                                                                                                                     tAAAAGTCCCTAACAACCTGAGGTCACTACGCGTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTg        >  2:455943/1‑141 (MQ=255)
                                                                                                                                      aaaaGTGCCCTAACATCCTGAGGGCACGACGCTTCCGGTCGGCTTCGCTTTAGCGCATGCATAAAAGCCCCGAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTg        <  1:902306/141‑1 (MQ=37)
                                                                                                                                      aaaaGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGa       <  1:461816/141‑1 (MQ=255)
                                                                                                                                      aaaaGTCCCCAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGa       <  1:647008/141‑1 (MQ=255)
                                                                                                                                      aaaaGCCCCTAACAACCTGCGGTCACTACGCGTCCGTTCGCCTGAGATTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGa       <  2:614274/141‑1 (MQ=37)
                                                                                                                                      aaaaGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGa       <  2:447421/141‑1 (MQ=255)
                                                                                                                                        aaGCCCCTAACCACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGAcc     <  2:42697/141‑1 (MQ=255)
                                                                                                                                           tCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCACATAGTATTTTTTAAATTGATTTTAAATATTCTTGACCttt  <  1:672379/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
TTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTATTTTTTAAATTGATTTTAAATATTCTTGACCTTT  >  USA300TCH1516_ALE/2504850‑2505129

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