Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,273,186 G→C 5.1% P339P (CCG→CCC)  hslU → ATP‑dependent protease ATPase subunit HslU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,273,1860GC 5.1% 225.8 / 3.2 69P339P (CCG→CCC) hslUATP‑dependent protease ATPase subunit HslU
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (2/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TCAAGATGTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGAAGCATTGC  >  USA300TCH1516_ALE/1273048‑1273316
                                                                                                                                          |                                                                                                                                  
tCAAGATGTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAg                                                                                                                                  <  1:529454/141‑1 (MQ=255)
  aaGATGTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTg                                                                                                                                >  2:315560/1‑141 (MQ=255)
     aTGTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCAGAGAGACt                                                                                                                             >  2:222484/1‑141 (MQ=255)
       gTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTAGCTAAGCCGAGTGACTTg                                                                                                                           <  1:272384/141‑1 (MQ=255)
        tctcAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACGGGGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGa                                                                                                                          <  2:73090/141‑1 (MQ=255)
         ctcAAGACAAGGTGTTCAAAGAGATATTTTCCAAATACTGGAAGGTAGCGTTAGTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGAt                                                                                                                         >  1:500047/1‑141 (MQ=255)
         ctcAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTTAGCCGAGTGACTTGAt                                                                                                                         <  1:278942/141‑1 (MQ=255)
         ctcAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTTAGCCGAGTGACTTGAt                                                                                                                         <  1:313582/141‑1 (MQ=255)
         ctcAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTCAGCCGAGTGACTTGAt                                                                                                                         <  1:278893/141‑1 (MQ=255)
              gACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAg                                                                                                                    <  1:17466/141‑1 (MQ=255)
                 aaGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTTGACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTACGACGGGTGAGTTGGTACAAGAAt                                                                                                                 >  1:95106/1‑141 (MQ=255)
                 aaGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAAt                                                                                                                 >  1:468142/1‑141 (MQ=255)
                 aaGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAAt                                                                                                                 >  1:91740/1‑141 (MQ=255)
                  aGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGTTACCAGAAtt                                                                                                                >  2:559878/1‑141 (MQ=255)
                   ggTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTg                                                                                                               >  2:99582/1‑141 (MQ=255)
                   ggTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGGTATTCAGACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCCAGTGACTTGATACCAGAAttt                                                                                                               >  2:364292/1‑140 (MQ=255)
                   ggTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGGTATTCAGACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCCAGTGACTTGATACCAGAATTg                                                                                                               >  2:364302/1‑141 (MQ=255)
                   ggTGATCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGGATTg                                                                                                               >  1:381396/1‑141 (MQ=255)
                      gtTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGCTACCAGAATTGCaa                                                                                                            >  1:230586/1‑141 (MQ=255)
                        tCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAgg                                                                                                          <  1:476287/141‑1 (MQ=255)
                         cAACGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGt                                                                                                         <  1:430594/141‑1 (MQ=255)
                            agagATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACTAGAATTGCAAGGTCGt                                                                                                      >  1:427582/1‑141 (MQ=255)
                             gagaTATTTTACCTACACTTGAAGGTAGCGTTATTCATACCAACTATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTTCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGGCGtt                                                                                                     <  1:395357/141‑1 (MQ=255)
                                atatTTTACCTATACTTGACGGTAGCGTTATTCAAACCAAATATGGTACTGTCAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGCTACCAGAATTGCAAGGTCGTTTTc                                                                                                  <  1:358725/141‑1 (MQ=255)
                                atacTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTc                                                                                                  <  1:2401/141‑1 (MQ=255)
                                 tatTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTcc                                                                                                 >  1:487266/1‑141 (MQ=255)
                                   ttttACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGa                                                                                               >  2:141799/1‑141 (MQ=255)
                                          tataCTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGTGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGt                                                                                        <  2:276815/141‑1 (MQ=255)
                                          cataCTTGAAGGTAGCTTTATTCAAACCAAATATCGTACTGTGAATACTGAACATATGCTGTTTATCGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGt                                                                                        <  1:222484/140‑1 (MQ=255)
                                               ttGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAAc                                                                                   <  2:119445/141‑1 (MQ=255)
                                               ttGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAAc                                                                                   <  1:182828/141‑1 (MQ=255)
                                                  aaGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTg                                                                                >  1:424312/1‑141 (MQ=255)
                                                  aaGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTg                                                                                >  1:30195/1‑141 (MQ=255)
                                                    ggTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGAt                                                                              >  1:211008/1‑141 (MQ=255)
                                                       aGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGt                                                                           >  2:321661/1‑141 (MQ=255)
                                                        gCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGtt                                                                          <  2:492682/141‑1 (MQ=255)
                                                                cAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTa                                                                  >  1:144261/1‑141 (MQ=255)
                                                                  aaCCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTgga                                                                >  1:262761/1‑141 (MQ=255)
                                                                     cAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTagaaga                                                             <  1:315560/141‑1 (MQ=255)
                                                                       aaTATGGTACTGTGAAAACTGAACATAGGCTGTTTAAAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACCTGATACCAGAATTGCAAGGTCGGTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGAtt                                                           <  1:111996/141‑1 (MQ=255)
                                                                           tGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGt                                                       <  2:70791/141‑1 (MQ=255)
                                                                              tACTGTGAATGCCGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGTGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAg                                                    <  1:6009/141‑1 (MQ=255)
                                                                                cTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGaa                                                  >  1:237818/1‑141 (MQ=255)
                                                                                  gtgAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAAtt                                                >  1:147058/1‑141 (MQ=255)
                                                                                             aCATATGCTGTTTCTAGGTGCTGGCGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGCTTAGAGTTGACCTGTATCGTTTATCGGTCGAAGATTTTGTAAGAATTTTGACAGAAcc                                     >  2:161281/1‑141 (MQ=255)
                                                                                              cATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGCAGATTTTGTAAGAATTTTGACAGAACCa                                    <  1:250062/141‑1 (MQ=255)
                                                                                              cATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTTTAAGAATTTGGACAGAACCa                                    >  1:17229/1‑141 (MQ=255)
                                                                                               atatGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCaa                                   <  1:141799/141‑1 (MQ=255)
                                                                                                tatGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCaaa                                  >  1:465479/1‑141 (MQ=255)
                                                                                                   gCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCGCATTTCAGCTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAAtt                               >  1:408557/1‑141 (MQ=255)
                                                                                                   gCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACGGAACCGAAAtt                               >  2:185048/1‑141 (MQ=255)
                                                                                                   gCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAAtt                               <  2:164358/141‑1 (MQ=255)
                                                                                                     tGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGt                             >  2:524089/1‑141 (MQ=255)
                                                                                                      gTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTc                            <  1:220964/141‑1 (MQ=255)
                                                                                                       tttATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCa                           <  1:227334/141‑1 (MQ=255)
                                                                                                       tttATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCa                           >  1:500785/1‑141 (MQ=255)
                                                                                                         tataGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCatt                         <  2:211008/141‑1 (MQ=255)
                                                                                                            aGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGCTACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCattaat                      >  1:242776/1‑141 (MQ=255)
                                                                                                             ggTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCattaatt                     >  1:483308/1‑141 (MQ=255)
                                                                                                                        tttCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTCATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGa          >  1:379477/1‑141 (MQ=255)
                                                                                                                         ttCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTACTTAAACAATATGaa         >  1:70462/1‑141 (MQ=255)
                                                                                                                         ttCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGaa         >  1:5125/1‑141 (MQ=255)
                                                                                                                         ttCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGaa         >  2:186079/1‑141 (MQ=255)
                                                                                                                          tCCATGTATCTAAGCCCAGTGACTTGATACCAGAATTGTAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAATATTTTGTAAGAATTTTGACAGAACCAAAATTGTCGTTAATTAAACAATATGAAg        <  1:364302/141‑1 (MQ=255)
                                                                                                                          tCCATGTATCTAAGCCCAGTGACTTGATACCAGAATTGTAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAAATTTTGTAAGAATTTTGACAGAACCAAAATTGTATTTAATTAAACAATATGAAg        <  1:364292/141‑1 (MQ=255)
                                                                                                                             atttATCTAAGCCGATTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGAAGCAt     <  2:17163/138‑1 (MQ=255)
                                                                                                                             aTGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAAGTGTCATTAATTAAACAATATGAAGCAt     <  2:144261/141‑1 (MQ=255)
                                                                                                                                tATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGAAGCATTGc  >  2:223503/1‑141 (MQ=255)
                                                                                                                                tATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGAAGCATTGc  >  1:102084/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                  
TCAAGATGTCTCAAGACAAGGTGTTCAAAGAGATATTTTACCTATACTTGAAGGTAGCGTTATTCAAACCAAATATGGTACTGTGAATACTGAACATATGCTGTTTATAGGTGCTGGAGCTTTCCATGTATCTAAGCCGAGTGACTTGATACCAGAATTGCAAGGTCGTTTTCCGATTAGAGTTGAACTTGATAGTTTATCGGTAGAAGATTTTGTAAGAATTTTGACAGAACCAAAATTGTCATTAATTAAACAATATGAAGCATTGC  >  USA300TCH1516_ALE/1273048‑1273316

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