Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,049,713 G→T 5.1% Q457H (CAG→CAT menD → 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,049,7130GT 5.1% 698.1 / 20.2 196Q457H (CAG→CATmenD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base G (107/79);  new base T (4/6);  total (111/85)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 201 total aligned reads displayed.

TGTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTG  >  USA300TCH1516_ALE/1049623‑1049803
                                                                                          |                                                                                          
tgtgCATAAACGAATAACATTATTGATAGGTGATTTATCCTTTTCTCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAg                                                                                            >  2:2472319/1‑91 (MQ=255)
tgtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAg                                                                                            <  1:2833471/91‑1 (MQ=255)
 gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                           >  1:742160/1‑91 (MQ=255)
 gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                           >  1:547401/1‑91 (MQ=255)
  tgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGac                                                                                          >  2:3847037/1‑90 (MQ=255)
   gCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATg                                                                                         <  2:1348482/91‑1 (MQ=255)
      ataaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAat                                                                                      >  2:1796244/3‑91 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  2:2043436/91‑1 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  2:1556199/91‑1 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  1:3348555/91‑1 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  1:1244731/91‑1 (MQ=255)
         aCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATAtt                                                                                   <  1:2421161/91‑1 (MQ=255)
           gAATGACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAttat                                                                                 <  2:212198/91‑3 (MQ=255)
           gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATTTTCAGATGAATATTGt                                                                                 <  1:710761/91‑1 (MQ=255)
           gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAttat                                                                                 <  1:3863131/91‑3 (MQ=255)
           gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGAt                                                                                          <  2:169369/82‑1 (MQ=255)
           gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGt                                                                                 >  2:1908377/1‑91 (MQ=255)
           gAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGt                                                                                 <  2:1893893/91‑1 (MQ=255)
            aaTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAttatt                                                                                <  1:394456/91‑4 (MQ=255)
            aaTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGta                                                                                >  2:2237739/1‑91 (MQ=255)
            aaTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGta                                                                                <  1:2327089/91‑1 (MQ=255)
            aaTAACATTATTGATAGGTGATTTATCATTTTATCACGATATGAATGGACTATTAATGTCAAAATTAGATAATATTCATCTGAATAttatt                                                                                <  2:3507767/91‑4 (MQ=255)
             aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAttattt                                                                               <  1:3159448/91‑5 (MQ=255)
             aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtat                                                                               >  1:2343184/1‑91 (MQ=255)
                 cATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTattatt                                                                           >  1:444145/1‑91 (MQ=255)
                  attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                          <  2:1733415/91‑1 (MQ=255)
                  attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                          <  2:3765589/91‑1 (MQ=255)
                  attattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTg                                                                          >  1:1027204/1‑91 (MQ=255)
                     attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                       >  2:1353405/1‑91 (MQ=255)
                       tGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaa                                                                     <  1:108564/91‑1 (MQ=255)
                         aTAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACg                                                                   <  2:2707758/91‑1 (MQ=255)
                          tAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGa                                                                  <  2:641906/91‑1 (MQ=255)
                           aGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAt                                                                 >  2:1583545/1‑91 (MQ=255)
                            ggTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtg                                                                >  1:2373677/1‑91 (MQ=255)
                                aTTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggtgg                                                            <  2:844786/91‑1 (MQ=255)
                                 tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                           >  2:2533922/1‑91 (MQ=255)
                                  ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                           <  2:313364/90‑1 (MQ=255)
                                  ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                          >  1:1171624/1‑91 (MQ=255)
                                   tatcattttatcatGATATGAATGGTCTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCgg                                                         <  1:3607781/91‑1 (MQ=255)
                                   tatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCgg                                                         <  1:1819793/91‑1 (MQ=255)
                                       attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACGACGATGGTGGCGGTAtt                                                     >  2:458003/1‑91 (MQ=255)
                                       attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtt                                                     >  1:3655962/1‑91 (MQ=255)
                                       attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtt                                                     >  2:457909/1‑91 (MQ=255)
                                         tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                   >  1:3438257/1‑91 (MQ=255)
                                         tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                   >  1:3438566/1‑91 (MQ=255)
                                         tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                   >  1:838008/1‑91 (MQ=255)
                                          ttatcatGATATGCATGGACTATTAATGTCAAAATTAAATAATGTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAttttt                                                  <  1:3882684/91‑1 (MQ=255)
                                           taaaatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                 >  2:2929832/5‑91 (MQ=255)
                                            atcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTc                                                >  1:407160/1‑91 (MQ=255)
                                             tcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCa                                               >  1:1429592/1‑91 (MQ=255)
                                             tcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCa                                               >  2:3590169/1‑91 (MQ=255)
                                              catcatatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATATTTTCat                                              >  2:3443619/5‑91 (MQ=255)
                                              catGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCat                                              >  2:1218041/1‑91 (MQ=255)
                                              catGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCat                                              <  1:3844115/91‑1 (MQ=255)
                                              catGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCat                                              <  2:3033723/91‑1 (MQ=255)
                                               atGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCata                                             >  2:3405704/1‑91 (MQ=255)
                                                 gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                 <  1:2732766/85‑1 (MQ=255)
                                                 gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAtt                                           <  1:2627192/91‑1 (MQ=255)
                                                 gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAtt                                           <  1:612633/91‑1 (MQ=255)
                                                  atatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAttt                                          >  2:545212/1‑91 (MQ=255)
                                                  atatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAttt                                          >  1:3108231/1‑91 (MQ=255)
                                                   tatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCatataaa                                         <  2:2005902/91‑4 (MQ=255)
                                                   tatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTa                                         >  1:2551725/1‑91 (MQ=255)
                                                    atGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAc                                        >  2:2054527/1‑91 (MQ=255)
                                                    atGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAc                                        <  1:1486504/91‑1 (MQ=255)
                                                          ggACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  <  2:1091551/91‑1 (MQ=255)
                                                           gACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaaa                                 <  1:62278/91‑1 (MQ=255)
                                                           gACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaaa                                 >  1:1387251/1‑91 (MQ=255)
                                                           gACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaaa                                 <  1:226378/91‑1 (MQ=255)
                                                             cTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaag                               >  1:2466400/1‑91 (MQ=255)
                                                              tATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaaga                              >  1:1454527/1‑91 (MQ=255)
                                                              tATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaaga                              <  2:712065/91‑1 (MQ=255)
                                                               aTTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaa                             >  2:2190125/1‑91 (MQ=255)
                                                                ttAATGTCAAAATTAAATAATATTCAGATGAATATTGTCTTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaaa                            >  2:3686647/1‑91 (MQ=255)
                                                                ttAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaaa                            >  2:1210148/1‑91 (MQ=255)
                                                                ttAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaaa                            >  2:2616362/1‑91 (MQ=255)
                                                                     gTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                       >  1:157916/1‑91 (MQ=255)
                                                                     gTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                       >  1:2856641/1‑91 (MQ=255)
                                                                       cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  <  2:2733455/78‑1 (MQ=255)
                                                                       cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  <  2:1666557/78‑1 (MQ=255)
                                                                       cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  >  1:1666557/1‑78 (MQ=255)
                                                                        aaaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACt                    >  1:1540240/1‑91 (MQ=255)
                                                                           aTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGAc                 <  1:3441133/91‑1 (MQ=255)
                                                                           aTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAACACAAAAAGAAAGTGCAACTGAc                 <  2:248005/91‑1 (MQ=255)
                                                                            ttAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACt                >  2:2627508/1‑91 (MQ=255)
                                                                              aaaTAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                       >  1:484685/1‑82 (MQ=255)
                                                                               tacaataTTCATCTGAATATTGTATTATTGAACAACGATGGTGGCAGTATTTTTTCATTTTTACCACAAAAAGAAAGTGCAACTGACTAtt             >  1:3295502/4‑91 (MQ=255)
                                                                               aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtt             <  1:2237739/91‑1 (MQ=255)
                                                                               aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtt             >  1:2496360/1‑91 (MQ=255)
                                                                               aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtg             >  2:3849700/1‑90 (MQ=255)
                                                                                ataataTTCGGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt            <  1:3714212/91‑1 (MQ=255)
                                                                                ataataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt            >  1:2306727/1‑91 (MQ=255)
                                                                                 taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTTTTTTTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  2:3259352/1‑91 (MQ=255)
                                                                                 taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  1:1946116/1‑91 (MQ=255)
                                                                                 caataTTCATCTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  1:2906081/2‑91 (MQ=255)
                                                                                    tatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGaa        >  2:3487946/1‑91 (MQ=255)
                                                                                      ttCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACg      <  2:2373677/91‑1 (MQ=255)
                                                                                        cAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGgt    >  2:370956/1‑91 (MQ=255)
                                                                                        cAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGgt    >  2:1418867/1‑91 (MQ=255)
                                                                                          gATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGgttg  <  2:2169643/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTG  >  USA300TCH1516_ALE/1049623‑1049803

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: