Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,765 T→A 5.3% intergenic (+32/‑119) btr → / → yxeP_1 HTH‑type transcriptional activator Btr/putative hydrolase YxeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE115,7650TA 5.3% 813.0 / 24.7 240intergenic (+32/‑119)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base T (101/125);  new base A (9/4);  total (111/129)
Fisher's exact test for biased strand distribution p-value = 9.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 243 total aligned reads displayed.

TTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCAT  >  USA300TCH1516_ALE/115675‑115855
                                                                                           |                                                                                          
 taaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATa                                                                                           <  1:3221844/91‑1 (MQ=255)
  aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                          >  2:1442506/1‑91 (MQ=255)
  aaaaTGATTACGAATGGGATTCACTTAACTGAAGTAAAACGTTACCCAAGTTCATAAACTGATCACAACTCACAAATCTTGATATACATaa                                                                                          >  2:1567242/1‑91 (MQ=255)
   aaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAt                                                                                         <  2:1112612/91‑1 (MQ=255)
   aaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAt                                                                                         >  1:2646974/1‑91 (MQ=255)
     aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACGCAAGTTCATCAAATGATCACAAATCACAATTTTTGATATACATAAttt                                                                                       <  2:3796744/91‑1 (MQ=255)
     aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAttt                                                                                       >  1:2006886/1‑91 (MQ=255)
      tGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtg                                                                                      <  1:247418/91‑1 (MQ=255)
      ggATGACGAAATGGATTCTCTTAATTGAAGTAATAGGGAACCCAAATTCATTAAATGAGCACAAATCACAAATTATGTATATCATAATTtg                                                                                      <  2:310274/90‑1 (MQ=12)
       gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                                     <  2:499822/91‑1 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                    <  1:1812110/91‑1 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                    >  2:3267883/1‑91 (MQ=255)
        aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                                    <  1:1177130/91‑1 (MQ=255)
         ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                   <  2:3853140/91‑1 (MQ=255)
         ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                   >  1:2752064/1‑91 (MQ=255)
         ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                   >  1:1081252/1‑91 (MQ=255)
         ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                                   >  2:2909005/1‑91 (MQ=255)
             gAATGGGATTAACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                               <  2:389682/91‑1 (MQ=255)
                tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                            >  1:1575646/1‑91 (MQ=255)
                 gggATTCACTTAATTGAAGTAAAACGTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtata                                                                          <  2:2505757/91‑1 (MQ=255)
                    aTTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATAtt                                                                        <  1:2493110/91‑1 (MQ=255)
                     ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                                       >  1:229053/1‑91 (MQ=255)
                     ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                                       >  1:2945141/1‑91 (MQ=255)
                      tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTTATATACATAATTTGTGATTTTTTATATTCa                                                                      >  2:2939551/1‑91 (MQ=255)
                       cACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCaa                                                                     <  2:3082750/91‑1 (MQ=255)
                         cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                                   <  1:2352412/91‑1 (MQ=255)
                         cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                                   <  1:1659883/91‑1 (MQ=255)
                           tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                                 <  2:918149/91‑1 (MQ=255)
                           tAATTGAAGAAAAACGCTACCCAAGGTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                                 >  2:270799/1‑91 (MQ=255)
                            aaTTGTAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                                <  2:3634993/91‑1 (MQ=255)
                              ttGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaa                                                              >  1:2443182/1‑91 (MQ=255)
                                gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaaaa                                                            <  2:1243496/91‑2 (MQ=255)
                                  aGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTg                                                          >  2:3403880/1‑91 (MQ=255)
                                  aGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTg                                                          >  2:1290546/1‑91 (MQ=255)
                                        aCGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaa                                                    >  1:3058862/1‑91 (MQ=255)
                                         cGTTAGCCAAGTTCATAAAATGATCACAAATCACGACTTTTGATATACGTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                                   <  2:224176/91‑1 (MQ=255)
                                         agTTACCCAAGTTCATAAAAAGATCAAAAATCACAAATTTTGATATACATAATTTGTAATATTTAATAATAAAAGAAAAAATAaaaaaaaa                                                   <  2:2049151/90‑1 (MQ=255)
                                          gTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                                  >  2:1058247/1‑91 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               <  1:3436090/91‑1 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               >  2:456809/1‑91 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               >  1:272983/1‑91 (MQ=255)
                                             aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                               >  2:3069762/1‑91 (MQ=255)
                                               ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                             <  2:2125955/91‑1 (MQ=255)
                                               ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                             >  2:2585675/1‑91 (MQ=255)
                                                cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTATAATTGCAAAAAATTAATgg                                            <  2:3461263/91‑1 (MQ=255)
                                                cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATgg                                            <  2:1149689/91‑1 (MQ=255)
                                                  aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                          >  1:2720526/1‑91 (MQ=255)
                                                   gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTa                                         <  2:2218184/91‑1 (MQ=255)
                                                    ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                                        >  2:1261954/1‑91 (MQ=255)
                                                     tCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAAc                                       >  2:588926/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    >  1:1061980/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    >  2:2424747/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    >  2:75741/1‑91 (MQ=255)
                                                        tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                                    <  2:3474056/91‑1 (MQ=255)
                                                         aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                   <  1:957417/91‑1 (MQ=255)
                                                         aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                   <  2:3701227/91‑1 (MQ=255)
                                                         aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                                   <  2:1125695/91‑1 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  >  1:1257161/1‑91 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  <  1:564515/91‑1 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  >  1:2851765/1‑91 (MQ=255)
                                                          aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                                  <  1:2865283/91‑1 (MQ=255)
                                                           aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                                 <  1:1875445/91‑1 (MQ=255)
                                                           aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                                 <  1:936213/91‑1 (MQ=255)
                                                           aaTGATCACAAATCACAAAGTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                                 <  1:3903004/91‑1 (MQ=255)
                                                            aTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtg                                <  1:752075/91‑1 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  1:2020922/1‑91 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            <  2:2278904/91‑1 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  2:3135198/1‑91 (MQ=255)
                                                                tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  1:1626305/1‑91 (MQ=255)
                                                                taaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                            >  2:2368792/5‑91 (MQ=21)
                                                                 cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                           >  2:3573358/1‑91 (MQ=255)
                                                                 cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                           >  2:2332064/1‑91 (MQ=255)
                                                                  acaaatcacaaatTTTGATATACATAATCTGCGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATAAATGGTTAACATCTCTGttgctg                          >  2:3019483/1‑91 (MQ=255)
                                                                  aaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGttgttg                          >  1:3346561/3‑91 (MQ=21)
                                                                   caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg                         <  1:1761096/91‑1 (MQ=255)
                                                                    aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCTAAGCTAAAGTTGAAAAAAATTATTGGTTAACATCTCTGTTGTTGGc                        <  1:3155091/91‑1 (MQ=255)
                                                                    aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGc                        >  2:1593482/1‑91 (MQ=255)
                                                                     aatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCa                       >  2:1814365/1‑91 (MQ=25)
                                                                      atcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCaa                      <  1:367097/91‑1 (MQ=255)
                                                                       tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat                     <  1:1532003/91‑1 (MQ=255)
                                                                       tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat                     <  2:2033917/91‑1 (MQ=255)
                                                                         acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat                   >  1:1027658/1‑91 (MQ=255)
                                                                         acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAGTAATGGTTAACATCTCTGTTGTTGGCAatat                   <  1:3533177/91‑1 (MQ=255)
                                                                            aatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaa                <  1:3043176/91‑1 (MQ=255)
                                                                             atTTTGATATACATAATTTGTGATTTTTTTTTTTTTTTGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat               >  2:3735445/1‑91 (MQ=255)
                                                                             atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat               >  1:2715978/1‑91 (MQ=255)
                                                                                ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaatgat            <  1:1619652/91‑4 (MQ=255)
                                                                                ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAaataat            >  1:315553/1‑91 (MQ=255)
                                                                                 tGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATg           <  2:3945670/91‑1 (MQ=255)
                                                                                  gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa          <  2:1776994/91‑1 (MQ=255)
                                                                                  gATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGa          >  1:3847680/1‑91 (MQ=255)
                                                                                   atatACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAt         <  2:2467319/91‑1 (MQ=255)
                                                                                    tataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGAtt        <  2:535774/91‑1 (MQ=255)
                                                                                     ataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTa       <  1:2547892/91‑1 (MQ=255)
                                                                                     ataCATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTa       <  2:724163/91‑1 (MQ=255)
                                                                                         cATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATc    <  2:2188540/91‑1 (MQ=255)
                                                                                          aTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTatcat   <  2:2135021/91‑5 (MQ=255)
                                                                                          aTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCa   <  1:3518720/91‑1 (MQ=255)
                                                                                          aTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCa   <  1:408189/91‑1 (MQ=255)
                                                                                          aTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCa   <  1:649553/91‑1 (MQ=255)
                                                                                           |                                                                                          
TTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAATGATTAATCAT  >  USA300TCH1516_ALE/115675‑115855

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: