Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 4,449 C→A 6.2% Q171K (CAA→AAA)  recF → DNA replication and repair protein RecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE4,4490CA 6.2% 362.4 / 2.4 129Q171K (CAA→AAA) recFDNA replication and repair protein RecF
Reads supporting (aligned to +/- strand):  ref base C (63/55);  new base A (3/7);  total (67/62)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 132 total aligned reads displayed.

AAAAGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGA  >  USA300TCH1516_ALE/4310‑4588
                                                                                                                                           |                                                                                                                                           
aaaaGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCa                                                                                                                                            >  1:144555/1‑141 (MQ=255)
 aaaGGCTCTCCTCAAATACGACAACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaa                                                                                                                                           <  1:712817/141‑1 (MQ=255)
   aGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCaaaa                                                                                                                                         <  2:717797/141‑1 (MQ=255)
    ggCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTCCTTAAAGCAGTTACAATTAGGCCaaaaa                                                                                                                                        >  2:241823/1‑141 (MQ=255)
         tcctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGa                                                                                                                                   >  1:1012497/1‑141 (MQ=255)
          cctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATCCAAACGTAGTTGAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAACAAGGAc                                                                                                                                  >  1:504309/1‑141 (MQ=255)
          cctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGAc                                                                                                                                  >  1:403407/1‑141 (MQ=255)
          cctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGAc                                                                                                                                  >  1:1054547/1‑141 (MQ=255)
          cctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAATTACAATTAGGCCAAAAAAAGGAc                                                                                                                                  >  1:1023526/1‑141 (MQ=255)
           ctcAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAGAAAGGACt                                                                                                                                 >  2:725486/1‑141 (MQ=255)
             cAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTa                                                                                                                               >  1:323047/1‑141 (MQ=255)
               aaTAAGACGACGCTTTATAGATATGGTGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTaaa                                                                                                                             >  2:657872/1‑140 (MQ=255)
                  aaGACGACGCTTTATAGAAATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAATAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAAt                                                                                                                          <  2:616651/141‑1 (MQ=255)
                   aGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATg                                                                                                                         <  1:534354/141‑1 (MQ=255)
                    gacgacGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGt                                                                                                                        >  1:1004005/1‑141 (MQ=255)
                    cacgacgCTTTATAGATATGGAGTTGCGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGt                                                                                                                        <  1:516672/140‑1 (MQ=255)
                     acgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTaaacatgat                                                                                                                       >  2:898191/1‑134 (MQ=255)
                     acgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGtt                                                                                                                       >  2:457624/1‑141 (MQ=255)
                     acgacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAACCACGATTATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGtt                                                                                                                       >  1:1040718/1‑141 (MQ=255)
                       gacgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTgg                                                                                                                     <  1:677415/141‑1 (MQ=255)
                         cgCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGaa                                                                                                                   <  1:585194/141‑1 (MQ=255)
                            cttATAGATAAGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTCTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTa                                                                                                                <  1:340597/140‑1 (MQ=255)
                                taGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCGAAAAAAAGGACTTAACAATGTTGGAAGTATTaa                                                                                                            >  2:245109/1‑141 (MQ=255)
                                taGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaa                                                                                                            >  1:557542/1‑141 (MQ=255)
                                 aGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTaaa                                                                                                           <  1:51614/141‑1 (MQ=255)
                                      tgaaGTTGGTCCAAATTTCTTCTGTTTACTGCAAAGATTTAGCTCAATACCCACGAATTTTAAAGCAAAAGAATAATTCCTTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagc                                                                                                      <  2:616948/138‑1 (MQ=255)
                                      tGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagc                                                                                                      >  1:478200/1‑141 (MQ=255)
                                      tGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagc                                                                                                      <  2:945670/141‑1 (MQ=255)
                                      tGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAAGGTTGGAAGTATTAAATcagc                                                                                                      >  1:577900/1‑141 (MQ=255)
                                      gggAGTGGGGCCAAGTTTCTGCTGTTTACTTAAATGATGAAGCTCAATACCAAGGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAATTAGGACAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagc                                                                                                      <  2:481019/140‑1 (MQ=255)
                                       ggAATTAGGCCAAATTTCTGCTGATTACTTAAAAGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagca                                                                                                     <  2:916888/141‑1 (MQ=255)
                                        gAGTTGGGCTAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagcag                                                                                                    <  1:263085/141‑1 (MQ=255)
                                        gAGTTGGCCCAAATTTCTGCTGTTAAATTAAATAATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATcagcag                                                                                                    <  2:362069/141‑1 (MQ=255)
                                         ggTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  <  2:132244/140‑1 (MQ=255)
                                         aGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                                                                   <  1:367448/141‑1 (MQ=255)
                                         aGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGt                                                                                                   <  2:429726/141‑1 (MQ=255)
                                          gTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  <  1:129429/141‑1 (MQ=255)
                                          gTTGGGCCAAATTTCAGTTGTTTACTGAAATGATCTCGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGTAGTTACAATTAGGACAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGtt                                                                                                  <  1:643237/141‑1 (MQ=255)
                                            tGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTg                                                                                                <  2:444661/141‑1 (MQ=255)
                                            tGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTCTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTg                                                                                                <  2:583849/141‑1 (MQ=255)
                                             gggCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCGGTTACCATTAGGCCAAAACAAGGACTTAACAGCGTTGGCAGTATTAAGTGAGCAGTTTGc                                                                                               >  1:161467/1‑141 (MQ=255)
                                               gCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTg                                                                                             >  1:399665/1‑141 (MQ=255)
                                                ccAAATTTCTGCTGTTTACAGAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGa                                                                                            >  1:689207/1‑141 (MQ=255)
                                                       tCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCa                                                                                     >  1:279217/1‑141 (MQ=255)
                                                       tCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCa                                                                                     >  1:959880/1‑141 (MQ=255)
                                                        ctgctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTCCAATTAGGCAAACAAAAGGACTTAACAATGTTGAAAGTATTAAATGAGCAGTTTTCTGAATTTGCaa                                                                                    >  1:402644/1‑141 (MQ=255)
                                                          ggtgtgTACTTAAATGACTTAGCTCAAGACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATg                                                                                  <  2:893125/139‑1 (MQ=255)
                                                          gctgTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCACTg                                                                                  >  1:232289/1‑141 (MQ=255)
                                                                  cTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAAc                                                                          >  1:276803/1‑141 (MQ=255)
                                                                           tttAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACg                                                                 >  1:621397/1‑141 (MQ=255)
                                                                            ttAGCTCAATACCAACGTATTTTAAAGCACAAGCATTTTTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGt                                                                >  1:672550/1‑141 (MQ=255)
                                                                            ttAGCTCAATAACAACGTATTTTAAAGCAAAACAATAATTACTTAAAGCAGTTACAATAAGGCAAAAAAAAGGACTTCACAACCTTGGAAGTATTAAAGCAGCAGTTTACTGAATATGCAATGAACGCACCTGATAGACGt                                                                >  2:638636/1‑141 (MQ=255)
                                                                              aGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCGGTTCTAATTAGGCCAAAACAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGc                                                              <  1:993665/141‑1 (MQ=255)
                                                                              aGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCGGCAGTTTGATGAATATACAATGAAATTAAATGATAAACGTGc                                                              >  2:205866/1‑141 (MQ=255)
                                                                              aGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAAGGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGc                                                              >  1:522684/1‑141 (MQ=255)
                                                                              aGCTCAATACCAACGTAATTTAAAGAAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTAGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGc                                                              <  2:691982/141‑1 (MQ=255)
                                                                              aGCTAAATACCAACGTATTTTAAAGCAAAAGAGTAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTGACTGATAACCGTGc                                                              >  1:759818/1‑141 (MQ=255)
                                                                                cTCAATAACAACTTTATTTAAAGCAAAAGAATAATTACGTAAAGCAGTTAAAAAAAAAAAAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGcac                                                            <  2:602281/141‑1 (MQ=255)
                                                                                  cAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCCCAt                                                          >  1:301509/1‑141 (MQ=255)
                                                                                  cAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAACAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACTTGCACAt                                                          >  1:726960/1‑141 (MQ=255)
                                                                                      aCCAACGTATTTTAAAGCAAAAGAACAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTa                                                      <  1:615863/141‑1 (MQ=255)
                                                                                       caaaCGTGTTTTAAAGCAAAAAAATAAGTACTTAAAGCAGTTACAATTAAGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAt                                                     <  2:176673/139‑1 (MQ=255)
                                                                                        cAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAtt                                                    <  2:689207/141‑1 (MQ=255)
                                                                                        cAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGACCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTAtt                                                    >  2:748947/1‑141 (MQ=255)
                                                                                        cAACGTAGTTCAAAGCAACAGAATAATTACCTCAAGCAGTTACAAGTAGGCCAAAAAAAGGACTTAACAATGTTGGACGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGGTAAACGTGCACATTTTAtt                                                    <  2:108204/141‑1 (MQ=255)
                                                                                           cGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCaa                                                 >  2:512160/1‑141 (MQ=255)
                                                                                             tATTTTAAAGCAAAAGAATAATTACTTAAAGCAGGTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAaga                                               <  2:975940/141‑1 (MQ=255)
                                                                                              aTTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAagag                                              <  2:1004005/141‑1 (MQ=255)
                                                                                               ttttGAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAAGAATTAACAATGTAGGCAGTATTAAAACAGCAGTTTGCTGAGTATGCAATGAAAGTAACTGTAAAACGTGCACCTTTTATTccaaagg                                             >  2:177784/1‑135 (MQ=255)
                                                                                               ttttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGc                                             <  2:504309/141‑1 (MQ=255)
                                                                                                 ttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGGCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTa                                           <  1:457624/141‑1 (MQ=255)
                                                                                                 ttAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTa                                           <  2:395204/141‑1 (MQ=255)
                                                                                                  tAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTag                                          >  1:1007976/1‑141 (MQ=255)
                                                                                                   aaaGCAAAAGAATAATCACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTaga                                         <  1:861543/141‑1 (MQ=255)
                                                                                                   aaaGCAAAAGAATAAATACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACTTGCACATTTTATTCAAGAGCTaga                                         <  2:66343/141‑1 (MQ=255)
                                                                                                     aggaaaaGAATAATTACTTAAAGCAGTTACAATTAGGTCAAAAAATGGACTTAACAATGTTTGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGt                                       <  2:59733/138‑1 (MQ=255)
                                                                                                         aaaGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGtt                                   >  1:116356/1‑141 (MQ=255)
                                                                                                           aGAATCATTACTTAAAGCAGTTAAAATTAGGCAAAAAAAAGGACTTAACAATGTTGGAAGGATTAAATCAGCAGTTTGCTGAATATTCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAg                                 <  2:26640/141‑1 (MQ=255)
                                                                                                               taatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAACGAATGCTACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTaa                             >  1:1038214/1‑141 (MQ=255)
                                                                                                               aaaatACTTACAGCAGTTACAATTAAGACAAAAAAAGGACTTAACAATGTTGGAAGTATGAAATCAGCAGTTTGCTGAATATCCAACGAAAGAAACTGATAAACGTCCACATTTTATTCAAGAGCTAGAGTCGTTAGCTaa                             <  1:836853/137‑1 (MQ=255)
                                                                                                                aatTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTaaa                            >  1:426139/1‑141 (MQ=255)
                                                                                                                 atTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATACAATGAAAATAACTGATAAACGTGAACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAAc                           >  2:114304/1‑141 (MQ=255)
                                                                                                                  ttACTTAAAGCTGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGACAGTAACTGATAAACGTGCACATTTTATTCAGGAGCGATTTTCGTTAGCTAAAcc                          <  2:1025378/141‑1 (MQ=255)
                                                                                                                  ttACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAAcc                          <  1:949265/141‑1 (MQ=255)
                                                                                                                    aCTTAGAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGa                        <  1:638246/141‑1 (MQ=255)
                                                                                                                        aaaGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGGCGTTAGCTAAACCGATTCa                    >  1:1031708/1‑141 (MQ=255)
                                                                                                                         aaGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCAt                   <  2:403407/141‑1 (MQ=255)
                                                                                                                           gCAGTTACAATTAGGTCAAAAAAAGTACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTATACCGATTCATGc                 <  2:129660/141‑1 (MQ=255)
                                                                                                                           gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAACCCGATTCATGc                 >  2:626385/1‑141 (MQ=255)
                                                                                                                           gCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGc                 <  2:129734/141‑1 (MQ=255)
                                                                                                                                 aCAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTAt           >  2:770348/1‑141 (MQ=255)
                                                                                                                                  cAATTAGGCCAAAAAATGGCCTTAAGAATGTTGGAAGTATGAAAGCAGGAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACTTGAACTTTTTATTAAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTAta          <  2:29438/141‑2 (MQ=255)
                                                                                                                                   aaTTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATca         <  2:734713/141‑1 (MQ=255)
                                                                                                                                   aaTTAGGACAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATca         <  2:226485/141‑1 (MQ=255)
                                                                                                                                     ttAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATcaca       <  2:279217/141‑1 (MQ=255)
                                                                                                                                        ggCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTGTTCAAGAGCTAGGGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAAt    >  1:302544/1‑141 (MQ=255)
                                                                                                                                        agCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATGCATGCTGGTATCACAAAt    <  2:553469/140‑1 (MQ=255)
                                                                                                                                          ccAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAGGAGCTAGAGGCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGa  >  1:291344/1‑141 (MQ=255)
                                                                                                                                          ccAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGCTTCATGCTGGTATCACAAATGa  >  1:813167/1‑141 (MQ=255)
                                                                                                                                          ccAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGa  >  1:535388/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AAAAGGCTCTCCTCAAATAAGACGACGCTTTATAGATATGGAGTTGGGCCAAATTTCTGCTGTTTACTTAAATGATTTAGCTCAATACCAACGTATTTTAAAGCAAAAGAATAATTACTTAAAGCAGTTACAATTAGGCCAAAAAAAGGACTTAACAATGTTGGAAGTATTAAATCAGCAGTTTGCTGAATATGCAATGAAAGTAACTGATAAACGTGCACATTTTATTCAAGAGCTAGAGTCGTTAGCTAAACCGATTCATGCTGGTATCACAAATGA  >  USA300TCH1516_ALE/4310‑4588

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: