Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,322 G→C 15.8% 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,3220GC15.8% 421.5 / 66.4 147intergenic (+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 (53/70);  new base C (5/18);  total (59/88)
Fisher's exact test for biased strand distribution p-value = 6.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 157 total aligned reads displayed.

ATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAC  >  USA300TCH1516_ALE/740234‑740407
                                                                                        |                                                                                      
ataactatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      >  1:3280148/1‑91 (MQ=255)
ataactatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      >  1:3280319/1‑91 (MQ=255)
ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      <  1:3750119/91‑1 (MQ=255)
ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      <  1:4124147/91‑1 (MQ=255)
ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      >  1:2166284/1‑91 (MQ=255)
ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      >  1:3661836/1‑91 (MQ=255)
ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                      >  1:1326980/1‑91 (MQ=255)
 taaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgatg                                                                                     <  1:16783/91‑1 (MQ=255)
  aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                    <  1:3413290/91‑1 (MQ=255)
  aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                    <  1:1186825/91‑1 (MQ=255)
  aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCAtgtg                                                                                   <  1:1355587/91‑1 (MQ=255)
  aaatatGTATGACAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                    >  1:3185119/1‑91 (MQ=255)
   aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgtg                                                                                   <  1:350793/91‑1 (MQ=255)
    atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                  >  1:1911679/1‑91 (MQ=255)
    atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGag                                                                                 <  1:871064/91‑2 (MQ=255)
       tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                              >  1:1474819/1‑91 (MQ=255)
       tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAtt                                                                              <  1:1481238/91‑1 (MQ=255)
       tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                              <  1:508215/91‑5 (MQ=38)
         tatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGc                                                                            <  1:583048/91‑1 (MQ=38)
          atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                           <  1:124209/91‑1 (MQ=255)
          atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                           <  1:3122852/91‑1 (MQ=255)
          atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCt                                                                           <  1:3458930/91‑1 (MQ=38)
           tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTa                                                                          <  1:3180990/91‑1 (MQ=37)
            gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                         <  1:2426744/91‑1 (MQ=37)
             aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGt                                                                        <  1:1823070/91‑1 (MQ=255)
             aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGt                                                                        <  1:1823111/91‑1 (MQ=255)
             aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                        <  1:2782893/91‑1 (MQ=37)
             aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTCTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                        <  1:3293667/91‑1 (MQ=37)
               aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                      <  1:904744/91‑3 (MQ=25)
               aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                      <  1:2709021/91‑3 (MQ=25)
                   aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTatatat                                                                  <  1:1445919/91‑1 (MQ=255)
                   aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTatatat                                                                  <  1:254240/91‑1 (MQ=255)
                   aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTatatat                                                                  <  1:2313241/91‑1 (MQ=255)
                     gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                                >  1:1740776/1‑91 (MQ=255)
                     gggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCa                                                                >  1:2237408/1‑91 (MQ=255)
                      ggTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAg                                                               <  1:3103577/91‑1 (MQ=255)
                        tGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTa                                                             <  1:3268982/91‑1 (MQ=255)
                          agaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACa                                                           >  1:867257/1‑91 (MQ=255)
                           gaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACaa                                                          >  1:3578999/1‑91 (MQ=255)
                           gaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACaa                                                          >  1:2924210/1‑91 (MQ=255)
                            aTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAat                                                         <  1:1451865/91‑1 (MQ=255)
                              aaTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattt                                                       <  1:3019890/91‑1 (MQ=255)
                               aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                      >  1:1052624/1‑91 (MQ=255)
                               aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                      >  1:1053167/1‑91 (MQ=255)
                                ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                     <  1:1343772/91‑1 (MQ=255)
                                 tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                    >  1:3095531/1‑91 (MQ=255)
                                  ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttattt                                                   >  1:2431347/1‑91 (MQ=255)
                                   tttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                  <  1:1729643/91‑1 (MQ=255)
                                   tttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                  <  1:760849/91‑1 (MQ=255)
                                    ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                 <  1:2427523/91‑1 (MQ=255)
                                    ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                 <  1:2047406/91‑1 (MQ=255)
                                    ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                 >  1:409023/1‑91 (MQ=255)
                                     tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                >  1:4063857/1‑91 (MQ=255)
                                      gTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatt                                               >  1:2926992/1‑91 (MQ=255)
                                      gTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatt                                               <  1:1584508/91‑1 (MQ=255)
                                       tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatta                                              >  1:1014524/1‑91 (MQ=255)
                                        aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                             <  1:2460895/91‑1 (MQ=255)
                                        aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                             <  1:2460850/91‑1 (MQ=255)
                                        aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                             <  1:1505665/91‑1 (MQ=255)
                                        aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                             <  1:2096406/91‑1 (MQ=255)
                                           aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                          >  1:1318388/1‑91 (MQ=255)
                                           aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                          <  1:1057515/91‑1 (MQ=255)
                                            aaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGa                                         <  1:1957492/91‑1 (MQ=255)
                                              aTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATg                                       <  1:2766517/91‑1 (MQ=255)
                                              aTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATg                                       >  1:2442591/1‑91 (MQ=255)
                                               tGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                      >  1:28316/1‑91 (MQ=255)
                                               tGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                      <  1:674571/91‑1 (MQ=255)
                                               tGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGa                                      <  1:674384/91‑1 (MQ=255)
                                                  aaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATg                                   >  1:2204321/1‑91 (MQ=255)
                                                   aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                  >  1:113464/1‑91 (MQ=255)
                                                   aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                  >  1:3134636/1‑91 (MQ=255)
                                                   aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                  >  1:2025674/1‑91 (MQ=255)
                                                   aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                  <  1:1937868/91‑1 (MQ=255)
                                                    aTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCa                                 <  1:695592/91‑1 (MQ=255)
                                                     tttttttCAAAAATAAAAAACTAGCAATCTCACATGACGTGAATTGCTAGTCTATATCAGTACAACTTATTTAATTAATGGATGAGTGCAt                                >  1:3855402/1‑91 (MQ=255)
                                                        ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                             <  1:588746/91‑1 (MQ=255)
                                                        ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                             <  1:1165537/91‑1 (MQ=255)
                                                           tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                          <  1:2595443/91‑1 (MQ=255)
                                                            cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                         <  1:3696760/91‑1 (MQ=255)
                                                            cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                         <  1:4114393/91‑1 (MQ=255)
                                                            cAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATAGATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGa                         >  1:4154030/1‑91 (MQ=255)
                                                                aaTAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACt                     <  1:3875604/91‑1 (MQ=255)
                                                                  tAAAAAACTGGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTc                   <  1:3971699/91‑1 (MQ=255)
                                                                  tAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTc                   <  1:2511083/91‑1 (MQ=255)
                                                                    aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                 >  1:1044386/1‑91 (MQ=255)
                                                                      aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa               >  1:3960139/1‑91 (MQ=255)
                                                                      aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa               >  1:2062230/1‑91 (MQ=255)
                                                                      aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa               >  1:1082307/1‑91 (MQ=255)
                                                                       aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg              >  1:1419995/1‑91 (MQ=255)
                                                                       aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg              >  1:2423838/1‑91 (MQ=255)
                                                                         cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt            <  1:2786332/91‑1 (MQ=255)
                                                                         cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt            <  1:3255568/91‑1 (MQ=255)
                                                                           aGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg          >  1:2009592/1‑91 (MQ=255)
                                                                            gCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGa         >  1:3642872/1‑91 (MQ=255)
                                                                             cAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGaa        >  1:4110615/1‑91 (MQ=255)
                                                                              caattcACATCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGaa        >  1:2867762/5‑91 (MQ=38)
                                                                              aaTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt       <  1:36183/91‑1 (MQ=255)
                                                                                   cacaTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc  >  1:1156080/1‑91 (MQ=255)
                                                                                    acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc  >  1:609504/1‑91 (MQ=255)
                                                                                    acaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc  >  1:3256854/1‑91 (MQ=255)
                                                                                        |                                                                                      
ATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAC  >  USA300TCH1516_ALE/740234‑740407

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