Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,322 G→C 21.4% intergenic (+549/+39) 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,3220GC21.4% 657.5 / 188.0 287intergenic (+549/+39)pitA_1/sle1_2Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand):  ref base G (90/135);  new base C (23/38);  total (114/173)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 318 total aligned reads displayed.

ATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAG  >  USA300TCH1516_ALE/740232‑740412
                                                                                          |                                                                                           
 tataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatga                                                                                            >  2:2669013/1‑91 (MQ=255)
 tataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatga                                                                                            <  2:2696779/91‑1 (MQ=255)
  ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                           <  2:866273/91‑1 (MQ=255)
  ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                           <  2:1358663/91‑1 (MQ=255)
  ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACatcatg                                                                                          <  1:550338/91‑1 (MQ=255)
   taaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgatg                                                                                          >  1:273017/1‑91 (MQ=255)
    aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                         <  2:1487012/91‑1 (MQ=255)
    aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCAtgtg                                                                                        <  2:3410215/91‑1 (MQ=255)
    aaatatGTATGATAGGGAAGGGAGAGATAATTTTTTGTACCAAAATGAAAATTAATTAAAAAAATAAAAAACTAGCAATTCACATCAtgtg                                                                                        <  2:3196504/91‑1 (MQ=255)
     aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                       <  1:1878913/91‑1 (MQ=255)
     aatagGTATGATAGGGAAGGGGGAGATAATTTTTTGTACCAAAGTGAAAATTTTTGTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                       <  2:2790676/91‑1 (MQ=255)
      atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                       >  2:3899502/1‑91 (MQ=255)
       tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                     >  1:1816767/1‑91 (MQ=255)
       tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                     >  2:3115832/1‑91 (MQ=255)
        atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                    <  2:1558468/91‑1 (MQ=255)
        atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagat                                                                                    <  1:3350024/91‑4 (MQ=39)
        atgtatgATAGGGAAGGGTGAGATAATTTTGTGTACCAAAATGAAAATTATTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                    <  2:3575473/91‑1 (MQ=255)
         tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCTTAATTAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                                   >  2:3528040/1‑91 (MQ=255)
         tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                                   <  1:1928475/91‑1 (MQ=255)
         tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                                   <  1:1606894/91‑5 (MQ=38)
         tgtaggATAGGGAAGGGTTAGATAATTTTCCGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                                   <  2:963204/91‑5 (MQ=37)
         tggatgATAGTGAAGGGTAAGATAAGTGTAAGTACCAAAATAAAAAGTTATTAAAAAAATAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                                   <  2:1264715/88‑1 (MQ=255)
          gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                  >  2:3813783/1‑91 (MQ=255)
          gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                  <  2:856651/91‑1 (MQ=255)
          gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                  >  1:3841594/1‑91 (MQ=255)
            atgATAGGGAATGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTGCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                                <  2:3295160/91‑1 (MQ=37)
            atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  1:586963/91‑1 (MQ=255)
            atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  2:1597886/91‑1 (MQ=255)
            atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                >  1:157061/1‑91 (MQ=255)
            atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAAGTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                                <  2:2350454/91‑1 (MQ=37)
            atgATAGAGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  2:1824145/91‑1 (MQ=255)
            ataaTAGGGAAGGGTGAGAAAATTTTTTGAAACAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  2:2368705/88‑1 (MQ=255)
             tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                               <  2:906650/91‑1 (MQ=255)
              gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                              <  1:1470865/91‑1 (MQ=255)
              gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                              <  2:651079/91‑1 (MQ=37)
                 aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                           <  1:880013/91‑1 (MQ=255)
                  gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTTAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtata                                                                          >  2:1905149/1‑91 (MQ=255)
                  gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                          <  2:2936912/91‑4 (MQ=25)
                  gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTGTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                          <  2:3463600/91‑4 (MQ=25)
                   ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCGCACATGATGTGAATTGCTAGtatat                                                                         >  1:3666721/1‑91 (MQ=255)
                   ggAAGAGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatat                                                                         <  2:3592254/91‑1 (MQ=255)
                      aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                                      >  1:1061270/1‑91 (MQ=255)
                        ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                                    <  2:700640/91‑1 (MQ=255)
                         gTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGt                                                                   >  1:2088234/1‑91 (MQ=255)
                          tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTa                                                                  <  1:1279931/91‑1 (MQ=255)
                             gaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACaa                                                               >  1:202166/1‑91 (MQ=255)
                               tAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                             <  1:3626535/91‑1 (MQ=255)
                               tAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                             <  1:1185867/91‑1 (MQ=255)
                               tAATTATTTGTACCAAAATGAAAAATTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                             <  2:549829/91‑1 (MQ=255)
                                   tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                         <  1:2537017/91‑1 (MQ=255)
                                    ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttattt                                                        <  2:1057549/91‑1 (MQ=255)
                                    ttatGTACCAAAATGAAAATTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                       <  2:3046873/88‑1 (MQ=255)
                                     tttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                       <  1:213969/91‑1 (MQ=255)
                                      ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                      >  2:2133492/1‑91 (MQ=255)
                                           ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                 >  2:461308/1‑91 (MQ=255)
                                           ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                 >  1:1710850/1‑91 (MQ=255)
                                             aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                               >  1:3754306/1‑91 (MQ=255)
                                              aaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                              <  1:1624736/91‑1 (MQ=255)
                                               aaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGAt                                             <  1:3135558/91‑1 (MQ=255)
                                                 tGAAAATTTTTTTCAAAATTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                           >  2:2838238/1‑91 (MQ=255)
                                                  gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                          <  1:448293/91‑1 (MQ=255)
                                                    aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                        <  1:934749/91‑1 (MQ=255)
                                                     aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                       >  2:251353/1‑91 (MQ=255)
                                                       tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCAt                                     <  2:778898/91‑1 (MQ=255)
                                                         tttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAg                                   >  2:1736667/1‑91 (MQ=255)
                                                          ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                  <  1:508201/91‑1 (MQ=255)
                                                          ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                  <  2:3176614/91‑1 (MQ=255)
                                                           tttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCt                                 <  2:1648151/91‑1 (MQ=255)
                                                             tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                               <  2:3489920/91‑1 (MQ=255)
                                                              cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                              <  2:2394916/91‑1 (MQ=255)
                                                              cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                              >  1:3819265/1‑91 (MQ=255)
                                                              aaaaaaTAAAAGACTAGCAATCTTACATGATGTGAATTGCTAGTATATATCAGTACAATGTATTTAATTAATGGATGAATGCATAGCTAGa                              <  2:3017823/90‑1 (MQ=255)
                                                               aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                             >  2:1209295/1‑91 (MQ=255)
                                                                aaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaaa                            >  1:2211353/1‑91 (MQ=255)
                                                                     aaaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCt                       <  1:240639/91‑1 (MQ=255)
                                                                      aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATCGATGAATGCATAGCTAGAAACTTCTg                      <  2:2249643/91‑1 (MQ=255)
                                                                        aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                    <  1:2585325/91‑1 (MQ=255)
                                                                         aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                   <  1:3533813/91‑1 (MQ=255)
                                                                         aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                   <  1:2897509/91‑1 (MQ=255)
                                                                         aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                   <  1:755246/91‑1 (MQ=255)
                                                                          aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                  <  1:285652/91‑1 (MQ=255)
                                                                           cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                 <  2:2041835/91‑1 (MQ=255)
                                                                           cTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                 >  1:1429728/1‑91 (MQ=37)
                                                                            tAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTg                <  2:1983948/91‑1 (MQ=255)
                                                                             aGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg               >  1:1150749/1‑91 (MQ=255)
                                                                             aGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg               >  2:573699/1‑91 (MQ=38)
                                                                                 aTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAtt           >  2:1780123/1‑91 (MQ=255)
                                                                                   ctcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGt         <  1:3393298/91‑1 (MQ=255)
                                                                                    tcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        <  1:400846/91‑1 (MQ=255)
                                                                                    tcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        >  2:709371/1‑91 (MQ=255)
                                                                                     cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        >  2:3390518/1‑91 (MQ=255)
                                                                                     cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        >  1:65698/1‑91 (MQ=255)
                                                                                     cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        >  1:698020/1‑91 (MQ=255)
                                                                                      acaTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGACTTgaaaa      >  2:1188588/1‑87 (MQ=255)
                                                                                      acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc       >  2:1945857/1‑91 (MQ=255)
                                                                                      acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc       >  1:3428182/1‑91 (MQ=255)
                                                                                       caTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACgg     <  2:2205467/91‑1 (MQ=255)
                                                                                        atgatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGt    <  2:3723176/91‑1 (MQ=255)
                                                                                          gatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg  <  1:81691/91‑1 (MQ=255)
                                                                                          gatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg  <  1:1012009/91‑1 (MQ=255)
                                                                                          |                                                                                           
ATATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAG  >  USA300TCH1516_ALE/740232‑740412

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: