Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 683,202 T→A 6.0% Y23F (TAT→TTT)  ftsL_1 ← Cell division protein FtsL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE683,2020TA 6.0% 288.3 / 8.6 101Y23F (TAT→TTT) ftsL_1Cell division protein FtsL
Reads supporting (aligned to +/- strand):  ref base T (44/50);  new base A (3/3);  total (48/53)
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.01e-01

Only 100 of 102 total aligned reads displayed.

CATCAGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACTTTAAAA  >  USA300TCH1516_ALE/683065‑683342
                                                                                                                                         |                                                                                                                                              
catcaGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCGAATCTTCAATTTGATTTTGTTTCATTGTAACTTGTTGTTGAGGAATCTGATatc                                                                                                                                             >  1:85260/1‑139 (MQ=255)
catcaGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGGTacc                                                                                                                                             >  1:329926/1‑139 (MQ=255)
  tcaGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCa                                                                                                                                           <  2:206542/141‑1 (MQ=255)
  tcaGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCa                                                                                                                                           <  1:672885/141‑1 (MQ=255)
     taaaTAACTGTTTTTTATTGTCTTTTTACTAATTTATTTATCTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAAtt                                                                                                                                        <  2:452194/140‑1 (MQ=255)
     gAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAAtt                                                                                                                                        <  1:594850/141‑1 (MQ=255)
     gAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAAtt                                                                                                                                        <  1:787440/141‑1 (MQ=255)
     gAAATAACTGTTTTTGATTGTCTCTTTACTAAATCATGTAACCCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAAtt                                                                                                                                        <  2:504934/141‑1 (MQ=255)
     aaaaTAACTGTTTTTGATTGTACCTTTACTTATTACAGTATCTCTCCATTTTTATCTGTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACGTGTTGTTTAGGAATCTGATAACCAAtt                                                                                                                                        <  2:724282/140‑1 (MQ=255)
       aaTAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTaa                                                                                                                                      >  2:640719/1‑141 (MQ=255)
          aaCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTTATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTaaaca                                                                                                                                   <  1:411232/141‑3 (MQ=255)
          aaCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAAc                                                                                                                                   <  2:439315/141‑1 (MQ=255)
          aaCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAAc                                                                                                                                   >  1:817286/1‑141 (MQ=255)
           aCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACa                                                                                                                                  >  1:221197/1‑141 (MQ=255)
           aCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACa                                                                                                                                  >  1:663894/1‑141 (MQ=255)
           aCTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTAATTTTGTTTCATTTTAACTTGTTTTTTAGTAATCTGATAACCTATTAAAACa                                                                                                                                  >  2:415456/1‑141 (MQ=255)
           aCTGTTTTTGATTGTCTCTTTACTAATTTATGTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTCGTTTAGGAATCTGATAACCAATTAAAACa                                                                                                                                  >  2:863524/1‑141 (MQ=255)
            cTGTTTTTGATTGTCTCTTTACTAATTTATTTAGCTCTGCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAg                                                                                                                                 <  2:26357/141‑1 (MQ=255)
             tGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGc                                                                                                                                <  2:261877/141‑1 (MQ=255)
             tGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAAGAATCGTATAACCAATTAAAACatc                                                                                                                                <  2:473404/141‑3 (MQ=255)
                   tGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACaaaa                                                                                                                          <  2:143320/141‑1 (MQ=255)
                      ttGTCTCTTTACTAATTTATTTTACTATCCTTTTTTATCTCTTAATATACATTGTTCTTCCTGCAAATCTTCAATTTGATTTTCTTTCATTTTAATTTGTTGTTTAGGAATCTGATAACCAACTAAAACAGCACAAact‑aa                                                                                                                      >  2:127811/1‑141 (MQ=255)
                                  tAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTTTTATAAATTTATATGTTCAATTTGATTTTGTTACAATTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTAc                                                                                                          <  2:114880/141‑1 (MQ=255)
                                  tAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTAc                                                                                                          >  2:788051/1‑141 (MQ=255)
                                  tAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTAc                                                                                                          <  2:608346/141‑1 (MQ=255)
                                  tAATTTATTTAACTCTCCATATTTATCTCTTAATATACGTCGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTATTTAGAAATTTAATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTAc                                                                                                          >  2:519898/1‑141 (MQ=255)
                                   tatttatttAACTCTCCATTTTTATCTCTTAATAGACGTTGTTCTTCCTGCAAATCTGCAATTTGATTTTGTTTCATTTTAACTGGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAGAT‑AATATTACTGTTACa                                                                                                         <  1:269446/140‑1 (MQ=255)
                                     tttatttaACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAag                                                                                                       >  2:438073/1‑141 (MQ=255)
                                     tttatttaACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAag                                                                                                       >  2:1027847/1‑141 (MQ=255)
                                     tttatttaACTCTCCATTTTTATCTCCTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAag                                                                                                       >  2:493292/1‑141 (MQ=255)
                                       tatttaACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGTAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAagag                                                                                                     <  2:373548/141‑1 (MQ=255)
                                         tttaACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTa                                                                                                   >  1:252641/1‑141 (MQ=255)
                                         tttaACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTa                                                                                                   >  1:579367/1‑141 (MQ=255)
                                           taACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTGATGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTaaa                                                                                                 >  1:315682/1‑141 (MQ=255)
                                           taACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTaaa                                                                                                 <  2:957554/141‑1 (MQ=255)
                                           taACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTaaa                                                                                                 >  1:665291/1‑141 (MQ=255)
                                             aCTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACa                                                                                               >  2:961034/1‑141 (MQ=255)
                                                 tcCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACgtaat                                                                                           >  1:216669/1‑141 (MQ=255)
                                                 tcCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACataat                                                                                           <  1:1006546/141‑1 (MQ=255)
                                                  ccATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGACTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACataata                                                                                          <  1:715032/141‑1 (MQ=255)
                                                   cATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACAtaataa                                                                                         <  2:120775/141‑1 (MQ=255)
                                                   cATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACAtaataa                                                                                         >  1:236173/1‑141 (MQ=255)
                                                     tttttATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑tt                                                                                      <  1:733309/141‑1 (MQ=255)
                                                     tttttATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑tt                                                                                      <  1:640719/141‑1 (MQ=255)
                                                      ttttATCTCTTAATATACGTTGTTCTTCCTGCAAATATTCAATTTGGTTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAACAT‑AATCTTACGGGTACAAGAGTAAACATAATAA‑ttt                                                                                     <  2:18780/141‑1 (MQ=255)
                                                        ttATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑ttttt                                                                                   <  1:160811/141‑1 (MQ=255)
                                                        ttATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑ttttt                                                                                   <  1:840123/141‑1 (MQ=255)
                                                        ttATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAAATTATAATTTATGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑ttttt                                                                                   <  2:624144/141‑1 (MQ=255)
                                                         tATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑tatttt                                                                                  >  1:906416/1‑141 (MQ=255)
                                                          aTCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTATAACAGAACAAAAT‑AATATTACTGTTACACGAGTAAACATAATAA‑ttttttt                                                                                 >  2:154934/1‑140 (MQ=255)
                                                          aTCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑ACTATTACTGTTACAAGAGTCCACATAAAAG‑TTTTTTg                                                                                 >  2:613288/1‑141 (MQ=255)
                                                          aTCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTg                                                                                 <  1:264827/141‑1 (MQ=255)
                                                          aTCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTATTATTTAGGAATATGATAATCAATTAAAACAGCACAAAAT‑AATATTACTTTTACAAAATTAAACAACATAA‑tatttta                                                                                 >  2:489075/1‑140 (MQ=255)
                                                           tctcTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACTAGAGTAAACATAATAA‑tttttttc                                                                                >  2:853510/1‑139 (MQ=255)
                                                           tctcTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGc                                                                                >  2:538638/1‑141 (MQ=255)
                                                             tctTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCtt                                                                              >  1:403080/1‑141 (MQ=255)
                                                             tctTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCtt                                                                              >  1:43967/1‑141 (MQ=255)
                                                             tctTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCtt                                                                              >  1:45160/1‑141 (MQ=255)
                                                              tttAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTc                                                                             <  2:765758/140‑1 (MQ=255)
                                                              tttAATATACGTTGTTCTTACTTGAAATATTAAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTc                                                                             <  2:114121/140‑1 (MQ=255)
                                                                  atatACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCattt                                                                         <  2:465529/141‑1 (MQ=255)
                                                                     taCGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCatttatt                                                                      >  1:629433/1‑141 (MQ=255)
                                                                         ttgttCTTCTTGCATAACAACAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCa                                                                  <  2:249347/141‑1 (MQ=255)
                                                                          tgttCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCAc                                                                 >  2:127621/1‑141 (MQ=255)
                                                                          tgttCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATACGA‑AACAAAATTGATCAACAGTATTTCTATTTTTCTGTTACAAGAGTAACCATAATAA‑TTTTTTGCTTCATTTATTTTCAc                                                                 >  2:221223/1‑141 (MQ=255)
                                                                            ttcttcCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTc                                                               >  1:656778/1‑141 (MQ=255)
                                                                              cttcCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTTTTTTTTTTTTTTTCactact                                                             >  2:837063/1‑141 (MQ=255)
                                                                              cttcCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCt                                                             >  1:516059/1‑141 (MQ=255)
                                                                               ttcCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCtt                                                            >  1:171637/1‑141 (MQ=255)
                                                                                    gCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTCCTTCATTTATTTTCACTCCTTATCAt                                                       >  1:80738/1‑141 (MQ=255)
                                                                                     cAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATg                                                      >  1:4043/1‑141 (MQ=255)
                                                                                        aTCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAAt                                                   <  2:817286/141‑1 (MQ=255)
                                                                                        aTCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAAt                                                   <  1:830865/141‑1 (MQ=255)
                                                                                        aTCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAAt                                                   >  2:579908/1‑141 (MQ=255)
                                                                                          cTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCTTTTATTTTCACTACTTATCATGATTaa                                                 >  2:302871/1‑141 (MQ=255)
                                                                                             cAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAAc                                              <  2:752770/141‑1 (MQ=255)
                                                                                             cAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAAc                                              <  1:841167/141‑1 (MQ=255)
                                                                                                    aTTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATcc                                       <  2:49662/141‑1 (MQ=255)
                                                                                                     ttttGTTTCATTTTAACTTGTTGTTTAGGAATCTGATATCCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCg                                      <  1:438073/141‑1 (MQ=255)
                                                                                                     ttttGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCg                                      <  1:493292/141‑1 (MQ=255)
                                                                                                     ttttGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCg                                      <  1:1027847/141‑1 (MQ=255)
                                                                                                         gTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCgtcag                                  >  1:614804/1‑141 (MQ=255)
                                                                                                                     cTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATg                      <  1:579908/141‑1 (MQ=255)
                                                                                                                       tgttgtTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCGCTCCTTATCCTGAATAAAACTATATCCGTCAGTCATTTTCTATGTa                    >  1:538619/1‑141 (MQ=255)
                                                                                                                       tgttgtTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTa                    >  2:829723/1‑141 (MQ=255)
                                                                                                                           gttTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGtt                >  1:80383/1‑141 (MQ=255)
                                                                                                                           gttTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGtt                >  1:59236/1‑141 (MQ=255)
                                                                                                                              tAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCAt             <  1:106996/141‑1 (MQ=255)
                                                                                                                                tgAATCTAATAACCAATTAAAACAGTAAAAAAT‑AATATTACTTTTACAAGATTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTa           <  2:693531/140‑1 (MQ=255)
                                                                                                                                ataaTCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTa           <  2:664337/139‑1 (MQ=255)
                                                                                                                                 gAATCTGATAACCAATTAAAAAAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAATTTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTa           <  2:243661/141‑1 (MQ=255)
                                                                                                                                   aTCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACt        <  1:872435/141‑1 (MQ=255)
                                                                                                                                    tCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACtt       >  1:350290/1‑141 (MQ=255)
                                                                                                                                    tCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACtt       <  2:236173/141‑1 (MQ=255)
                                                                                                                                    tCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACtt       <  1:1011837/141‑1 (MQ=255)
                                                                                                                                     cTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACttt      <  2:45160/141‑1 (MQ=255)
                                                                                                                                     cTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACttt      <  2:43967/141‑1 (MQ=255)
                                                                                                                                     cTGAAAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACttt      <  2:403080/141‑1 (MQ=255)
                                                                                                                                        aTAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTAATTTTAAATGTAAGTTCATTAAATTTaaa   >  1:346998/1‑141 (MQ=255)
                                                                                                                                         tAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTCCATAGAATAA‑GTGTTTGCTTCATTTATTGTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACTTTaaaa  >  1:287196/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CATCAGAAATAACTGTTTTTGATTGTCTCTTTACTAATTTATTTAACTCTCCATTTTTATCTCTTAATATACGTTGTTCTTCCTGCAAATCTTCAATTTGATTTTGTTTCATTTTAACTTGTTGTTTAGGAATCTGATAACCAATTAAAACAGCACAAAAT‑AATATTACTGTTACAAGAGTAAACATAATAA‑TTTTTTGCTTCATTTATTTTCACTCCTTATCATGAATAAAACTATATCCGTCAGTCATTTTCAATGTACGTTCATTAACTTTAAAA  >  USA300TCH1516_ALE/683065‑683342

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: ‑