Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,352,660 A→G 7.2% intergenic (‑514/+31) USA300HOU_RS12000 ← / ← rpsI hypothetical protein/30S ribosomal protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,352,6600AG 7.2% 311.6 / 11.0 99intergenic (‑514/+31)USA300HOU_RS12000/rpsIhypothetical protein/30S ribosomal protein S9
Reads supporting (aligned to +/- strand):  ref base A (32/59);  new base G (2/5);  total (34/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 109 total aligned reads displayed.

CATCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGT  >  USA300TCH1516_ALE/2352571‑2352748
                                                                                         |                                                                                        
cATCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAatat                                                                                         >  1:2760083/1‑91 (MQ=255)
 aTCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATc                                                                                        <  1:3049449/91‑1 (MQ=255)
 aTCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATc                                                                                        >  1:2029042/1‑91 (MQ=255)
  tCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCg                                                                                       <  1:922832/91‑1 (MQ=255)
  tCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCg                                                                                       <  1:3802416/91‑1 (MQ=255)
   cTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGa                                                                                      <  1:592130/91‑1 (MQ=255)
    ttttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAg                                                                                     <  1:2940999/91‑1 (MQ=255)
    ttttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAg                                                                                     <  1:3697098/91‑1 (MQ=255)
     tttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGATATTATGTCGAgg                                                                                    <  1:1390037/91‑1 (MQ=38)
     tttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAgg                                                                                    <  1:1564703/91‑1 (MQ=255)
     tttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAgg                                                                                    <  1:1564521/91‑1 (MQ=255)
     tttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAgg                                                                                    <  1:1097326/91‑1 (MQ=255)
      ttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGATATTATGTCGAGgt                                                                                   <  1:1919250/91‑1 (MQ=38)
      ttATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGgt                                                                                   >  1:1432404/1‑91 (MQ=255)
       tATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGATATTATGTCGAGgtg                                                                                  <  1:2947394/91‑1 (MQ=38)
       tATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGgtg                                                                                  >  1:670351/1‑91 (MQ=255)
        aTGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGgtgt                                                                                 <  1:747464/91‑1 (MQ=255)
        aTGGGGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGgtgt                                                                                 <  1:3476304/91‑1 (MQ=255)
         tggtggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGATATTATGTCGAGGTGtt                                                                                <  1:2401845/91‑1 (MQ=38)
         tggtggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGATATTATGTCGAGGTGtt                                                                                <  1:2401447/91‑1 (MQ=38)
         tggtggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACCTAATATCGAGGTGtt                                                                                >  1:1352438/1‑91 (MQ=255)
         tggtggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGtt                                                                                >  1:2256535/1‑91 (MQ=255)
          ggtggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGttt                                                                               >  1:1262277/1‑91 (MQ=255)
            tggGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTg                                                                             >  1:870672/1‑91 (MQ=255)
             gggAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGt                                                                            >  1:307766/1‑91 (MQ=255)
             gggAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGt                                                                            >  1:121945/1‑91 (MQ=255)
              ggAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACACAAAAAAACACCTCGACATAATATCGAGGTGTTTTGta                                                                           >  1:722999/1‑91 (MQ=255)
              ggAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACACAAAAAAACACCTCGACATAATATCGAGGTGTTTTGta                                                                           >  1:723524/1‑91 (MQ=255)
                 aGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGtatat                                                                        <  1:3587714/91‑1 (MQ=255)
                 aGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGtatat                                                                        <  1:3587470/91‑1 (MQ=255)
                 aGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGtatat                                                                        <  1:3774194/91‑1 (MQ=255)
                  ggTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGtatata                                                                       >  1:1632556/1‑91 (MQ=255)
                   gTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTatatat                                                                      >  1:1569302/1‑91 (MQ=255)
                   gTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTatatat                                                                      <  1:1995255/91‑1 (MQ=255)
                   gTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTatatat                                                                      >  1:701828/1‑91 (MQ=255)
                      aaaaTTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGt                                                                   <  1:4136311/91‑1 (MQ=255)
                      aaaaTTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGt                                                                   <  1:4135820/91‑1 (MQ=255)
                         aTTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCg                                                                <  1:1991478/91‑1 (MQ=255)
                            cATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACa                                                             >  1:1029674/1‑91 (MQ=255)
                            cATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACa                                                             >  1:3696199/1‑91 (MQ=255)
                              tttATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAAt                                                           >  1:2670204/1‑91 (MQ=255)
                               ttATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACTTAATATCGAGGTGTTTTGTATATATCGTCCGACAAtt                                                          <  1:3505304/91‑1 (MQ=255)
                               ttATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACGTAATATCGAGGTGTTTTGTATATATCGTCCGACAAtt                                                          <  1:3505617/91‑1 (MQ=255)
                               ttATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAAtt                                                          <  1:1081314/91‑1 (MQ=255)
                               ttATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAAtt                                                          >  1:1184349/1‑91 (MQ=255)
                                tataTTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTa                                                         <  1:3156467/91‑1 (MQ=255)
                                tataTTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTa                                                         <  1:1193683/91‑1 (MQ=255)
                                tataTGCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTa                                                         <  1:3156282/91‑1 (MQ=255)
                                 atatTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTaa                                                        >  1:1190735/1‑91 (MQ=255)
                                  tatTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAAc                                                       <  1:2514375/91‑1 (MQ=255)
                                  tatTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAAc                                                       <  1:1747312/91‑1 (MQ=255)
                                   atTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACg                                                      <  1:3428995/91‑1 (MQ=255)
                                    ttCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGt                                                     <  1:3216319/91‑1 (MQ=255)
                                     tCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGtt                                                    <  1:3664113/91‑1 (MQ=255)
                                     tCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGtt                                                    <  1:1467775/91‑1 (MQ=255)
                                     tCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGtt                                                    <  1:783714/91‑1 (MQ=255)
                                      cTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGttt                                                   <  1:2909928/91‑1 (MQ=255)
                                        acacATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGTGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTg                                                 <  1:3499507/91‑1 (MQ=255)
                                        acacATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTg                                                 <  1:2249775/91‑1 (MQ=255)
                                        acacATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTg                                                 <  1:1510606/91‑1 (MQ=255)
                                        acacATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTg                                                 <  1:2249788/91‑1 (MQ=255)
                                        acacATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGCGTTTTGTATATATCGTCCGACAATTAACGTTTTg                                                 >  1:1001003/1‑91 (MQ=255)
                                         cacaTTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTga                                                <  1:2918134/91‑1 (MQ=255)
                                         cacaTTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTga                                                <  1:320235/91‑1 (MQ=255)
                                             ttCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAAt                                            >  1:3442628/1‑91 (MQ=255)
                                                 atatTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAg                                        >  1:1734412/1‑91 (MQ=255)
                                                     tCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGa                                    <  1:1369441/91‑1 (MQ=255)
                                                     tCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGa                                    <  1:1369382/91‑1 (MQ=255)
                                                        ccgcAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAacg                                 >  1:438156/1‑91 (MQ=255)
                                                          gcaaaaacgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAacgac                               <  1:254898/91‑1 (MQ=255)
                                                            aaaaacgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAAcgacga                             <  1:3362767/91‑1 (MQ=255)
                                                            aaaaacgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAAcgacga                             >  1:1998189/1‑91 (MQ=255)
                                                             aaaacgcgaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAg                            <  1:2098592/91‑1 (MQ=255)
                                                             aaaacgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAg                            <  1:262204/91‑1 (MQ=255)
                                                             aaaacgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAg                            <  1:3257374/91‑1 (MQ=255)
                                                                acgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAgctg                         >  1:2194328/1‑91 (MQ=255)
                                                                acgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAgctg                         <  1:1515511/91‑1 (MQ=255)
                                                                 cgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAgctgc                        <  1:1166592/91‑1 (MQ=255)
                                                                 cgcaaaaaAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAgctgc                        >  1:2522206/1‑91 (MQ=255)
                                                                   atacaaaaCACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCtt                      <  1:822944/87‑1 (MQ=255)
                                                                     aaaaaaCACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTa                    <  1:3923123/91‑1 (MQ=255)
                                                                         aaCACCTCGATATTATGTCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAc                >  1:4004451/1‑91 (MQ=14)
                                                                         aaCACCTCGATATTATGTCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAc                >  1:4004660/1‑91 (MQ=14)
                                                                          acacCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAcc               <  1:2844011/91‑1 (MQ=255)
                                                                          acacCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAcc               <  1:794461/91‑1 (MQ=255)
                                                                          acacCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAcc               <  1:124454/91‑1 (MQ=255)
                                                                          acacCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAcc               <  1:3567396/91‑1 (MQ=255)
                                                                          acacCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGAcc               <  1:3569321/91‑1 (MQ=255)
                                                                             ccTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTgg            >  1:1749172/1‑91 (MQ=255)
                                                                               tCGACATAATTTCGAGGTGTTTTGTATATATCGACCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGtt          <  1:2923970/91‑1 (MQ=255)
                                                                               tCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGtt          >  1:1743034/1‑91 (MQ=255)
                                                                               tCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGtt          <  1:3474395/91‑1 (MQ=255)
                                                                               tCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGtt          <  1:918898/91‑1 (MQ=255)
                                                                                    ataataTCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTa     <  1:3099359/91‑1 (MQ=255)
                                                                                    ataataTCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTa     <  1:3435729/91‑1 (MQ=255)
                                                                                       atatCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGt  <  1:192407/91‑1 (MQ=255)
                                                                                       atatCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGt  >  1:4093559/1‑91 (MQ=255)
                                                                                       atatCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGt  <  1:3248098/91‑1 (MQ=255)
                                                                                       atatCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGt  >  1:2478122/1‑91 (MQ=255)
                                                                                       atatCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGt  >  1:2478121/1‑91 (MQ=255)
                                                                                         |                                                                                        
CATCTTTTATGGTGGGAAGGTAAAAATTCATTTATATTCTACACATTCCATATTCGCCGCAAAAACGCAAAAAAACACCTCGACATAATATCGAGGTGTTTTGTATATATCGTCCGACAATTAACGTTTTGAGAATTGAGGTGAACGACGAGCTGCTTTAAGACCTGGTTTTTTACGT  >  USA300TCH1516_ALE/2352571‑2352748

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