Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,328:1 +G 15.2% intergenic (+555/+33) 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,3281.G15.2% 973.1 / 209.5 282intergenic (+555/+33)pitA_1/sle1_2Low‑affinity inorganic phosphate transporter 1/N‑acetylmuramoyl‑L‑alanine amidase sle1
Reads supporting (aligned to +/- strand):  ref base . (94/145);  new base G (27/16);  total (121/161)
Fisher's exact test for biased strand distribution p-value = 6.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 303 total aligned reads displayed.

AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCAT  >  USA300TCH1516_ALE/740237‑740417
                                                                                            |                                                                                         
aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                            <  1:477331/91‑1 (MQ=255)
aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                            <  1:1864221/91‑1 (MQ=255)
aatagGTATGATAGGGAAGGGGGAGATAATTTTTTGTACCAAAGTGAAAATTTTTGTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                            <  2:2790676/91‑1 (MQ=255)
 atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                            >  2:3899502/1‑91 (MQ=255)
 atatGTATGATAGGGAAGGGTGAGATAAATATGTGTACCAAAATGAAAATATTTTGCAAAAATAAAAAACTAGCAATTCACATCATGTGag                                                                                           <  2:1835226/91‑2 (MQ=255)
  tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                          <  2:1271369/91‑1 (MQ=255)
   atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                         <  2:1558468/91‑1 (MQ=255)
   atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagat                                                                                         <  1:3350024/91‑4 (MQ=39)
    tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCTTAATTAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                                        >  2:3528040/1‑91 (MQ=255)
    tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                                        <  1:1606894/91‑5 (MQ=38)
    tgtaggATAGGGAAGGGTTAGATAATTTTCCGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                                        <  2:963204/91‑5 (MQ=37)
     gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                       >  1:3840998/1‑91 (MQ=255)
     gtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                       <  2:856609/91‑1 (MQ=255)
       atgATAGGGAATGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTGCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                                     <  2:3295160/91‑1 (MQ=37)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                     <  2:1597886/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                     <  1:586963/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                     <  1:2132497/91‑1 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                     >  1:157061/1‑91 (MQ=255)
       atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAAGTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                                     <  2:2350454/91‑1 (MQ=37)
        tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTa                                                                                    >  1:1069154/1‑91 (MQ=255)
         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                                   <  2:3032962/91‑1 (MQ=255)
         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                                   <  1:1470865/91‑1 (MQ=255)
         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                                   <  2:651079/91‑1 (MQ=37)
         gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                                   <  2:1465081/91‑1 (MQ=37)
          aTAGGGATGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                                  <  1:1972572/91‑1 (MQ=37)
          aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGt                                                                                  <  1:2046921/91‑1 (MQ=255)
           tAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtt                                                                                 <  1:2404601/91‑2 (MQ=37)
            aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                                <  1:1169987/91‑3 (MQ=25)
             gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTGTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                               <  2:3463600/91‑4 (MQ=25)
              ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatat                                                                              >  2:285652/1‑91 (MQ=255)
              ggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCGCACATGATGTGAATTGCTAGtatat                                                                              >  1:3666721/1‑91 (MQ=255)
                  gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                                          >  2:2077598/1‑91 (MQ=255)
                  gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                                          <  1:832983/91‑1 (MQ=255)
                   ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                                         <  2:700640/91‑1 (MQ=255)
                     tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTa                                                                       <  2:893130/91‑1 (MQ=255)
                     tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTa                                                                       <  1:1279931/91‑1 (MQ=255)
                      gagaTAATTTTTTGTAGCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGTAATCTCACATGATGTGAATTGCTAGTATATATCAGTAc                                                                      <  2:3780049/91‑1 (MQ=255)
                          tAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                                  <  1:3626535/91‑1 (MQ=255)
                          tAATTATTTGTACCAAAATGAAAAATTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatt                                                                  <  2:549829/91‑1 (MQ=255)
                               ttttGGACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATAGGTTCAatttattg                                                             <  2:1057538/91‑2 (MQ=255)
                                tttGTACCAAAATGAAAATTTTTTACAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                            <  2:918652/91‑1 (MQ=255)
                                  tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                          >  2:1125598/1‑91 (MQ=255)
                                    tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatta                                                        <  2:759524/91‑1 (MQ=255)
                                     aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                                       >  1:2169993/1‑91 (MQ=255)
                                      ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                      >  2:3148865/1‑91 (MQ=255)
                                      ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                      >  1:1710850/1‑91 (MQ=255)
                                      ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                      >  2:461308/1‑91 (MQ=255)
                                         aaaTGCAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                                   <  2:3202513/91‑1 (MQ=255)
                                         aaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                                   <  1:1624736/91‑1 (MQ=255)
                                          aaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGAt                                                  <  1:3135558/91‑1 (MQ=255)
                                            tGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                                <  2:213369/91‑1 (MQ=255)
                                               aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                             >  1:2833690/1‑91 (MQ=255)
                                                aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                            >  2:251353/1‑91 (MQ=255)
                                                    tttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAg                                        >  2:1736667/1‑91 (MQ=255)
                                                      tttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCt                                      <  2:1648151/91‑1 (MQ=255)
                                                       ttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTa                                     <  2:2977329/91‑1 (MQ=255)
                                                        tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                                    <  2:1815574/91‑1 (MQ=255)
                                                        tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                                    <  2:3489920/91‑1 (MQ=255)
                                                        tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                                    <  1:2539592/91‑1 (MQ=255)
                                                         cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                                   <  2:2394916/91‑1 (MQ=255)
                                                          aaaaaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGaa                                  <  1:3455299/91‑1 (MQ=255)
                                                               tAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTc                             <  2:3577810/91‑1 (MQ=255)
                                                                 agaaaCTAGCAATCTCACATGATGTTAATTTCTAGTGTAGATCAGTACAATTTGTTTAATGAATGGATGAATGCAGAGCTAGAAACTTcgg                           <  2:3565910/89‑3 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                           <  1:1923776/91‑1 (MQ=255)
                                                                 aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATCGATGAATGCATAGCTAGAAACTTCTg                           <  2:2249643/91‑1 (MQ=255)
                                                                    aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                        <  1:3533813/91‑1 (MQ=255)
                                                                    aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                        <  1:755246/91‑1 (MQ=255)
                                                                     aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                       >  1:3939159/1‑91 (MQ=255)
                                                                      cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                      <  1:3393742/91‑1 (MQ=255)
                                                                      cTAGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                      >  1:1429728/1‑91 (MQ=37)
                                                                       tAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTg                     <  2:1983948/91‑1 (MQ=255)
                                                                        aGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg                    <  2:120195/91‑1 (MQ=255)
                                                                        aGCAATTCACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg                    >  2:573699/1‑91 (MQ=38)
                                                                           aaTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt                 >  1:2865324/1‑91 (MQ=255)
                                                                           aaTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt                 >  1:424493/1‑91 (MQ=255)
                                                                            aTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAtt                >  2:1780123/1‑91 (MQ=255)
                                                                             attcACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAtt                >  2:3064543/3‑91 (MQ=39)
                                                                               tcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa             <  1:528162/91‑1 (MQ=255)
                                                                               tcACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa             >  2:709371/1‑91 (MQ=255)
                                                                               tcACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGATGGAATTGt              >  1:3783609/1‑91 (MQ=255)
                                                                                cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa             >  1:65698/1‑91 (MQ=255)
                                                                                cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa             >  2:3390518/1‑91 (MQ=255)
                                                                                cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa             >  1:698020/1‑91 (MQ=255)
                                                                                 acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc            >  2:3709661/1‑91 (MQ=255)
                                                                                 acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc            >  2:1945857/1‑91 (MQ=255)
                                                                                  caTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACgg          <  2:2205467/91‑1 (MQ=255)
                                                                                   atcatgTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACgg          >  1:2439305/4‑91 (MQ=255)
                                                                                    acatgTGAGATCGCTAGTATATACAAATACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGATGGAATTATACGGt         >  2:2736088/3‑91 (MQ=255)
                                                                                     gatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg       >  1:899935/1‑91 (MQ=255)
                                                                                     gatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg       >  1:1055412/1‑91 (MQ=255)
                                                                                      atgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGt      <  1:2997127/91‑1 (MQ=255)
                                                                                      atgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGt      <  1:1003307/91‑1 (MQ=255)
                                                                                      atgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGt      <  1:810619/91‑1 (MQ=255)
                                                                                      atgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGt      <  2:2030201/91‑1 (MQ=255)
                                                                                       tgtgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGtt     >  2:2379396/1‑91 (MQ=255)
                                                                                        gtgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTc    <  1:2588532/91‑1 (MQ=255)
                                                                                        gtgAATTGCTAGTATATATCAGGACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTc    <  2:3180945/91‑1 (MQ=255)
                                                                                         tgAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCa   <  1:981179/91‑1 (MQ=255)
                                                                                          gAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCat  >  2:126003/1‑91 (MQ=255)
                                                                                          gAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCat  >  1:1188156/1‑91 (MQ=255)
                                                                                            |                                                                                         
AATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAGTTCAT  >  USA300TCH1516_ALE/740237‑740417

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: