Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 482,712 C→A 38.9% intergenic (+509/‑174) USA300HOU_RS02285 → / → USA300HOU_RS02295 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE482,7120CA38.9% 256.9 / 306.2 252intergenic (+509/‑174)USA300HOU_RS02285/USA300HOU_RS02295hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (83/71);  new base A (54/44);  total (137/115)
Fisher's exact test for biased strand distribution p-value = 8.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 260 total aligned reads displayed.

AATCTAACTAATGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTGATAT  >  USA300TCH1516_ALE/482622‑482802
                                                                                          |                                                                                          
aaTCTAACTAATGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAc                                                                                            <  2:3106929/91‑1 (MQ=255)
aaTCTAACTAATGTAACGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAc                                                                                            <  2:3036787/91‑1 (MQ=255)
 ctCTAACTAATGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACt                                                                                           <  2:3926434/90‑1 (MQ=255)
   ctaactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGtaataa                                                                                         <  2:1976306/91‑1 (MQ=255)
   ctaactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTaa                                                                                         >  2:1004063/1‑91 (MQ=255)
   ctaactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTaa                                                                                         >  2:3954150/1‑91 (MQ=255)
   ctaactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTaa                                                                                         <  1:2085062/91‑1 (MQ=255)
    taactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTaaa                                                                                        <  1:547957/91‑1 (MQ=255)
    taactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTaaa                                                                                        <  1:2911263/91‑1 (MQ=255)
     aactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAAt                                                                                       >  1:3441406/1‑91 (MQ=255)
     aactaaTGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAAt                                                                                       >  1:2237297/1‑91 (MQ=255)
        ttatgtaattAGTAACTCTATATAGACAAAAATATCTCTATTGAAACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGa                                                                                    <  2:3167508/89‑1 (MQ=255)
        taatgtaatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGa                                                                                    <  2:2264906/91‑1 (MQ=255)
           tgtaatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGaaaa                                                                                 <  2:2291816/91‑1 (MQ=255)
            gtaatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAAt                                                                                <  1:1641716/91‑1 (MQ=255)
            gtaatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAAt                                                                                >  2:2402241/1‑91 (MQ=255)
              aatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAg                                                                              <  1:945714/91‑1 (MQ=255)
              aatgATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAg                                                                              >  1:2604268/1‑91 (MQ=255)
                   ttAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCaaaa                                                                         >  1:1014340/1‑91 (MQ=255)
                       cGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATAtt                                                                     <  2:1826954/91‑1 (MQ=255)
                       cGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATAtt                                                                     >  2:411851/1‑91 (MQ=255)
                         cTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGt                                                                   <  2:1878842/91‑1 (MQ=255)
                          tatataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGtt                                                                  <  1:2423298/91‑1 (MQ=255)
                          tatataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGtt                                                                  >  2:2900708/1‑91 (MQ=255)
                          tatataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGtt                                                                  >  2:3905041/1‑91 (MQ=255)
                           atataGACAAAACTAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATAttatta                                                                 >  1:3057017/1‑91 (MQ=255)
                           atataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTa                                                                 >  2:961557/1‑91 (MQ=255)
                           atataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTa                                                                 >  2:351651/1‑91 (MQ=255)
                           atataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTa                                                                 >  2:86786/1‑91 (MQ=255)
                           atataGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTa                                                                 >  1:1340375/1‑91 (MQ=255)
                                gACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCa                                                            >  1:2576176/1‑91 (MQ=255)
                                  cAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAAt                                                          <  1:2193906/91‑1 (MQ=255)
                                  cAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAAt                                                          <  2:518467/91‑1 (MQ=255)
                                   aaaaaTAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAAtt                                                         <  2:2633372/91‑1 (MQ=255)
                                     aaaTAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATtat                                                       >  2:1508585/1‑91 (MQ=255)
                                       aTAACTCATTTGATACGCCTTACCCCCTCTCTATAATAAATATGCATCGTACTAAATTGAAGATAGCAAATTATTGTGATTTCAATTATAc                                                     <  2:1587112/91‑1 (MQ=255)
                                        tAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACa                                                    >  2:2199898/1‑91 (MQ=255)
                                        tAACTCTATTGATACGCCTTACCCCCTCTGTAGAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACa                                                    >  1:603571/1‑91 (MQ=255)
                                         aaCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACaa                                                   >  2:632142/1‑91 (MQ=255)
                                         aaCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATACAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACaa                                                   >  2:2590955/1‑91 (MQ=255)
                                          aCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAAt                                                  >  2:2397773/1‑91 (MQ=255)
                                          aCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAAt                                                  >  2:2063021/1‑91 (MQ=255)
                                          aCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAAt                                                  >  2:1743170/1‑91 (MQ=255)
                                          aCTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAAt                                                  >  2:1471973/1‑91 (MQ=255)
                                            tctATTGATACGCCTTACCCCCTCTGTATTATAAATATAACTCGTACTAAATTGAAGATAGCAAAATATTGTTATTTCAATTATACAATGt                                                <  2:1786598/91‑1 (MQ=255)
                                            tctATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGt                                                <  1:2180949/91‑1 (MQ=255)
                                            tctATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGt                                                >  1:673215/1‑91 (MQ=255)
                                             ctATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGtt                                               <  1:764438/91‑1 (MQ=255)
                                              tATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGttt                                              <  2:2530444/91‑1 (MQ=255)
                                               aTTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGtttt                                             <  2:2084499/91‑1 (MQ=255)
                                                ttGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATAGTTATTTCACTTATACAATGttttt                                            >  2:484566/1‑90 (MQ=255)
                                                ttGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTa                                            >  1:1272623/1‑91 (MQ=255)
                                                 tGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTAt                                           <  1:394361/91‑1 (MQ=255)
                                                 tGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTAt                                           <  1:3404230/91‑1 (MQ=255)
                                                  gATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTAtt                                          <  1:428511/91‑1 (MQ=255)
                                                   aTACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTAttt                                         <  1:3800080/91‑1 (MQ=255)
                                                   aTACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTAttt                                         <  1:2995847/91‑1 (MQ=255)
                                                    tACGCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTg                                        >  1:1819221/1‑91 (MQ=255)
                                                       gCCTTACCCCCTCTGTATAATAAATATAATTCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCaa                                     <  2:1612471/91‑1 (MQ=255)
                                                       gCCTTACCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCaa                                     <  2:671434/91‑1 (MQ=255)
                                                        ccTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAat                                    >  1:3506364/1‑91 (MQ=255)
                                                         cTTACCGCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAata                                   >  1:3075059/1‑91 (MQ=255)
                                                         cTTACCCCCTCTGTATAATACATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAata                                   >  1:1410921/1‑91 (MQ=255)
                                                         cTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAata                                   <  1:2265649/91‑1 (MQ=255)
                                                           tACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAAtata                                 <  2:1672784/91‑1 (MQ=255)
                                                            aCCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATAc                                >  2:1886425/1‑91 (MQ=255)
                                                            aCCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATAc                                >  1:95746/1‑91 (MQ=255)
                                                            aCCCCCTCTGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATAc                                >  2:1075204/1‑91 (MQ=255)
                                                             cccccTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACa                               <  1:1980078/91‑1 (MQ=255)
                                                              ccccTCTGTATAATAAATATACATCGTACTAAATTGACAATAGCACAATATCGTTATTTTAATTATACAATGTTTTATTTGCAATATGCAt                              >  2:2627209/1‑91 (MQ=255)
                                                               cccTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATa                             <  2:1701714/91‑1 (MQ=255)
                                                                 ctctGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAAc                           >  2:478882/1‑91 (MQ=255)
                                                                 ctctGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAAc                           >  1:2970358/1‑91 (MQ=255)
                                                                 ctctGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAAc                           <  2:3841040/91‑1 (MQ=255)
                                                                 ctctGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAAc                           <  2:1180015/91‑1 (MQ=255)
                                                                  tctGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACa                          <  1:3625524/91‑1 (MQ=255)
                                                                  tctGTATAATAAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACa                          >  1:1169757/1‑91 (MQ=255)
                                                                   ctGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACaa                         >  1:3404888/1‑91 (MQ=255)
                                                                       ataataAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATg                     >  2:2745929/1‑91 (MQ=255)
                                                                       ataataAATATAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATg                     >  1:2846198/1‑91 (MQ=255)
                                                                         aataagtaTAAATCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGaa                   <  2:3646057/91‑1 (MQ=255)
                                                                          ataaatataaatCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGaaa                  >  2:3272119/1‑90 (MQ=255)
                                                                          ataaatataaatCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAAc                  >  2:761420/1‑91 (MQ=255)
                                                                           taaatataaatCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACt                 <  1:2737361/91‑1 (MQ=255)
                                                                            aaatataaatCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACtt                <  2:3280980/91‑1 (MQ=255)
                                                                             aatataaatCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACttt               <  2:3749012/91‑1 (MQ=255)
                                                                              atataaatCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTta              >  2:3688159/1‑91 (MQ=255)
                                                                               tataaatCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTtat             >  2:2119207/1‑91 (MQ=255)
                                                                               tataaatCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTtat             >  1:2532005/1‑91 (MQ=255)
                                                                                  aaatCGTACTAAATTGAAACTAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAAt          >  2:3461247/1‑91 (MQ=255)
                                                                                    atCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAAttt        <  2:2955580/91‑1 (MQ=255)
                                                                                    atCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAAttt        <  2:2955849/91‑1 (MQ=255)
                                                                                    atCGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAAttt        <  2:1290456/91‑1 (MQ=255)
                                                                                      cGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTg      >  2:2230899/1‑91 (MQ=255)
                                                                                      cGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTg      >  2:145527/1‑91 (MQ=255)
                                                                                      cGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTg      <  1:498215/91‑1 (MQ=255)
                                                                                      cGTAATAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTg      <  1:2381713/91‑1 (MQ=255)
                                                                                       gTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTGa     >  1:1775503/1‑91 (MQ=255)
                                                                                        tACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTGat    >  1:1504172/1‑91 (MQ=255)
                                                                                        tACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTGat    >  1:813262/1‑91 (MQ=255)
                                                                                          |                                                                                          
AATCTAACTAATGTAATGATTAACGCTATATAGACAAAAATAACTCTATTGATACGCCTTACCCCCTCTGTATAATAAATATAAATCGTACTAAATTGAAAATAGCAAAATATTGTTATTTCAATTATACAATGTTTTATTTGCAATATACATAACAATATGAACTTTATAATTTTGATAT  >  USA300TCH1516_ALE/482622‑482802

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: