Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,315 Δ1 bp 16.0% intergenic (+542/+46) pitA_1 → / ← sle1_2 Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE740,3150C.16.0% 905.4 / 163.3 287intergenic (+542/+46)pitA_1/sle1_2Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand):  ref base C (96/144);  new base . (5/41);  total (102/185)
Fisher's exact test for biased strand distribution p-value = 8.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 299 total aligned reads displayed.

ATACAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGT  >  USA300TCH1516_ALE/740225‑740405
                                                                                          |                                                                                           
aTACAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAtc                                                                                             <  2:3539406/91‑1 (MQ=255)
aTACAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAtc                                                                                             <  1:207642/91‑1 (MQ=255)
  aCAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAtctc                                                                                           <  1:2812568/91‑1 (MQ=255)
  aCAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAAtctc                                                                                           >  1:1760465/1‑91 (MQ=255)
   cAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTca                                                                                          >  2:3414650/1‑91 (MQ=255)
    aaaaTATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTcac                                                                                         >  2:1008898/1‑91 (MQ=255)
      aatataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACatc                                                                                      <  2:734514/91‑2 (MQ=255)
        aataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatga                                                                                     <  2:2696752/90‑1 (MQ=255)
         ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                    <  2:866273/91‑1 (MQ=255)
         ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                    <  2:2549669/91‑1 (MQ=255)
         ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                    <  2:2949419/91‑1 (MQ=255)
         ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                    <  1:1852429/91‑1 (MQ=255)
           aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                  <  2:1487012/91‑1 (MQ=255)
           aaatatGTATGATAGGGAAGGGAGAGATAATTTTTTGTACCAAAATGAAAATTAATTAAAAAAATAAAAAACTAGCAATTCACATCAtgtg                                                                                 <  2:3196504/91‑1 (MQ=255)
            aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgtg                                                                                 <  1:3473398/91‑1 (MQ=255)
            aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                <  1:1878913/91‑1 (MQ=255)
            aatatATATGATAGGGAAGGGTGAGATAATTTTTTGAACCAAAATGAAAATTTATTACAAAAAAAAAAAACTAGCAATCTCACATGAtgtg                                                                                 <  2:2446952/91‑1 (MQ=255)
            aagatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                <  2:2790716/88‑1 (MQ=255)
             atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                >  2:3899492/1‑91 (MQ=255)
             atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGag                                                                               <  1:3644150/91‑2 (MQ=255)
             atatGTATGATAGGGAAGGGTGAGAGAAGTTTAAGAACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGag                                                                               <  2:35514/91‑2 (MQ=255)
              tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                              <  2:1271369/91‑1 (MQ=255)
               atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                             <  2:3575578/91‑1 (MQ=255)
               atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                             <  2:1558468/91‑1 (MQ=255)
               atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagat                                                                             <  1:2159062/91‑4 (MQ=39)
               atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagat                                                                             <  1:3350024/91‑4 (MQ=39)
                tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                            <  1:1928475/91‑1 (MQ=255)
                tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                            <  1:1606894/91‑5 (MQ=38)
                tgtaggATAGGGAAGGGTTAGATAATTTTCCGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                            <  2:963204/91‑5 (MQ=37)
                 gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                           <  2:856651/91‑1 (MQ=255)
                 gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                           >  1:3841211/1‑91 (MQ=255)
                 gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                           >  1:588194/1‑91 (MQ=255)
                   atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                         <  2:1663954/91‑1 (MQ=255)
                   atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                         <  1:2465830/91‑1 (MQ=255)
                   atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                         <  1:2132497/91‑1 (MQ=255)
                   atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                         <  2:3295159/91‑1 (MQ=38)
                   atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAAGTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                         <  2:2350454/91‑1 (MQ=37)
                   atgATAGGGAAGGGTGAGATAATATTTTGTACCAAAATAAAAATAATAAAAAAAAAAAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                         <  2:3624519/91‑1 (MQ=255)
                    tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                        >  1:2552963/1‑91 (MQ=255)
                    tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                        >  1:1069154/1‑91 (MQ=255)
                    tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                        <  2:906650/91‑1 (MQ=255)
                     gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                       <  1:1469223/91‑1 (MQ=255)
                     gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                       <  1:3809577/91‑1 (MQ=255)
                     gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                       <  2:651079/91‑1 (MQ=37)
                      aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                      <  1:56159/91‑1 (MQ=37)
                        aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                    <  1:1169987/91‑3 (MQ=25)
                        aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                    <  2:2385811/91‑3 (MQ=25)
                         gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTGTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                   <  2:3463600/91‑4 (MQ=25)
                          ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatat                                                                  >  2:285652/1‑91 (MQ=255)
                           gAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatata                                                                 <  1:3578149/91‑1 (MQ=255)
                             aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                               >  1:1060873/1‑91 (MQ=255)
                             aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                               <  2:1811101/91‑1 (MQ=255)
                             aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTATAAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                               <  2:1826160/91‑1 (MQ=255)
                              gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                              <  1:832983/91‑1 (MQ=255)
                              gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                              >  2:1194957/1‑91 (MQ=255)
                               ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                             <  2:700640/91‑1 (MQ=255)
                                gagagaTAATGTTTTGTCATAAAATGAAAATTTTATTAAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGt                                                            <  2:1933291/89‑1 (MQ=255)
                                gTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGt                                                            <  2:3494451/91‑1 (MQ=255)
                                  gagaTAATTTTTTGTAGCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGTAATCTCACATGATGTGAATTGCTAGTATATATCAGTAc                                                          <  2:3780049/91‑1 (MQ=255)
                                    gaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACaa                                                        >  1:202166/1‑91 (MQ=255)
                                      tAATTATTTGTACCAAAATGAAAAATTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                      <  2:549829/91‑1 (MQ=255)
                                        aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                    <  1:118930/91‑1 (MQ=255)
                                         ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAAATAACAATCTCACATGATGTGAATTGCTAGTATATATAAGTACAatttat                                                   >  2:3374561/1‑91 (MQ=255)
                                           ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttattt                                                 <  2:1057549/91‑1 (MQ=255)
                                           ttatGTACCAAAATGAAAATTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                <  2:3046873/88‑1 (MQ=255)
                                            tttGTTCCAAAATGAAAATTTTTATAACAAATAAAAAACTAGCAATCTCACATGCTGTGAATTGCTAGTATATATCAGTACAAtttattta                                                >  2:2133132/1‑91 (MQ=255)
                                              tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                              <  1:3726219/91‑1 (MQ=255)
                                                tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatta                                            <  2:759524/91‑1 (MQ=255)
                                                  ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                          >  2:461308/1‑91 (MQ=255)
                                                  ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                          >  1:1710850/1‑91 (MQ=255)
                                                  ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                          >  2:3148865/1‑91 (MQ=255)
                                                     aaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                       <  1:1624736/91‑1 (MQ=255)
                                                     aaaTGAAAAATTTATTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                       <  2:1365473/91‑1 (MQ=255)
                                                      aaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGAt                                      <  1:3135558/91‑1 (MQ=255)
                                                        tGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                    >  1:2871970/1‑91 (MQ=255)
                                                         gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                   <  1:448293/91‑1 (MQ=255)
                                                         gAAAATTTTTTTCAAAAATAAAAAACTAGCAACCTCACATGATGTGAATTGCTAGTATATATCAGTACCATTTATTTAATTAATGGATGaa                                   >  2:3189759/1‑91 (MQ=255)
                                                           aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                 <  1:934749/91‑1 (MQ=255)
                                                            aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                <  2:1642004/91‑1 (MQ=255)
                                                               ttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATa                             >  1:660557/1‑91 (MQ=255)
                                                                 ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                           <  1:508201/91‑1 (MQ=255)
                                                                     aaaaaaTAAAAGACTAGCAATCTTACATGATGTGAATTGCTAGTATATATCAGTACAATGTATTTAATTAATGGATGAATGCATAGCTAGa                       <  2:3017823/90‑1 (MQ=255)
                                                                      aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                      >  1:2722792/1‑91 (MQ=255)
                                                                      aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                      >  2:1209295/1‑91 (MQ=255)
                                                                       aaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                     >  1:2211353/1‑91 (MQ=255)
                                                                          aTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACtt                  <  2:3915812/91‑1 (MQ=255)
                                                                          aTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACtt                  <  2:1767387/91‑1 (MQ=255)
                                                                            aaaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCt                <  1:240639/91‑1 (MQ=255)
                                                                             aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg               >  1:1415150/1‑91 (MQ=255)
                                                                             aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATCGATGAATGCATAGCTAGAAACTTCTg               <  2:2249643/91‑1 (MQ=255)
                                                                               aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa             >  2:27407/1‑91 (MQ=255)
                                                                                aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg            <  1:755246/91‑1 (MQ=255)
                                                                                aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg            >  2:2031211/1‑91 (MQ=255)
                                                                                 aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc           >  1:3939159/1‑91 (MQ=255)
                                                                                 aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc           <  1:285652/91‑1 (MQ=255)
                                                                                 aCTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc           >  2:2176624/1‑91 (MQ=37)
                                                                                  cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt          <  1:3393742/91‑1 (MQ=255)
                                                                                   tAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTg         <  2:1981956/91‑1 (MQ=255)
                                                                                    aGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg        >  2:573699/1‑91 (MQ=38)
                                                                                        aTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAtt    <  1:1550188/91‑1 (MQ=255)
                                                                                          |                                                                                           
ATACAAAATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGT  >  USA300TCH1516_ALE/740225‑740405

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: