Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,619 Δ1 bp 5.7% intergenic (+55/+85) 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,6190T. 5.7% 777.3 / 34.9 229intergenic (+55/+85)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base T (109/107);  new base . (7/6);  total (116/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 233 total aligned reads displayed.

AAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTG  >  USA300TCH1516_ALE/1421529‑1421709
                                                                                           |                                                                                          
  aTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtg                                                                                          >  1:3435496/1‑91 (MQ=255)
  aTTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtg                                                                                          >  1:2292093/1‑91 (MQ=255)
   ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                         >  1:1242757/1‑91 (MQ=255)
   ttAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgt                                                                                         <  2:2443959/91‑1 (MQ=255)
    tAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATtgtg                                                                                        <  2:3330748/91‑1 (MQ=255)
     aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                       >  1:3016727/1‑91 (MQ=255)
     aaCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTga                                                                                       >  1:254383/1‑91 (MQ=255)
      aCAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTgag                                                                                      <  2:1313774/91‑1 (MQ=255)
         gATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAAt                                                                                   >  1:1261902/1‑91 (MQ=255)
          aTGAAATTGGGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                                  <  1:3396927/91‑1 (MQ=255)
          aTGAAATTGAGCAAGCAAATAACTAGGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                                  <  1:3600656/91‑1 (MQ=255)
          aTGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATc                                                                                  <  1:3799758/91‑1 (MQ=255)
           tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                                 <  2:1483849/91‑1 (MQ=255)
           tGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCt                                                                                 >  2:1426929/1‑91 (MQ=255)
            gAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCtt                                                                                <  2:1206168/91‑1 (MQ=255)
             aaaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTg                                                                               >  2:1181162/1‑91 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              <  2:2560412/91‑1 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              >  2:1162201/1‑91 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              <  2:1569213/91‑1 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              <  2:2021413/91‑1 (MQ=255)
              aaTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTgg                                                                              <  2:426081/91‑1 (MQ=255)
               aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                             <  1:1392549/91‑1 (MQ=255)
               aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGt                                                                             >  1:1585678/1‑91 (MQ=255)
               aTTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGt                                                                             <  1:2001303/91‑1 (MQ=21)
                ttGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtt                                                                            <  2:1059870/91‑1 (MQ=255)
                 tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGttt                                                                           >  1:2708730/1‑91 (MQ=255)
                 tGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGttt                                                                           <  1:2764566/91‑1 (MQ=18)
                  gAGCAAGCAAATTACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                          >  1:2974723/1‑91 (MQ=255)
                  gAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttt                                                                          >  1:961022/1‑91 (MQ=255)
                    gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                        <  2:742843/91‑1 (MQ=255)
                    gcaagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGtttttt                                                                        >  2:98825/1‑91 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  2:2654077/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  2:1716847/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  1:425936/91‑1 (MQ=255)
                      aagcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATAGAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTg                                                                      <  2:1685493/91‑1 (MQ=255)
                       agcaaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGa                                                                     >  1:656038/1‑91 (MQ=255)
                         caaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAgg                                                                   >  2:5843/1‑91 (MQ=255)
                         caaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAgg                                                                   >  2:2046154/1‑91 (MQ=255)
                         caaATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAgg                                                                   <  1:755986/91‑1 (MQ=255)
                             tAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGttt                                                               >  1:481916/1‑91 (MQ=255)
                               aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                             >  1:1319732/1‑91 (MQ=255)
                               aCTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTAt                                                             <  2:2856201/91‑1 (MQ=255)
                                cTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATg                                                            <  1:286554/91‑1 (MQ=255)
                                  aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                          <  1:2650486/91‑1 (MQ=255)
                                  aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                          <  2:7832/91‑1 (MQ=255)
                                  aaGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaa                                                          <  1:2911124/91‑1 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGTTTTATGaat                                                         <  1:1955864/91‑1 (MQ=255)
                                   aGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaat                                                         <  1:3236930/91‑1 (MQ=255)
                                    gATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaata                                                        >  1:1567740/1‑91 (MQ=255)
                                     aTGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGaataa                                                       >  2:1052695/1‑91 (MQ=255)
                                       gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                                     >  1:2599556/1‑91 (MQ=255)
                                       gCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAataata                                                     <  1:1450448/91‑1 (MQ=255)
                                           gaCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGaaa                                                 >  1:3695457/1‑91 (MQ=255)
                                            aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAat                                                >  2:918527/1‑91 (MQ=255)
                                            aCAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAat                                                <  1:352016/91‑1 (MQ=255)
                                             cAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAata                                               >  1:2417407/1‑91 (MQ=255)
                                              aaTTTATAGTAGAACAGTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                              <  2:3704203/91‑1 (MQ=255)
                                              aaTTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatat                                              >  2:3640783/1‑91 (MQ=255)
                                                tttATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                            <  1:222365/91‑1 (MQ=255)
                                                tttATAGTAGAACAATATAAACCAAGATTTTCACATCACTAATTGTGAGAATCATGGTTTTTTTGAGGGTTTATGAATAATAGAAatatat                                            >  2:1250310/1‑91 (MQ=255)
                                                  tataGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAatatatat                                          <  2:299389/91‑1 (MQ=255)
                                                   ataGTAGAACAATATAAACCAAGATTTTCACATAACTGATTGTGAGCATCTTGGTTTTTCTGAGGGTTTATGAATAATAGAAATATATAtt                                         >  1:3095873/1‑91 (MQ=255)
                                                    tagtagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATAttt                                        >  2:2130272/1‑91 (MQ=255)
                                                       tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                                     >  1:1970473/1‑91 (MQ=255)
                                                       tagAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATc                                     >  1:3480152/1‑91 (MQ=255)
                                                         gAACAATGTAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                                   <  2:3456622/91‑1 (MQ=255)
                                                         gAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGt                                   <  2:1381486/91‑1 (MQ=255)
                                                           aCAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTaa                                 <  1:3181104/91‑1 (MQ=255)
                                                             aaTATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAgg                               <  1:2244618/91‑1 (MQ=255)
                                                              atatAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAggg                              >  1:1911140/1‑91 (MQ=255)
                                                               tataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGa                             <  1:2109694/91‑1 (MQ=255)
                                                                ataAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTATATGAATAATAGAAATATATATTTATCGTAAGGGAt                            >  2:2691066/1‑91 (MQ=255)
                                                                  aaaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAttt                          <  2:2177965/91‑1 (MQ=255)
                                                                   aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                         >  2:2092047/1‑91 (MQ=255)
                                                                   aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                         <  1:2560467/91‑1 (MQ=255)
                                                                   aaCCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                         >  1:2511447/1‑91 (MQ=255)
                                                                   aaCCAAGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGAtttt                         >  1:770046/1‑91 (MQ=21)
                                                                     ccAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTtat                       >  1:2063291/1‑91 (MQ=255)
                                                                       aaGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg                     <  2:1266983/91‑1 (MQ=255)
                                                                       aaGATTCTCACAATTAGTAATGTGAAAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAg                     >  2:2392057/1‑91 (MQ=21)
                                                                        aGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGt                    >  1:824081/1‑91 (MQ=255)
                                                                         gATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGta                   >  2:949956/1‑91 (MQ=255)
                                                                           ttttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaat                 <  2:3788036/91‑1 (MQ=255)
                                                                            tttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaata                >  2:3687224/1‑91 (MQ=255)
                                                                             ttCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGtaataa               >  1:1686565/1‑91 (MQ=255)
                                                                              tCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAAc              >  1:9006/1‑91 (MQ=255)
                                                                                tcaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGCTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            <  2:3842939/90‑1 (MQ=255)
                                                                                acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            >  1:1592152/1‑91 (MQ=255)
                                                                                acaTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtt            >  2:3558765/1‑91 (MQ=255)
                                                                                 caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACttt           <  1:981237/91‑1 (MQ=255)
                                                                                 caTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACtta           <  2:2201414/91‑2 (MQ=255)
                                                                                   aTTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACttaa          <  1:918527/91‑3 (MQ=255)
                                                                                    ttAGTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTaa         <  1:3008246/91‑1 (MQ=255)
                                                                                    ttACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTaa         >  1:3024239/1‑91 (MQ=255)
                                                                                        tAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTa     >  2:3316871/1‑91 (MQ=255)
                                                                                         aaTTGTGAGAATCTTGGTTTTTTTTTTTTTTTTTGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAt    >  2:3800641/1‑91 (MQ=255)
                                                                                         aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTATGTTAt    >  2:664224/1‑91 (MQ=255)
                                                                                          aTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGTTAtt   <  1:1041243/91‑1 (MQ=255)
                                                                                           ttGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGTTAtta  <  2:3104357/91‑2 (MQ=255)
                                                                                           ttGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGTTAtta  <  2:2253618/91‑2 (MQ=255)
                                                                                           |                                                                                          
AAATTAACAGATGAAATTGAGCAAGCAAATAACTAAGATGCGAGACAATTTATAGTAGAACAATATAAACCAAGATTTTCACATTACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTG  >  USA300TCH1516_ALE/1421529‑1421709

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: