Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,435,507 A→C 9.3% L331W (TTG→TGG)  yifK ← putative transport protein YifK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,435,5070AC 9.3% 354.5 / 22.7 107L331W (TTG→TGG) yifKputative transport protein YifK
Reads supporting (aligned to +/- strand):  ref base A (42/55);  new base C (7/3);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 116 total aligned reads displayed.

ATACTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAA  >  USA300TCH1516_ALE/2435417‑2435597
                                                                                          |                                                                                          
aTACTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACa                                                                                            <  1:789937/91‑1 (MQ=255)
aTACTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACa                                                                                            <  1:121273/91‑1 (MQ=255)
 tACTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACaa                                                                                           >  1:2561800/1‑91 (MQ=255)
  aCTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAAt                                                                                          >  1:2190715/1‑91 (MQ=255)
   ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                         <  1:365608/91‑1 (MQ=255)
   ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                         <  1:1574718/91‑1 (MQ=255)
   ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                         >  1:333922/1‑91 (MQ=255)
   ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                         >  1:1791424/1‑91 (MQ=255)
     ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                       <  1:677051/91‑1 (MQ=255)
     ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                       >  1:3673857/1‑91 (MQ=255)
     ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                       >  1:1373809/1‑91 (MQ=255)
        cGCTCCATCGATATATAAATGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACCATTGTaca                                                                                    <  1:1870312/91‑1 (MQ=255)
        cGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACCATTGTaca                                                                                    <  1:1870539/91‑1 (MQ=255)
          cTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACAcc                                                                                  <  1:3644909/91‑1 (MQ=255)
            ccATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCAt                                                                                >  1:2672203/1‑91 (MQ=255)
             cATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACCATTGTacactgtt                                                                               <  1:2062948/91‑5 (MQ=255)
               tCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTtc                                                                             <  1:612763/91‑1 (MQ=255)
               tCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTtc                                                                             <  1:3182717/91‑1 (MQ=255)
                 gATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTtctc                                                                           <  1:2310693/91‑1 (MQ=255)
                  atatatAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCa                                                                          <  1:2258896/91‑1 (MQ=255)
                   tatataAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCAt                                                                         >  1:3604089/1‑91 (MQ=255)
                     tataAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGa                                                                       <  1:1163458/91‑1 (MQ=255)
                         aaGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGAtttt                                                                   >  1:2943114/1‑91 (MQ=255)
                           ggTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAg                                                                 >  1:1428492/1‑91 (MQ=255)
                           ggTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAg                                                                 >  1:2367219/1‑91 (MQ=255)
                             tAAAATTACGTTTAATAAAGCACCAATTAAAATACCAGGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTa                                                               <  1:2518654/91‑1 (MQ=255)
                             tAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTa                                                               >  1:557337/1‑91 (MQ=255)
                              aaaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaa                                                              >  1:2818265/1‑91 (MQ=255)
                              aaaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaa                                                              >  1:2818095/1‑91 (MQ=255)
                              aaaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaa                                                              >  1:2817812/1‑91 (MQ=255)
                              aaaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaa                                                              >  1:2817647/1‑91 (MQ=255)
                               aaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaaa                                                             <  1:3938682/91‑1 (MQ=255)
                               aaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaaa                                                             >  1:3811811/1‑91 (MQ=255)
                                aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                            <  1:2706908/91‑1 (MQ=255)
                                aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                            <  1:301035/91‑1 (MQ=255)
                                aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                            <  1:299475/91‑1 (MQ=255)
                                aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                            >  1:2044651/1‑91 (MQ=255)
                                aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                            >  1:4136539/1‑91 (MQ=255)
                                     cGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATAtttt                                                       >  1:3808767/1‑91 (MQ=255)
                                       tttAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCg                                                     <  1:2391670/91‑1 (MQ=255)
                                        ttAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCgg                                                    <  1:3571127/91‑1 (MQ=255)
                                           ataaAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAt                                                 >  1:3845726/1‑91 (MQ=255)
                                            taaAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAtt                                                <  1:1654492/91‑1 (MQ=255)
                                             aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                               <  1:513169/91‑1 (MQ=255)
                                             aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                               <  1:4154135/91‑1 (MQ=255)
                                             aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                               <  1:4154013/91‑1 (MQ=255)
                                             aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCACGATTTTAGTAAATATTTTCGGCAttt                                               <  1:489116/91‑1 (MQ=255)
                                              aaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTc                                              >  1:1864185/1‑91 (MQ=255)
                                              aaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATAGTTTCGGCATTTc                                              <  1:3188189/91‑1 (MQ=255)
                                               aGCACCAATTAAAATACTAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCa                                             <  1:2769353/91‑1 (MQ=255)
                                               aGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCa                                             <  1:2710865/91‑1 (MQ=255)
                                                  aCCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCt                                          <  1:70996/91‑1 (MQ=255)
                                                     aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                       >  1:684238/1‑91 (MQ=255)
                                                     aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                       >  1:688268/1‑91 (MQ=255)
                                                     aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                       <  1:787609/91‑1 (MQ=255)
                                                     aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                       <  1:788690/91‑1 (MQ=255)
                                                      aTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTAt                                      <  1:2570614/91‑1 (MQ=255)
                                                      aTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTAt                                      <  1:2570544/91‑1 (MQ=255)
                                                       ttAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATg                                     <  1:1398346/91‑1 (MQ=255)
                                                        tAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGa                                    <  1:2490122/91‑1 (MQ=255)
                                                        tAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGa                                    <  1:3317641/91‑1 (MQ=255)
                                                         aaaaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAg                                   >  1:3789970/1‑91 (MQ=255)
                                                         aaaaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAg                                   >  1:3789920/1‑91 (MQ=255)
                                                          aaaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGc                                  >  1:2696489/1‑91 (MQ=255)
                                                           aaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCa                                 <  1:3063976/91‑1 (MQ=255)
                                                           aaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCa                                 <  1:3063946/91‑1 (MQ=255)
                                                           aaTACCAAGAGATACTGCTACAACAGTGTACAATGATACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCa                                 <  1:1187166/91‑1 (MQ=255)
                                                             tACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAg                               >  1:2113993/1‑91 (MQ=255)
                                                                    agaTACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAACTATTTTCGGCATTTCACCTTTATGAGCAAGTGTATaa                        <  1:3498808/91‑1 (MQ=255)
                                                                      aTACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAAt                      >  1:1385600/1‑91 (MQ=255)
                                                                        aCTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCa                    <  1:502905/91‑1 (MQ=255)
                                                                        aCTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCa                    <  1:504039/91‑1 (MQ=255)
                                                                        aCTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCa                    <  1:3350655/91‑1 (MQ=255)
                                                                         cTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCAt                   >  1:1559083/1‑91 (MQ=255)
                                                                           gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                 <  1:4074028/91‑1 (MQ=255)
                                                                           gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                 >  1:687739/1‑91 (MQ=255)
                                                                           gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                 >  1:270994/1‑91 (MQ=255)
                                                                               caacaGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACt             <  1:2757744/91‑1 (MQ=255)
                                                                                   ggtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa         >  1:3380112/2‑91 (MQ=255)
                                                                                   aGTGTACAATGGTACACCATTTCTCATGATTTTAGTACATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa         >  1:3484284/1‑91 (MQ=255)
                                                                                   aGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa         >  1:208410/1‑91 (MQ=255)
                                                                                    gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc        >  1:1887136/1‑91 (MQ=255)
                                                                                    gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc        >  1:1880810/1‑91 (MQ=255)
                                                                                    gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc        >  1:1880733/1‑91 (MQ=255)
                                                                                    gtgtACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc        <  1:2062863/91‑1 (MQ=255)
                                                                                    gtgtACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc        <  1:2706257/91‑1 (MQ=255)
                                                                                     tgtgcAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACt       >  1:54563/5‑91 (MQ=255)
                                                                                     tgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACt       >  1:3044062/1‑91 (MQ=255)
                                                                                      gtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTg      >  1:3114841/1‑91 (MQ=255)
                                                                                      gtACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTg      <  1:1701856/91‑1 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa    <  1:3776541/91‑1 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa    >  1:2544874/1‑91 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa    >  1:237904/1‑91 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa    >  1:237847/1‑91 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa    >  1:3811143/1‑91 (MQ=255)
                                                                                        aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGAAAGTGTATAAATCATACGACTCGCACTGaa    <  1:1596305/91‑1 (MQ=255)
                                                                                         cAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaa   >  1:2368259/1‑91 (MQ=255)
                                                                                          aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa  >  1:697827/1‑91 (MQ=255)
                                                                                          aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa  >  1:697981/1‑91 (MQ=255)
                                                                                          aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa  >  1:699842/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATACTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAA  >  USA300TCH1516_ALE/2435417‑2435597

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