Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,378,337 G→T 6.4% intergenic (+196/+223) glcU_2 → / ← USA300HOU_RS12220 putative glucose uptake protein GlcU/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,378,3370GT 6.4% 386.5 / 11.7 110intergenic (+196/+223)glcU_2/USA300HOU_RS12220putative glucose uptake protein GlcU/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (48/54);  new base T (5/3);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CTCTATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTATCAA  >  USA300TCH1516_ALE/2378247‑2378426
                                                                                          |                                                                                         
ctctATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATATCACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAg                                                                                           >  1:4153299/1‑91 (MQ=255)
ctctATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAg                                                                                           <  1:137338/91‑1 (MQ=255)
ctctATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAg                                                                                           <  1:141646/91‑1 (MQ=255)
 tctATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGc                                                                                          <  1:1751371/91‑1 (MQ=255)
   tataTGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGa                                                                                        <  1:3261788/91‑1 (MQ=255)
    atatGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGAt                                                                                       >  1:3381210/1‑91 (MQ=255)
      atGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATAc                                                                                     <  1:2373061/91‑1 (MQ=255)
      atGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATAc                                                                                     <  1:2224614/91‑1 (MQ=255)
         tCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACttt                                                                                  <  1:3931971/91‑1 (MQ=255)
         tCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACttt                                                                                  <  1:2913135/91‑1 (MQ=255)
         tCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACtag                                                                                  >  1:554945/1‑89 (MQ=255)
          ccTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACtttt                                                                                 <  1:114331/91‑1 (MQ=255)
          ccTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACtttt                                                                                 <  1:113599/91‑1 (MQ=255)
            tAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaa                                                                               <  1:788127/91‑1 (MQ=255)
            tAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaa                                                                               <  1:706924/91‑1 (MQ=255)
            tAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaa                                                                               <  1:1932311/91‑1 (MQ=255)
             aGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaaa                                                                              >  1:1401647/1‑91 (MQ=255)
              ggAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaaaa                                                                             >  1:4120114/1‑91 (MQ=255)
              ggAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaaaa                                                                             >  1:3958145/1‑91 (MQ=255)
              ggAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTaaaa                                                                             >  1:2273166/1‑91 (MQ=255)
               gAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAAt                                                                            >  1:620508/1‑91 (MQ=255)
                 aaTACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCt                                                                          >  1:516395/1‑91 (MQ=255)
                 aaTACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCt                                                                          <  1:2329123/91‑1 (MQ=255)
                  aTACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGATTTTAAAAGTATCGCTAATTTTAAAGTCtt                                                                         <  1:2615523/91‑1 (MQ=21)
                  aTACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCtt                                                                         >  1:623355/1‑91 (MQ=255)
                   tACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtg                                                                        >  1:2450101/1‑91 (MQ=255)
                   tACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtg                                                                        >  1:3126458/1‑91 (MQ=255)
                   tACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtg                                                                        >  1:3426299/1‑91 (MQ=255)
                    aCGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtga                                                                       >  1:2419425/1‑91 (MQ=255)
                     cGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtgat                                                                      >  1:1021462/1‑91 (MQ=255)
                     cGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtgat                                                                      <  1:1059832/91‑1 (MQ=255)
                      gTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtgatg                                                                     >  1:3235914/1‑91 (MQ=255)
                       tGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTtgatga                                                                    <  1:3897614/91‑1 (MQ=255)
                         gCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGata                                                                  <  1:4092058/91‑1 (MQ=255)
                           tctAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATAtt                                                                <  1:1478289/91‑1 (MQ=255)
                            ctAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTc                                                               <  1:59692/91‑1 (MQ=255)
                            ctAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTc                                                               <  1:61635/91‑1 (MQ=255)
                            ctAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTc                                                               >  1:2828475/1‑91 (MQ=255)
                            ctAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTc                                                               <  1:2874690/91‑1 (MQ=255)
                               aaaacaaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCata                                                            >  1:771059/1‑91 (MQ=255)
                               aaaacaaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCata                                                            >  1:3042251/1‑91 (MQ=255)
                                 aacaaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCatata                                                          <  1:1928819/91‑1 (MQ=255)
                                    aaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCa                                                       >  1:3864249/1‑91 (MQ=255)
                                    aaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCa                                                       <  1:1892305/91‑1 (MQ=255)
                                    aaaacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCa                                                       <  1:668092/91‑1 (MQ=255)
                                      aacGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAg                                                     >  1:2419358/1‑91 (MQ=255)
                                        cGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTa                                                   <  1:2957977/91‑1 (MQ=255)
                                        cGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTa                                                   >  1:2952900/1‑91 (MQ=255)
                                        cGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTa                                                   >  1:2821421/1‑91 (MQ=255)
                                           aaTAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgc                                                >  1:3606725/1‑91 (MQ=255)
                                           aaTAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgc                                                >  1:3606824/1‑91 (MQ=255)
                                           aaTAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgc                                                >  1:2862737/1‑91 (MQ=255)
                                           aaTAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgc                                                >  1:2863033/1‑91 (MQ=255)
                                            aTAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgcg                                               >  1:2192517/1‑91 (MQ=255)
                                             tAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATgcgc                                              >  1:821722/1‑91 (MQ=255)
                                              aaCACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGcc                                             <  1:2637946/91‑1 (MQ=255)
                                               acacATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCca                                            >  1:2809810/1‑91 (MQ=255)
                                                  catcatGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCcatac                                         >  1:1063934/1‑91 (MQ=255)
                                                   atcatGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCcataca                                        >  1:565458/1‑91 (MQ=255)
                                                    tcatGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCcatacat                                       >  1:811497/1‑91 (MQ=255)
                                                     catGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCcatacata                                      <  1:361301/91‑1 (MQ=255)
                                                       tGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATg                                    >  1:1649366/1‑91 (MQ=255)
                                                          cATTAATGATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 <  1:2082327/91‑1 (MQ=255)
                                                          cATTAATCATGCGTTTTAAGACTTTGAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 <  1:3847548/91‑1 (MQ=255)
                                                          cATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTATAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 <  1:3295627/91‑1 (MQ=255)
                                                          cATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 <  1:474237/91‑1 (MQ=255)
                                                          cATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 >  1:3546533/1‑91 (MQ=255)
                                                          cATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAg                                 <  1:355642/91‑1 (MQ=255)
                                                            ttAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTg                               >  1:2730428/1‑91 (MQ=255)
                                                              aaTCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGa                             >  1:3191369/1‑91 (MQ=255)
                                                              aaTCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGa                             >  1:1017538/1‑91 (MQ=255)
                                                              aaTCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGa                             <  1:1873776/91‑1 (MQ=255)
                                                              aaTCATGCGTTTTAAGACTTTAAAAGTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGa                             <  1:3864340/91‑1 (MQ=255)
                                                               aTCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGAt                            <  1:1777953/91‑1 (MQ=255)
                                                                    gCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTg                       >  1:2689042/1‑91 (MQ=255)
                                                                    gCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTg                       >  1:3437979/1‑91 (MQ=255)
                                                                      gTTGTAAGAGTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCa                     <  1:3728668/91‑1 (MQ=255)
                                                                       ttttAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCAt                    <  1:3235184/91‑1 (MQ=255)
                                                                       ttttAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCAt                    <  1:3478471/91‑1 (MQ=255)
                                                                       catcaaGATTTTAAAAGTATCGCTAATTTTAAAGTCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCAt                    >  1:918406/5‑91 (MQ=18)
                                                                        tttAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATa                   <  1:1502993/91‑1 (MQ=255)
                                                                          tAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAAc                 <  1:3160301/91‑1 (MQ=255)
                                                                          tAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAAc                 <  1:3158047/91‑1 (MQ=255)
                                                                           aaGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACg                <  1:1376976/91‑1 (MQ=255)
                                                                           aaGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACg                <  1:1542919/91‑1 (MQ=255)
                                                                             gattttAAAAGTATCGCTAATTTTAAAGTCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCa              >  1:816640/4‑91 (MQ=25)
                                                                             gattttAAAAGTATCGCTAATTTTAAAGTCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCa              >  1:816713/4‑91 (MQ=25)
                                                                             gACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCa              >  1:1561916/1‑91 (MQ=255)
                                                                             gACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCa              >  1:1188111/1‑91 (MQ=255)
                                                                              aCTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCAc             <  1:691936/91‑1 (MQ=255)
                                                                              aCTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCAc             <  1:692108/91‑1 (MQ=255)
                                                                                tttAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTg           <  1:528932/91‑1 (MQ=255)
                                                                                  tAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCa         <  1:307285/91‑1 (MQ=255)
                                                                                   aaaaTTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCAt        <  1:1230295/91‑1 (MQ=255)
                                                                                    aaaTTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTgagat       >  1:2385335/1‑87 (MQ=255)
                                                                                      aTTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTTCATTAt     <  1:3231739/91‑1 (MQ=255)
                                                                                      aTTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTAt     <  1:2530974/91‑1 (MQ=255)
                                                                                       ttATCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTATc    <  1:3337000/91‑1 (MQ=255)
                                                                                       ttATCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTATc    <  1:3337261/91‑1 (MQ=255)
                                                                                         aGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTATCaa  >  1:2001405/1‑91 (MQ=255)
                                                                                          |                                                                                         
CTCTATATGTCCTAGGAAATACGTGGCTCTAAAAACAAAACGCAATAACACATCATGACATTAATCATGCGTTTTAAGACTTTAAAATTAGCGATACTTTTAAAATCTTGATGATATTCATATATCAAGTATGCGCCATACATATGAAGTGGATAGCTGCATAACGCACTGCATTATCAA  >  USA300TCH1516_ALE/2378247‑2378426

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