Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 115,755 T→A 6.1% intergenic (+22/‑129) 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,7550TA 6.1% 726.1 / 25.8 223intergenic (+22/‑129)btr/yxeP_1HTH‑type transcriptional activator Btr/putative hydrolase YxeP
Reads supporting (aligned to +/- strand):  ref base T (106/103);  new base A (9/5);  total (115/108)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 226 total aligned reads displayed.

ATTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAAT  >  USA300TCH1516_ALE/115667‑115845
                                                                                        |                                                                                           
attaaacaataaaATGATTACGAAGGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAATTGATCACAAATCACAAATTTTg                                                                                           <  1:1143383/91‑1 (MQ=255)
 ttaaacattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACATTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGa                                                                                          >  2:3472748/1‑91 (MQ=255)
     acattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGAtata                                                                                      <  1:533247/91‑1 (MQ=255)
      cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                    <  2:2631052/91‑1 (MQ=255)
      cattaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATAc                                                                                    <  2:1803630/91‑1 (MQ=255)
       attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                   >  1:972411/1‑91 (MQ=255)
       attaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACa                                                                                   >  2:2785057/1‑91 (MQ=255)
        ttaaaATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACAt                                                                                  <  2:1121115/91‑1 (MQ=255)
          aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                >  2:1442506/1‑91 (MQ=255)
          aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                >  1:247423/1‑91 (MQ=255)
          aaaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATaa                                                                                >  1:2184673/1‑91 (MQ=255)
           aaaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAt                                                                               >  1:2646974/1‑91 (MQ=255)
            aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                              <  1:3301088/91‑1 (MQ=255)
            aaTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAAtt                                                                              >  1:171912/1‑91 (MQ=255)
             aTGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACGCAAGTTCATCAAATGATCACAAATCACAATTTTTGATATACATAAttt                                                                             <  2:3796744/91‑1 (MQ=255)
              ggATGACGAAATGGATTCTCTTAATTGAAGTAATAGGGAACCCAAATTCATTAAATGAGCACAAATCACAAATTATGTATATCATAATTtg                                                                            <  2:310274/90‑1 (MQ=12)
               gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                           <  1:2456425/91‑1 (MQ=255)
               gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                           <  2:499822/91‑1 (MQ=255)
               gATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgt                                                                           <  2:2360374/91‑1 (MQ=255)
                aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                          <  1:1812110/91‑1 (MQ=255)
                aTTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTtgtg                                                                          >  1:2760846/1‑91 (MQ=255)
                 ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                         >  1:1081252/1‑91 (MQ=255)
                 ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                         >  2:3078795/1‑91 (MQ=255)
                 ttACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGa                                                                         <  2:3853140/91‑1 (MQ=255)
                   aCGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtt                                                                       >  1:2015789/1‑91 (MQ=255)
                    cGAATGGTATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                      <  2:3521354/91‑1 (MQ=255)
                    cGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAttt                                                                      <  1:3606987/91‑1 (MQ=255)
                     gAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                     <  2:749764/91‑1 (MQ=255)
                     gAATGGGATTCACTTAATTGAAGGAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTGGATATACATAATTTGTGAtttt                                                                     <  2:568538/91‑1 (MQ=255)
                     gAATGGGATTAACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGAtttt                                                                     <  2:389682/91‑1 (MQ=255)
                        tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                  >  1:1576831/1‑91 (MQ=255)
                        tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                  <  2:1073117/91‑1 (MQ=255)
                        tGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTta                                                                  >  1:3465540/1‑91 (MQ=255)
                         gggATTCACTTAATTGAAGTAAAACGTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTtata                                                                <  2:2505757/91‑1 (MQ=255)
                             ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                             >  1:229053/1‑91 (MQ=255)
                             ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                             >  1:2945141/1‑91 (MQ=255)
                             ttCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTc                                                             <  1:825794/91‑1 (MQ=255)
                              tCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCa                                                            >  2:1978901/1‑91 (MQ=255)
                                 cTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAg                                                         <  1:1659883/91‑1 (MQ=255)
                                   tAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                       <  2:918149/91‑1 (MQ=255)
                                   tAATTGAAGAAAAACGCTACCCAAGGTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCt                                                       >  2:270799/1‑91 (MQ=255)
                                    aaTTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTa                                                      >  1:385859/1‑91 (MQ=255)
                                      ttGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTaaa                                                    >  1:2723350/1‑91 (MQ=255)
                                        gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                  <  1:3396856/91‑1 (MQ=255)
                                        gAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAAt                                                  <  2:1244120/91‑1 (MQ=255)
                                         aaGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATCCATAATTTGTGATTTTTTATATTCAAAGCTAAAAtt                                                 >  2:3117801/1‑91 (MQ=255)
                                            tAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCa                                              <  2:1021279/91‑1 (MQ=255)
                                                 cGTTAGCCAAGTTCATAAAATGATCACAAATCACGACTTTTGATATACGTAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                         <  2:224176/91‑1 (MQ=255)
                                                 cGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCaaaaaa                                         <  2:2049505/91‑1 (MQ=255)
                                                  gTTACCCAAGTTAATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAt                                        >  1:3048491/1‑91 (MQ=255)
                                                   ttACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATATGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAAtt                                       <  2:1146926/91‑1 (MQ=255)
                                                     aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                     >  1:272983/1‑91 (MQ=255)
                                                     aCCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                     >  2:456809/1‑91 (MQ=255)
                                                     aCCCAAGTAAATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTaa                                     <  2:1234667/91‑1 (MQ=255)
                                                       ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                   <  2:1605305/91‑1 (MQ=255)
                                                       ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                   <  2:2125955/91‑1 (MQ=255)
                                                       ccAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATg                                   >  2:2585675/1‑91 (MQ=255)
                                                        cAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTATAATTGCAAAAAATTAATgg                                  <  2:3461263/91‑1 (MQ=255)
                                                         aaGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGt                                 >  2:1531708/1‑91 (MQ=255)
                                                         aaGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATAGCAAAAAATTAATGGt                                 <  2:3147036/91‑1 (MQ=255)
                                                          aGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGtt                                <  1:3211623/91‑1 (MQ=255)
                                                           gTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATAGAAAAAAATTAATGGTTa                               <  2:2218081/91‑1 (MQ=255)
                                                            ttCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTaa                              <  2:1576831/91‑1 (MQ=255)
                                                               aTAAAATGATCACACATCACCAACTTTGATATACATAATTTGTGATTTATTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAt                           >  2:2436684/1‑91 (MQ=255)
                                                                tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                          >  1:3124568/1‑91 (MQ=255)
                                                                tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                          >  1:1214951/1‑91 (MQ=255)
                                                                tAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtc                          >  2:3333719/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)
                                                                 aaaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtct                         <  2:3701227/91‑1 (MQ=255)
                                                                  aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                        >  1:1257161/1‑91 (MQ=255)
                                                                  aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                        >  1:2851765/1‑91 (MQ=255)
                                                                  aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                        <  1:2865283/91‑1 (MQ=255)
                                                                  aaaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACAtctc                        <  1:564515/91‑1 (MQ=255)
                                                                   aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                       <  1:936213/91‑1 (MQ=255)
                                                                   aaTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                       <  1:1875445/91‑1 (MQ=255)
                                                                   aaTGATCACAAATCACAAAGTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATctct                       <  1:3903004/91‑1 (MQ=255)
                                                                    aTGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtg                      <  1:752075/91‑1 (MQ=255)
                                                                        tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                  <  2:3061715/91‑1 (MQ=255)
                                                                        tcacaaatcacaaaTTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCtgttgt                  >  1:1626305/1‑91 (MQ=255)
                                                                         cacaaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                 >  2:2332064/1‑91 (MQ=255)
                                                                         aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                 >  2:2926623/4‑91 (MQ=21)
                                                                         aaaaaatcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTgttgtt                 >  2:1291983/4‑91 (MQ=21)
                                                                           caaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTgg               <  1:1761096/91‑1 (MQ=255)
                                                                            aaatcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCTAAGCTAAAGTTGAAAAAAATTATTGGTTAACATCTCTGTTGTTGGc              <  1:3155091/91‑1 (MQ=255)
                                                                               tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat           <  2:2034032/91‑1 (MQ=255)
                                                                               tcacaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat           <  1:1532003/91‑1 (MQ=255)
                                                                               tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAATAATTAATGGTTAACATCTCTGTTGTTGGCAat           >  2:1275382/1‑91 (MQ=25)
                                                                               tcacaaatTATGTATATCAAAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAat           >  1:2567663/1‑91 (MQ=25)
                                                                                 acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat         >  2:360846/1‑91 (MQ=255)
                                                                                 acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat         >  2:2264500/1‑91 (MQ=255)
                                                                                 acaaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAatat         >  1:1027658/1‑91 (MQ=255)
                                                                                  caaatTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAAtata        <  1:2369162/91‑1 (MQ=255)
                                                                                     atTTTGATATACATAATTTGTGATTTTTTTTTTTTTTTGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat     >  2:3735445/1‑91 (MQ=255)
                                                                                     atTTTGATATACATAATTTGTGATTTTTTTTATTTTAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat     >  2:2962862/1‑91 (MQ=255)
                                                                                     atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat     <  1:1381099/91‑1 (MQ=255)
                                                                                     atTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaat     >  1:2715992/1‑91 (MQ=255)
                                                                                        ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATataaatgat  <  1:1619652/91‑4 (MQ=255)
                                                                                        ttGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAaataat  >  1:315553/1‑91 (MQ=255)
                                                                                        |                                                                                           
ATTAAACATTAAAATGATTACGAATGGGATTCACTTAATTGAAGTAAAACGTTACCCAAGTTCATAAAATGATCACAAATCACAAATTTTGATATACATAATTTGTGATTTTTTATATTCAAAGCTAAAATTGCAAAAAATTAATGGTTAACATCTCTGTTGTTGGCAATATAAATAAT  >  USA300TCH1516_ALE/115667‑115845

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: