Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,435,510 G→T 9.0% P330Q (CCA→CAA)  yifK ← putative transport protein YifK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,435,5100GT 9.0% 359.9 / 19.7 108P330Q (CCA→CAA) yifKputative transport protein YifK
Reads supporting (aligned to +/- strand):  ref base G (43/55);  new base T (7/3);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

CTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAATAC  >  USA300TCH1516_ALE/2435420‑2435600
                                                                                          |                                                                                          
ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                            >  1:333922/1‑91 (MQ=255)
ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                            <  1:1574718/91‑1 (MQ=255)
ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                            <  1:365608/91‑1 (MQ=255)
ctctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATg                                                                                            >  1:1791424/1‑91 (MQ=255)
  ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                          <  1:677051/91‑1 (MQ=255)
  ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                          >  1:1373809/1‑91 (MQ=255)
  ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                          >  1:1923540/1‑91 (MQ=255)
  ctTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGt                                                                                          >  1:3673857/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:3182717/91‑1 (MQ=255)
            tCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTtc                                                                                <  1:612763/91‑1 (MQ=255)
              gATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTtctc                                                                              <  1:2310693/91‑1 (MQ=255)
                tatataAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCAt                                                                            >  1:3604089/1‑91 (MQ=255)
                  tataAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGa                                                                          <  1:1163458/91‑1 (MQ=255)
                   ataAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGAt                                                                         >  1:1856972/1‑91 (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:3811811/1‑91 (MQ=255)
                            aaaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTaaa                                                                <  1:3938682/91‑1 (MQ=255)
                             aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                               >  1:2044651/1‑91 (MQ=255)
                             aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                               <  1:2706908/91‑1 (MQ=255)
                             aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                               >  1:4136539/1‑91 (MQ=255)
                             aaTTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAat                                                               <  1:301035/91‑1 (MQ=255)
                                 aCGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATAttt                                                           <  1:2211592/91‑1 (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:2657890/91‑1 (MQ=255)
                                        ataaAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAt                                                    >  1:3845726/1‑91 (MQ=255)
                                         taaAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAtt                                                   <  1:3180085/91‑1 (MQ=255)
                                         taaAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAtt                                                   <  1:1654492/91‑1 (MQ=255)
                                          aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                                  <  1:4154013/91‑1 (MQ=255)
                                          aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                                  <  1:4154135/91‑1 (MQ=255)
                                          aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCAttt                                                  <  1:513169/91‑1 (MQ=255)
                                          aaaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCACGATTTTAGTAAATATTTTCGGCAttt                                                  <  1:489116/91‑1 (MQ=255)
                                           aaGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATAGTTTCGGCATTTc                                                 <  1:3188189/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:788690/91‑1 (MQ=255)
                                                  aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                          <  1:787609/91‑1 (MQ=255)
                                                  aaTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTa                                          >  1:688268/1‑91 (MQ=255)
                                                   aTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTAt                                         <  1:2570544/91‑1 (MQ=255)
                                                   aTTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTAt                                         <  1:2570614/91‑1 (MQ=255)
                                                    ttAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATg                                        <  1:1398346/91‑1 (MQ=255)
                                                     tAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGa                                       <  1:2490027/91‑1 (MQ=255)
                                                     tAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGa                                       <  1:3317641/91‑1 (MQ=255)
                                                      aaaaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAg                                      >  1:3789920/1‑91 (MQ=255)
                                                      aaaaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAg                                      >  1:3789970/1‑91 (MQ=255)
                                                        aaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCa                                    <  1:3063946/91‑1 (MQ=255)
                                                        aaTACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCa                                    <  1:3063976/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:3350655/91‑1 (MQ=255)
                                                                     aCTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCa                       <  1:504039/91‑1 (MQ=255)
                                                                     aCTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCa                       <  1:502905/91‑1 (MQ=255)
                                                                      cTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCAt                      <  1:2736285/91‑1 (MQ=255)
                                                                        gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                    >  1:687739/1‑91 (MQ=255)
                                                                        gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                    >  1:270994/1‑91 (MQ=255)
                                                                        gCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATAc                    <  1:4074028/91‑1 (MQ=255)
                                                                            caacaGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACt                <  1:2757744/91‑1 (MQ=255)
                                                                                ggtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa            >  1:1449582/2‑91 (MQ=255)
                                                                                ggtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa            >  1:3380112/2‑91 (MQ=255)
                                                                                aGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCa            >  1:208410/1‑91 (MQ=255)
                                                                                 gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc           >  1:1880733/1‑91 (MQ=255)
                                                                                 gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc           >  1:1880810/1‑91 (MQ=255)
                                                                                 gtgtACCATTGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc           >  1:1887136/1‑91 (MQ=255)
                                                                                 gtgtACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc           <  1:2706257/91‑1 (MQ=255)
                                                                                 gtgtACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCAc           <  1:2062863/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:3465687/91‑1 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa       >  1:237847/1‑91 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa       >  1:237904/1‑91 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa       >  1:2544874/1‑91 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa       >  1:3811143/1‑91 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaa       <  1:3776541/91‑1 (MQ=255)
                                                                                     aCAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGAAAGTGTATAAATCATACGACTCGCACTGaa       <  1:1596305/91‑1 (MQ=255)
                                                                                      cAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaa      >  1:2368259/1‑91 (MQ=255)
                                                                                       aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa     <  1:2913743/91‑1 (MQ=255)
                                                                                       aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa     >  1:697827/1‑91 (MQ=255)
                                                                                       aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa     >  1:697981/1‑91 (MQ=255)
                                                                                       aaTGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGaaaa     >  1:699842/1‑91 (MQ=255)
                                                                                          ggTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAATAc  <  1:930985/91‑1 (MQ=255)
                                                                                          ggTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAATAc  <  1:3540156/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTCTTCGCTCCATCGATATATAAAGGTAAAATTACGTTTAATAAAGCACCAATTAAAATACCAAGAGATACTGCTACAACAGTGTACAATGGTACACCATTTCTCATGATTTTAGTAAATATTTTCGGCATTTCACCTTTATGAGCAAGTGTATAAATCATACGACTCGCACTGAAAATAC  >  USA300TCH1516_ALE/2435420‑2435600

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