Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,615 C→G 7.7% intergenic (+51/+89) pepF1_2 → / ← phoU Oligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,421,6150CG 7.7% 894.7 / 48.8 262intergenic (+51/+89)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base C (128/112);  new base G (15/5);  total (144/118)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 266 total aligned reads displayed.

TGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTT  >  USA300TCH1516_ALE/1421525‑1421705
                                                                                           |                                                                                          
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                           >  2:1573633/1‑91 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                           >  1:1461984/1‑91 (MQ=255)
   aaaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTaa                                                                                         >  2:3594511/1‑91 (MQ=255)
    aaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAAt                                                                                        >  1:1495062/1‑91 (MQ=255)
        tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                                    <  1:3901370/91‑1 (MQ=255)
         aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                   >  2:2871615/1‑91 (MQ=255)
         aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                   >  1:1740393/1‑91 (MQ=255)
          tcAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                  <  2:3927884/90‑1 (MQ=255)
          aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                  >  2:1339456/1‑91 (MQ=255)
          aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                  >  2:1849303/1‑91 (MQ=255)
          aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                  >  2:2016187/1‑91 (MQ=255)
            aGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGaa                                                                                >  2:190714/1‑91 (MQ=255)
             gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                               <  2:657978/91‑1 (MQ=255)
             gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                               >  1:1704375/1‑91 (MQ=255)
              aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                              <  2:83355/91‑1 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                             <  1:2117145/91‑1 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                             <  1:457169/91‑1 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCt                                                                             <  2:1745933/91‑1 (MQ=21)
               tGAAATTAAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATACTCACAATTAGTAATGTGAAAATCt                                                                             <  2:643981/91‑1 (MQ=21)
                  aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATAACTAATTGTGAGAATCTTgg                                                                          <  2:2782826/91‑1 (MQ=255)
                     tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                       >  2:3463679/1‑91 (MQ=255)
                        gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                    <  1:3432821/91‑1 (MQ=255)
                         caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                                   >  2:3564908/1‑91 (MQ=255)
                         caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttttt                                                                   <  2:488310/91‑1 (MQ=255)
                          aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                  <  1:1098936/91‑1 (MQ=255)
                          aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                  <  2:2146947/91‑1 (MQ=255)
                               aaTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATCTTCACATTACTAATTGTGAGAATCTTGGTTTATTTGAGGGt                                                             >  2:1171217/1‑91 (MQ=255)
                                  aaCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTa                                                          <  2:2426724/91‑1 (MQ=255)
                                   aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                         >  1:1185240/1‑91 (MQ=255)
                                   aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTGTGAGGGTTTAt                                                         >  1:3585281/1‑91 (MQ=255)
                                     tAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGa                                                       >  1:670385/1‑91 (MQ=255)
                                       aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                     >  2:3809068/1‑91 (MQ=255)
                                         aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                                   <  1:673997/91‑1 (MQ=255)
                                          tGCGAGACAATTTATAGTAGAACAATATAAACTAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataat                                                  <  1:3370748/91‑1 (MQ=255)
                                              agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                              <  2:136560/91‑1 (MQ=255)
                                              agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                              >  1:2577016/1‑91 (MQ=255)
                                              agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                              >  1:2154306/1‑91 (MQ=255)
                                              agaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaa                                              <  2:135970/91‑1 (MQ=255)
                                               gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                             <  1:3204035/91‑1 (MQ=255)
                                                  aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                          <  1:2252011/91‑1 (MQ=255)
                                                   aTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatata                                         >  1:1328877/1‑91 (MQ=255)
                                                     ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                                       >  1:50210/1‑91 (MQ=255)
                                                     ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                                       <  1:1535582/91‑1 (MQ=255)
                                                      tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTACGAATAATAGAAatatatat                                      >  2:109623/1‑91 (MQ=255)
                                                           tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCCTGGTTTTTTTGTGGGTTTTTGAATAATAGAAATAGATATTTATc                                 >  2:2543210/1‑91 (MQ=255)
                                                             gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                               <  1:3709743/91‑1 (MQ=255)
                                                              aaCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTa                              <  1:354487/91‑1 (MQ=255)
                                                               aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                             <  2:2267886/91‑1 (MQ=255)
                                                               aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                             <  1:2695533/91‑1 (MQ=255)
                                                                 aaTATAAACCAAGGTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                           <  2:2381812/91‑1 (MQ=255)
                                                                 aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                           <  2:310216/91‑1 (MQ=255)
                                                                 aaTATAAACCAAGATGTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                           <  1:1283752/91‑1 (MQ=255)
                                                                   tataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGa                         >  1:376951/1‑91 (MQ=255)
                                                                    ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTTGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                        >  2:2085210/1‑91 (MQ=255)
                                                                    ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATCCTTATCGTAAGGGAt                        >  2:646929/1‑91 (MQ=255)
                                                                    ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAt                        >  2:901977/1‑91 (MQ=255)
                                                                      aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                      >  2:732346/1‑91 (MQ=255)
                                                                      aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                      <  2:3738528/91‑1 (MQ=255)
                                                                      aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                      <  1:3620817/91‑1 (MQ=255)
                                                                      aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                      <  2:2122839/91‑1 (MQ=255)
                                                                       aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                     >  2:1781863/1‑91 (MQ=255)
                                                                       aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGACATATATATTTATCGTAAGGGATTTta                    >  1:2871837/1‑91 (MQ=18)
                                                                       aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                    >  1:1002238/1‑91 (MQ=18)
                                                                        aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                    >  2:3136506/1‑91 (MQ=255)
                                                                        aCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                    >  1:3258251/1‑91 (MQ=21)
                                                                         ccAAGATTTTCACATTACTAATTGTTAGAATCTTGGGTTTTTTGAGGGTTTATGAATAAGAGAAATATATATTTATCGTGAGGGATTTtat                   <  2:220172/91‑1 (MQ=255)
                                                                         ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat                   >  1:2513427/1‑91 (MQ=21)
                                                                         ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat                   >  1:1912027/1‑91 (MQ=21)
                                                                          cAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtata                  >  1:964629/1‑91 (MQ=21)
                                                                           aaGATTCTCACAATTAGTAATGTGAAAATCTTGGTTCTTTTGAGGGTTTATGAATAATAGAAATCTATATTTCTCGTAAGGGATTTTATAg                 >  1:1596596/1‑91 (MQ=18)
                                                                           aaGATTCTCACAATTAGTAATGTGAAAATCTTGGTCTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg                 >  2:2575198/1‑91 (MQ=21)
                                                                            aGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                <  2:3300074/91‑1 (MQ=255)
                                                                            aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGCGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                >  2:3003447/1‑91 (MQ=21)
                                                                            aGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                >  1:617642/1‑91 (MQ=21)
                                                                              aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa              <  1:1847301/91‑1 (MQ=255)
                                                                              aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa              <  1:2390139/91‑1 (MQ=255)
                                                                               ttttCACTTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             <  2:3510081/91‑1 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATGTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             <  2:2661886/91‑1 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGACTAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             <  2:3001485/91‑1 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             >  1:731550/1‑91 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             <  1:1515054/91‑1 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             >  1:3000722/1‑91 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat             >  1:3400666/1‑91 (MQ=255)
                                                                                tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata            >  2:3402215/1‑91 (MQ=255)
                                                                                 ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa           >  2:3395743/1‑91 (MQ=255)
                                                                                 ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa           >  1:3165704/1‑91 (MQ=255)
                                                                                  tCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAAc          <  1:650010/91‑1 (MQ=255)
                                                                                   cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt         >  1:2171524/1‑91 (MQ=255)
                                                                                   cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt         >  2:2822224/1‑91 (MQ=255)
                                                                                   cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt         >  1:1051044/1‑91 (MQ=255)
                                                                                    acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGATTAATAGTAATATATATTTATCGTAAGGGATTTTATAGTAATATCtt        >  2:2484256/1‑91 (MQ=255)
                                                                                    acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt        >  1:2344572/1‑91 (MQ=255)
                                                                                    acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt        >  1:1251581/1‑91 (MQ=255)
                                                                                     caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACttt       <  1:2928046/91‑1 (MQ=255)
                                                                                       aTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTa      >  1:2895441/1‑91 (MQ=255)
                                                                                       aTTACTAATTGTGAGAATCTTGGCTTTTTTGAGGGTTTATGAACAATAGATATATATATTTATCGTAAGGGATTTTATAGTAATAACttaa      <  2:3674579/91‑3 (MQ=255)
                                                                                         tACTAATTGTGAGAAGCTTGGTTTTGGTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAg    <  2:3422998/91‑1 (MQ=255)
                                                                                          aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt   >  1:775150/1‑91 (MQ=255)
                                                                                          aCTAATTGTGAGAATCTTGGTTCTTGTGAGGGTTTATGAAGAATAGAAATATATATTTAGCGTAAGGGATTTTATAGTAATAACTTAAAGt   <  2:3190155/91‑1 (MQ=255)
                                                                                           cTAATTGTGAGAATCTTTGTTTTTTTGAGGGTTTATGAATGATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGtt  <  2:1473777/91‑1 (MQ=255)
                                                                                           |                                                                                          
TGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTT  >  USA300TCH1516_ALE/1421525‑1421705

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