Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,615 C→G 6.3% 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 6.3% 814.2 / 28.6 232intergenic (+51/+89)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base C (114/103);  new base G (7/8);  total (121/111)
Fisher's exact test for biased strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

TGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGT  >  USA300TCH1516_ALE/1421525‑1421704
                                                                                           |                                                                                         
tGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                           <  1:2126890/91‑1 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                          >  2:3676975/1‑91 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                          >  1:3733308/1‑91 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                          >  1:2252769/1‑91 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                          >  2:3281496/1‑91 (MQ=255)
 gAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACt                                                                                          >  2:388111/1‑91 (MQ=255)
  aaaaaTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTa                                                                                         <  2:3011920/91‑1 (MQ=255)
      aTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtg                                                                                     >  1:3435496/1‑91 (MQ=255)
       ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                    <  1:1927127/91‑1 (MQ=255)
        tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                                   <  2:2038533/91‑1 (MQ=255)
        tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                                   <  2:3330544/91‑1 (MQ=255)
        tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAAtgtg                                                                                    <  1:2615851/91‑4 (MQ=38)
         aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                  >  1:3016502/1‑91 (MQ=255)
          aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                 <  2:1313774/91‑1 (MQ=255)
           cAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgaga                                                                                >  1:1056877/1‑91 (MQ=255)
             gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTAc                                                                                           <  2:628239/78‑1 (MQ=255)
              aTGAAATTGAGCAAGCAAATAACTAGGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                             <  1:3600656/91‑1 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                            <  2:1483849/91‑1 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                            >  2:2486104/1‑91 (MQ=255)
               tGAAATTGAGCAAGCAAATAACTAAGAGGCGAGACAATTTATAGTAGGACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCt                                                                            <  2:1112022/91‑1 (MQ=21)
                gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                           <  1:1852807/91‑1 (MQ=255)
                gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                           <  2:1206168/91‑1 (MQ=255)
                gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCtt                                                                           <  2:3791712/91‑1 (MQ=21)
                   aTTGAGCCAGCAAATAACTAAGATGCGAGACAATTTATAGAAGAACACTATAAACCAAGATTTTCAACTTACTACTAGTGAGAATCTTGGt                                                                        >  2:3609771/1‑91 (MQ=255)
                   aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                        >  1:1585206/1‑91 (MQ=255)
                   aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                        >  1:1585694/1‑91 (MQ=255)
                   aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGt                                                                        <  1:2001303/91‑1 (MQ=21)
                     tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                      >  1:2708730/1‑91 (MQ=255)
                     tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                      >  1:2915253/1‑91 (MQ=255)
                     tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGttt                                                                      <  2:1230173/91‑1 (MQ=21)
                      gAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                     <  2:1358494/91‑1 (MQ=255)
                       agcaagcaAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttttt                                                                    >  1:2070522/1‑91 (MQ=255)
                        gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                   <  2:742843/91‑1 (MQ=255)
                        gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                   >  2:98825/1‑91 (MQ=255)
                        gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                   <  1:2300568/91‑1 (MQ=255)
                         caagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGAttttttt                                                                  >  2:3447006/1‑91 (MQ=255)
                          aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                 <  2:1716847/91‑1 (MQ=255)
                          aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATAGAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                 <  2:1685493/91‑1 (MQ=255)
                           agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                                >  1:3413047/1‑91 (MQ=255)
                           agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                                >  1:3538580/1‑91 (MQ=255)
                           agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                                >  1:656038/1‑91 (MQ=255)
                           agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                                >  1:606343/1‑91 (MQ=255)
                             caaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAgg                                                              <  1:755986/91‑1 (MQ=255)
                               aaTAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                           <  1:1475352/76‑1 (MQ=255)
                                   aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                        <  2:2854131/91‑1 (MQ=255)
                                   aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                        >  2:265745/1‑91 (MQ=255)
                                   aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                        <  1:3018889/91‑1 (MQ=255)
                                    cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                       <  1:286575/91‑1 (MQ=255)
                                    cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                       >  2:1041243/1‑91 (MQ=255)
                                    cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                       >  2:1252698/1‑91 (MQ=255)
                                      aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                     <  2:7832/91‑1 (MQ=255)
                                      aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                     <  2:5679/91‑1 (MQ=255)
                                       aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                    <  1:3236930/91‑1 (MQ=255)
                                        gATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaata                                                   >  1:1399794/1‑91 (MQ=255)
                                           gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                 <  1:3833080/74‑1 (MQ=255)
                                               gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                            >  1:3695457/1‑91 (MQ=255)
                                               gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                            >  1:2558381/1‑91 (MQ=255)
                                                aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAat                                           >  1:279414/1‑91 (MQ=255)
                                                aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGATTTTTTGAGGGTTTATGAATAATAGAAat                                           >  1:3117928/1‑91 (MQ=255)
                                                 cAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAata                                          >  1:2417407/1‑91 (MQ=255)
                                                  aaTTTATAGTAGAACAGTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                         <  2:3704203/91‑1 (MQ=255)
                                                   aTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatata                                        <  1:444198/91‑1 (MQ=255)
                                                    tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                       <  1:968346/91‑1 (MQ=255)
                                                    tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                       <  1:3007994/91‑1 (MQ=255)
                                                    tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                       <  1:222365/91‑1 (MQ=255)
                                                     ttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatata                                      >  1:1999158/1‑91 (MQ=255)
                                                      tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                                     >  2:809640/1‑91 (MQ=255)
                                                       ataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAtt                                    <  2:1413989/91‑1 (MQ=255)
                                                        tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAttt                                   >  2:2255983/1‑91 (MQ=255)
                                                           tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                                >  1:3480152/1‑91 (MQ=255)
                                                           tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                                >  1:1268465/1‑91 (MQ=255)
                                                           tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                                >  1:1970473/1‑91 (MQ=255)
                                                             gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                              <  2:10569/91‑1 (MQ=255)
                                                               aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                            <  2:2621117/91‑1 (MQ=255)
                                                               aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                            <  1:3181104/91‑1 (MQ=255)
                                                                cAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAg                           <  2:1171572/91‑1 (MQ=255)
                                                                 aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                          <  1:2244618/91‑1 (MQ=255)
                                                                   tataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGa                        <  1:901788/91‑1 (MQ=255)
                                                                   tataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGa                        <  2:3695457/91‑1 (MQ=255)
                                                                     aaaaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtt                      >  1:1367695/2‑91 (MQ=18)
                                                                      aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                     <  2:2177965/91‑1 (MQ=255)
                                                                       aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGCGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                    >  1:3091109/1‑91 (MQ=255)
                                                                       aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                    >  1:2511447/1‑91 (MQ=255)
                                                                       aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                    <  1:2560467/91‑1 (MQ=255)
                                                                       aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                    >  1:770046/1‑91 (MQ=21)
                                                                        aCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTta                   >  2:2087492/1‑91 (MQ=255)
                                                                         ccAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat                  >  1:2676259/1‑91 (MQ=255)
                                                                         ccAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat                  >  1:844317/1‑91 (MQ=21)
                                                                           aaGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg                >  2:2392057/1‑91 (MQ=21)
                                                                             gATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGta              >  2:949956/1‑91 (MQ=255)
                                                                              aTTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaa             >  2:2100352/1‑91 (MQ=255)
                                                                               ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat            <  2:1911140/91‑1 (MQ=255)
                                                                                 ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa          >  1:1686565/1‑91 (MQ=255)
                                                                                   cacaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACt        <  2:3518321/91‑1 (MQ=255)
                                                                                    gcaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGATTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt       <  2:3842895/90‑1 (MQ=255)
                                                                                     caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtta      <  2:2201414/91‑2 (MQ=255)
                                                                                       aTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTa     <  2:1509418/91‑1 (MQ=255)
                                                                                       aTTACTAATTGTGAGAATCTGGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTa     <  1:960717/91‑1 (MQ=255)
                                                                                          aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt  <  2:2977271/91‑1 (MQ=255)
                                                                                          aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt  <  1:1252698/91‑1 (MQ=255)
                                                                                           |                                                                                         
TGAAAAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGT  >  USA300TCH1516_ALE/1421525‑1421704

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