Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 790,718 T→A 6.7% intergenic (+217/‑257) kipA_1 → / → ltaS KipI antagonist/Lipoteichoic acid synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE790,7180TA 6.7% 419.8 / 16.5 120intergenic (+217/‑257)kipA_1/ltaSKipI antagonist/Lipoteichoic acid synthase
Reads supporting (aligned to +/- strand):  ref base T (58/54);  new base A (3/5);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 126 total aligned reads displayed.

ATAAAGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTT  >  USA300TCH1516_ALE/790628‑790807
                                                                                           |                                                                                         
aTAAAGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAAtttt                                                                                           <  1:192321/91‑1 (MQ=255)
 tAAAGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTa                                                                                          <  1:3774731/91‑1 (MQ=255)
   aaGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTaaa                                                                                        >  1:1918501/1‑91 (MQ=255)
    aGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAAt                                                                                       <  1:1058466/91‑1 (MQ=255)
    aGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAAt                                                                                       >  1:2338612/1‑91 (MQ=255)
    aGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTAAAAt                                                                                       <  1:2939132/91‑1 (MQ=255)
    aGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTAAAAt                                                                                       <  1:346451/91‑1 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    <  1:3433456/91‑1 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    >  1:1075912/1‑91 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    <  1:625449/91‑1 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    <  1:625367/91‑1 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    <  1:625294/91‑1 (MQ=255)
       gtgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTg                                                                                    <  1:3770147/91‑1 (MQ=255)
        tgtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGt                                                                                   >  1:1504790/1‑91 (MQ=255)
         gtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGtt                                                                                  >  1:2729276/1‑91 (MQ=255)
         gtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGtt                                                                                  >  1:2729216/1‑91 (MQ=255)
         gtTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGtt                                                                                  >  1:1519220/1‑91 (MQ=255)
             ttACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTAACAAATTTAAAATTTGtaatc                                                                               <  1:1685283/91‑5 (MQ=38)
              tACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTg                                                                             <  1:1669729/91‑1 (MQ=255)
               aCTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGa                                                                            >  1:1014174/1‑91 (MQ=255)
                cTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGaa                                                                           >  1:2902238/1‑91 (MQ=255)
                cTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGaa                                                                           >  1:4072643/1‑91 (MQ=255)
                 ttGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAAt                                                                          >  1:1654834/1‑91 (MQ=255)
                 ttGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAAt                                                                          <  1:3947273/91‑1 (MQ=255)
                 ttGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAAt                                                                          <  1:2524910/91‑1 (MQ=255)
                  tGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATg                                                                         <  1:373708/91‑1 (MQ=255)
                    aaTTTTGTCTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTATTTGAATGat                                                                       <  1:2080194/91‑1 (MQ=255)
                    aaTTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGat                                                                       <  1:269742/91‑1 (MQ=255)
                    aaTTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGat                                                                       <  1:2035184/91‑1 (MQ=255)
                         tGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTGTTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAatat                                                                  <  1:4139603/91‑1 (MQ=255)
                         tGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAatat                                                                  >  1:3513160/1‑91 (MQ=255)
                          gACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATAtt                                                                 >  1:2679876/1‑91 (MQ=255)
                           aCTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATATGATTACAAATTTTAAATTTGTTAATTGAATGATAATATta                                                                <  1:3588868/91‑1 (MQ=255)
                            cTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATtaa                                                               >  1:1918759/1‑91 (MQ=255)
                              aaaaTTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATtaaat                                                             <  1:37252/91‑1 (MQ=255)
                               aaaTTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATtaaata                                                            >  1:1844327/1‑91 (MQ=255)
                               aaaTTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATtaaata                                                            >  1:797434/1‑91 (MQ=255)
                                 aTTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATtaaataaa                                                          <  1:3819248/91‑1 (MQ=255)
                                    aCTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGaa                                                       >  1:3446930/1‑91 (MQ=255)
                                     ctctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGaaa                                                      <  1:2461938/91‑1 (MQ=255)
                                      tctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAAc                                                     <  1:1603488/91‑1 (MQ=255)
                                      tctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAAc                                                     >  1:2563473/1‑91 (MQ=255)
                                      tctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAAc                                                     >  1:2442938/1‑91 (MQ=255)
                                      tctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAAc                                                     >  1:504583/1‑91 (MQ=255)
                                      tctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAAc                                                     >  1:2259546/1‑91 (MQ=255)
                                       ctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACt                                                    <  1:2323650/91‑1 (MQ=255)
                                       ctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACt                                                    >  1:1894700/1‑91 (MQ=255)
                                       ctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACt                                                    <  1:161707/91‑1 (MQ=255)
                                        tataTTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACtt                                                   >  1:2017233/1‑91 (MQ=255)
                                        tataTTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACtt                                                   >  1:182148/1‑91 (MQ=255)
                                         atatTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTa                                                  <  1:1007685/91‑1 (MQ=255)
                                         atatTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTa                                                  >  1:1613810/1‑91 (MQ=255)
                                         atatTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTa                                                  <  1:1009122/91‑1 (MQ=255)
                                          tatttattAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTac                                                 <  1:2346813/91‑1 (MQ=255)
                                           atttattAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTaca                                                <  1:3707080/91‑1 (MQ=255)
                                             ttattaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTAcaca                                              <  1:804394/91‑1 (MQ=255)
                                              tattaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACaa                                             >  1:1483870/1‑91 (MQ=255)
                                               attaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAg                                            >  1:3634248/1‑91 (MQ=255)
                                               attaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAg                                            >  1:3634222/1‑91 (MQ=255)
                                                ttaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGc                                           <  1:407753/91‑1 (MQ=255)
                                                 taATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCa                                          >  1:341320/1‑91 (MQ=255)
                                                  aaTTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCaa                                         >  1:766589/1‑91 (MQ=255)
                                                  aaTTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCaa                                         >  1:2930578/1‑91 (MQ=255)
                                                  aaTTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCaa                                         >  1:1820901/1‑91 (MQ=255)
                                                   aTTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCaaa                                        >  1:1412708/1‑91 (MQ=255)
                                                   aTTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCaaa                                        <  1:2097812/91‑1 (MQ=255)
                                                      gAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAatat                                     <  1:2639596/91‑1 (MQ=255)
                                                       aGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATg                                    <  1:2768424/91‑1 (MQ=255)
                                                       aGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATg                                    <  1:2768157/91‑1 (MQ=255)
                                                         cTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAg                                  <  1:1661886/91‑1 (MQ=255)
                                                           aTGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGtt                                >  1:1958846/1‑91 (MQ=255)
                                                            tGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTg                               >  1:2568915/1‑91 (MQ=255)
                                                            tGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTg                               <  1:1535840/91‑1 (MQ=255)
                                                              cTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTa                             >  1:308372/1‑91 (MQ=255)
                                                              cTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTa                             >  1:308176/1‑91 (MQ=255)
                                                                tattattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGc                           >  1:2284029/1‑91 (MQ=255)
                                                                tattattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGc                           >  1:2060882/1‑91 (MQ=255)
                                                                  ttattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGccc                         <  1:211109/91‑1 (MQ=255)
                                                                  ttattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGccc                         >  1:1606507/1‑91 (MQ=255)
                                                                   tattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCa                        <  1:2824829/91‑1 (MQ=255)
                                                                    attaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCaa                       <  1:1015524/91‑1 (MQ=255)
                                                                    attaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCaa                       >  1:135405/1‑91 (MQ=255)
                                                                     ttaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCaaa                      <  1:1358741/91‑1 (MQ=255)
                                                                      taCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCaaaa                     <  1:3315084/91‑1 (MQ=255)
                                                                       aCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAAt                    <  1:279361/91‑1 (MQ=255)
                                                                          aTTTGATTACACATTTTAAATTTGTTCAATGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACt                 >  1:2650755/1‑91 (MQ=255)
                                                                            ttGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTg               <  1:3354096/91‑1 (MQ=255)
                                                                               aTTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTa            >  1:3134479/1‑91 (MQ=255)
                                                                               aTTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTa            >  1:1152871/1‑91 (MQ=255)
                                                                                ttACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTaa           >  1:1245454/1‑91 (MQ=255)
                                                                                ttACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTaa           >  1:688302/1‑91 (MQ=255)
                                                                                  aaCAAATTTAAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATCCTTGTTaaa          >  1:3974024/2‑91 (MQ=255)
                                                                                  aaCAAATTTAAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTaaa          >  1:3974085/2‑91 (MQ=255)
                                                                                   aCAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAAt         <  1:1272098/91‑1 (MQ=255)
                                                                                   aCAAATTTAAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCCAATATGAGTTGTAGCCCAAAATACTTGTTAAAt         >  1:3685328/1‑91 (MQ=255)
                                                                                    cAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATc        <  1:2273895/91‑1 (MQ=255)
                                                                                      aaTTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCaa      >  1:474216/1‑91 (MQ=255)
                                                                                      aaTTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCaa      >  1:474233/1‑91 (MQ=255)
                                                                                          ttAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACTCAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGtt  <  1:3297388/91‑1 (MQ=255)
                                                                                          ttAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGtt  >  1:2433026/1‑91 (MQ=255)
                                                                                           |                                                                                         
ATAAAGTGTGTTTTTACTTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTT  >  USA300TCH1516_ALE/790628‑790807

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