Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 740,322 G→C 12.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,3220GC12.4% 768.6 / 84.8 249intergenic (+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 (107/111);  new base C (11/20);  total (118/131)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

Only 100 of 270 total aligned reads displayed.

TATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAG  >  USA300TCH1516_ALE/740233‑740412
                                                                                         |                                                                                           
tataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatga                                                                                            <  1:48391/91‑1 (MQ=255)
tataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatga                                                                                            >  2:2021244/1‑91 (MQ=255)
 ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                           <  1:153118/91‑1 (MQ=255)
 ataaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                           <  2:350996/91‑1 (MQ=255)
 ataaataaGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgat                                                                                           <  1:2753343/91‑1 (MQ=255)
  taaatatGTATGATAGGGAAGGGTGAGATAATTTTTTTTATTATAATTAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACatgatg                                                                                          >  2:3870248/1‑91 (MQ=255)
  taaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAAAAAAAAAAAAAAATAAAAAACTAGCAATCTCACatgatg                                                                                          <  2:921067/91‑1 (MQ=255)
   aaatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGAtgt                                                                                         <  1:1408163/91‑1 (MQ=255)
    aatatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGa                                                                                       <  2:812587/91‑1 (MQ=255)
     atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                       >  1:1760047/1‑91 (MQ=255)
     atatGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGa                                                                                       <  1:2050990/91‑1 (MQ=255)
      tatgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGaa                                                                                     <  2:418244/91‑1 (MQ=255)
       atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                    <  1:1297042/91‑1 (MQ=255)
       atgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAAt                                                                                    <  1:2687409/91‑1 (MQ=255)
        tgtatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGagatt                                                                                   <  1:2564452/91‑5 (MQ=38)
         gtatgATAGGGAAGGGTGAGATAAATTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTg                                                                                  <  2:2481152/91‑1 (MQ=255)
          tatgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGc                                                                                 >  2:2472698/1‑91 (MQ=255)
           atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  1:1310746/91‑1 (MQ=255)
           atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  2:219560/91‑1 (MQ=255)
           atgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCt                                                                                <  2:502416/91‑1 (MQ=255)
            tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAATAATATAAAACTAGCAATCTCACACGATGTGAATTGCTa                                                                               >  2:1086614/1‑91 (MQ=255)
            tgATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTa                                                                               <  1:3531372/91‑1 (MQ=37)
             gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAg                                                                              <  2:560575/91‑1 (MQ=255)
             gATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAg                                                                              <  1:3250587/91‑1 (MQ=37)
              atatGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGt                                                                             <  1:3057052/91‑1 (MQ=255)
              aTAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGt                                                                             <  1:654878/91‑1 (MQ=37)
                aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                           >  1:402713/1‑91 (MQ=255)
                aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                           <  1:2004877/91‑3 (MQ=25)
                aGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGttt                                                                           <  1:113524/91‑3 (MQ=25)
                 gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtata                                                                          >  1:2047377/1‑91 (MQ=255)
                 gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                          <  1:3863674/91‑4 (MQ=25)
                 gggAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATTCACATCATGTGAGATTGCTAGtttt                                                                          <  1:577671/91‑4 (MQ=25)
                   gAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtatata                                                                        >  2:1713711/1‑91 (MQ=255)
                    aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATTTAAAACTAGCAATCTCACATGATGTGAATTGCTAGTatatat                                                                       >  2:2563616/1‑91 (MQ=255)
                    aaGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTatatat                                                                       <  2:1767990/91‑1 (MQ=255)
                     aGGGTGAGATAATTTTTTGTATCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                                      >  2:1877862/1‑91 (MQ=255)
                     aGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATc                                                                      >  1:3259029/1‑91 (MQ=255)
                        gTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGt                                                                   <  2:1988455/91‑1 (MQ=255)
                          gagaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTAc                                                                 >  1:1832055/1‑91 (MQ=255)
                           agaTAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACa                                                                >  2:1786989/1‑91 (MQ=255)
                               aaTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattt                                                            >  1:2540473/1‑91 (MQ=255)
                                aTTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAattta                                                           >  1:3511419/1‑91 (MQ=255)
                                 ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                          >  1:2829220/1‑91 (MQ=255)
                                 ttttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttat                                                          <  2:3770660/91‑1 (MQ=255)
                                  tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                         >  2:1177127/1‑91 (MQ=255)
                                  tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                         <  1:2727889/91‑1 (MQ=255)
                                  tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                         >  2:1241379/1‑91 (MQ=255)
                                  tttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttatt                                                         <  1:258684/91‑1 (MQ=255)
                                   ttttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAatttattt                                                        >  1:3567931/1‑91 (MQ=255)
                                    tttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                       >  2:2300579/1‑91 (MQ=255)
                                    ttcgTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAAtttattta                                                       >  2:2961047/4‑91 (MQ=255)
                                     ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                      >  2:3190031/1‑91 (MQ=255)
                                     ttGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTACttaa                                                      >  2:1643488/1‑91 (MQ=255)
                                     ttGTACCAAAATGAAAATATTTTTAAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaa                                                      <  2:1737631/91‑1 (MQ=255)
                                      tGTACCAAAATGAAAATTTTTTTCAAAAATTAACAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTACttaat                                                     >  2:3145839/1‑91 (MQ=255)
                                      tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                     >  1:374143/1‑91 (MQ=255)
                                      tGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaat                                                     >  2:2832030/1‑91 (MQ=255)
                                        tACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatta                                                   >  2:2341526/1‑91 (MQ=255)
                                        tACCAAAATGAAAATTTAAAAAAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaatta                                                   <  2:766301/91‑1 (MQ=255)
                                         aCCAAAATGAAAATTTTTTTCAAAAATAAAAAGCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                                  <  1:1043044/91‑1 (MQ=255)
                                         aCCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATttaattaa                                                  >  2:1370848/1‑91 (MQ=255)
                                          ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                 >  1:2677545/1‑91 (MQ=255)
                                          ccAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTtaattaat                                                 >  2:3867774/1‑91 (MQ=255)
                                            aaaaTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATgg                                               <  1:1284652/91‑1 (MQ=255)
                                               aTGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATg                                            >  1:2204592/1‑91 (MQ=255)
                                                 gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                          <  2:367244/91‑1 (MQ=255)
                                                 gAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGaa                                          <  1:76967/91‑1 (MQ=255)
                                                  aaaaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAAt                                         >  1:2480514/1‑91 (MQ=255)
                                                    aaTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGc                                       <  1:1878627/91‑1 (MQ=255)
                                                     aTTTTTTTCAAAAATAAAAAACTAGCAATCTCACATCATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCa                                      <  2:3055291/91‑1 (MQ=255)
                                                      tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                           >  1:3618057/1‑53 (MQ=255)
                                                      tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGtat                                                                           <  2:3618057/53‑1 (MQ=255)
                                                      tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCAt                                     <  1:1422688/91‑1 (MQ=255)
                                                      tttttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCAt                                     <  2:3217178/91‑1 (MQ=255)
                                                         ttttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGc                                  <  1:2961085/91‑1 (MQ=255)
                                                           ttCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTa                                >  1:2935161/1‑91 (MQ=255)
                                                            tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                               <  1:1103359/91‑1 (MQ=255)
                                                            tCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAg                               <  2:3347991/91‑1 (MQ=255)
                                                                  ttttaaaaCTAGCAATCTCACATGATGTTAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGATTGCATAGCTAGAAACtt                         >  2:1253179/5‑91 (MQ=255)
                                                                   tAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTc                        <  2:728359/91‑1 (MQ=255)
                                                                    aaaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCt                       <  2:1473577/91‑1 (MQ=255)
                                                                     aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                      <  1:806979/91‑1 (MQ=255)
                                                                     aaaaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTg                      >  1:1642269/1‑91 (MQ=255)
                                                                       aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                    <  2:3828094/91‑1 (MQ=255)
                                                                       aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                    >  1:981240/1‑91 (MQ=255)
                                                                       aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                    >  2:2789589/1‑91 (MQ=255)
                                                                       aaaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGaa                    >  1:1177883/1‑91 (MQ=255)
                                                                        aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                   <  2:3259928/91‑1 (MQ=255)
                                                                        aaCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAg                   >  2:215381/1‑91 (MQ=255)
                                                                         aCTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGc                  <  2:383311/91‑1 (MQ=255)
                                                                          cTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCt                 <  1:92859/91‑1 (MQ=255)
                                                                           tAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTg                <  1:2267215/91‑1 (MQ=255)
                                                                            aGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTgg               <  2:1681655/91‑1 (MQ=255)
                                                                              cAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGaa             >  2:2264495/1‑91 (MQ=255)
                                                                               aaTCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAAt            <  2:654415/91‑1 (MQ=255)
                                                                                    cacaTGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTAc       <  1:1256149/91‑1 (MQ=255)
                                                                                    cacaTCATGTGAGATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTa        >  2:2922697/1‑91 (MQ=255)
                                                                                      caTGATGTGAATTGCTGGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACgg     >  1:2632000/1‑91 (MQ=255)
                                                                                        tgatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTa   <  2:3457055/91‑1 (MQ=255)
                                                                                         gatgTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAg  <  2:106136/91‑1 (MQ=255)
                                                                                         |                                                                                           
TATAAATATGTATGATAGGGAAGGGTGAGATAATTTTTTGTACCAAAATGAAAATTTTTTTCAAAAATAAAAAACTAGCAATCTCACATGATGTGAATTGCTAGTATATATCAGTACAATTTATTTAATTAATGGATGAATGCATAGCTAGAAACTTCTGAAGCTGGAATTGTACGGTAG  >  USA300TCH1516_ALE/740233‑740412

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