Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,530,791 C→T 6.4% intergenic (‑42/+231) yefM ← / ← bdbD Antitoxin YefM/Disulfide bond formation protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,530,7910CT 6.4% 632.1 / 19.3 194intergenic (‑42/+231)yefM/bdbDAntitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand):  ref base C (79/102);  new base T (8/5);  total (87/107)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTAC  >  USA300TCH1516_ALE/2530701‑2530881
                                                                                          |                                                                                          
tAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATc                                                                                            >  2:3721404/1‑91 (MQ=255)
   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                         <  2:1543412/91‑1 (MQ=255)
   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                         <  2:1659663/91‑1 (MQ=255)
    tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTa                                                                                        <  2:1528865/91‑1 (MQ=255)
    tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                        <  1:726/91‑1 (MQ=255)
     gCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTAc                                                                                       <  1:852490/91‑1 (MQ=255)
      ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTACCACACCTCAATTGTACAACTATCTGTACa                                                                                      >  1:3106891/1‑91 (MQ=255)
       cTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACaa                                                                                     <  1:2156027/91‑1 (MQ=255)
       cTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACaa                                                                                     >  1:2103721/1‑91 (MQ=255)
        tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                                    >  2:962999/1‑91 (MQ=255)
        tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                                    >  1:956725/1‑91 (MQ=255)
        tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                                    <  2:196827/91‑1 (MQ=255)
         ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                                   >  2:1069729/1‑91 (MQ=255)
         ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCGTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                                   <  1:3271889/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  <  2:2225864/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  >  1:1368234/1‑91 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  <  2:3554979/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  <  2:2225630/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                  <  1:2703386/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACATCAATTGTACAAATATCTGTACAATaa                                                                                  <  2:3221258/91‑1 (MQ=255)
          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATATATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTACAAttg                                                                                  >  2:2942303/1‑89 (MQ=39)
             aTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAAc                                                                               <  1:1791287/91‑1 (MQ=255)
               tCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACaa                                                                             >  1:3708111/1‑91 (MQ=255)
                 tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                           >  2:296218/1‑91 (MQ=255)
                 tGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGc                                                                           >  1:713300/1‑91 (MQ=255)
                  gTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCg                                                                          >  1:3491514/1‑91 (MQ=255)
                   tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                         >  2:2403636/1‑91 (MQ=255)
                      aGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGacac                                                                      <  2:901209/91‑1 (MQ=255)
                       gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                     <  2:302530/91‑1 (MQ=255)
                       gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                     <  2:3879381/91‑1 (MQ=255)
                         gTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGa                                                                   <  1:2809208/91‑1 (MQ=255)
                          tATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaa                                                                  <  2:1058666/91‑1 (MQ=255)
                              aataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaa                                                              >  1:1672953/1‑91 (MQ=255)
                                 aatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAact                                                           >  2:1386918/1‑91 (MQ=255)
                                    ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACTCAGa                                                                   <  1:3056996/80‑1 (MQ=255)
                                    ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactact                                                        <  2:1216710/91‑1 (MQ=255)
                                    ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactact                                                        <  2:1191022/91‑1 (MQ=255)
                                     tttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAActacta                                                       >  2:2646371/1‑91 (MQ=255)
                                      ttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                      <  1:2761508/91‑1 (MQ=255)
                                      ttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                      <  1:1519390/91‑1 (MQ=255)
                                       taatatTCATATTAATCGCCTCCTTTTAACACACCTCTATTGTACAAATATCTGTACCATAAAACAAGCGACACAGAATAAAACTCCTaaa                                                     >  2:573453/1‑90 (MQ=255)
                                         ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                                   <  1:2796270/91‑1 (MQ=255)
                                         ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                                   >  2:2083701/1‑91 (MQ=255)
                                         ataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCa                                                   <  1:1393854/91‑1 (MQ=255)
                                           aatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAg                                                 >  1:1156067/1‑91 (MQ=255)
                                           aatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAg                                                 >  1:1122345/1‑91 (MQ=255)
                                            atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGa                                                <  2:2781985/91‑1 (MQ=255)
                                             tCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAt                                               <  1:3540987/91‑1 (MQ=255)
                                              cATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAtt                                              >  2:616991/1‑91 (MQ=255)
                                                 atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTTTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                           <  1:3870785/91‑1 (MQ=255)
                                                 atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                           >  2:288773/1‑91 (MQ=255)
                                                 atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                           <  2:246304/91‑1 (MQ=255)
                                                 atTTATCACCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                           <  1:1949385/91‑1 (MQ=255)
                                                   ttATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCa                                         >  2:1146994/1‑91 (MQ=255)
                                                    tATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCAc                                        >  1:2558207/1‑91 (MQ=255)
                                                     aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                       <  1:1130919/91‑1 (MQ=255)
                                                     aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                       <  1:1516670/91‑1 (MQ=255)
                                                     aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                       <  1:3209782/91‑1 (MQ=255)
                                                       cGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACggg                                     >  2:1795248/1‑91 (MQ=255)
                                                        gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                                    <  1:967120/91‑1 (MQ=255)
                                                        gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                                    >  2:449831/1‑91 (MQ=255)
                                                        gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                                    >  2:3296639/1‑91 (MQ=255)
                                                         cctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGga                                   >  2:1391560/1‑91 (MQ=255)
                                                          ctcctTTTAACACACTTCAATTGTAGTAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaa                                  <  2:3367308/91‑1 (MQ=255)
                                                           tcctTTTAACACACCTCAATTGTACAAATATCTGTATAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                                 <  2:3052973/91‑1 (MQ=255)
                                                           tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                                 <  1:1901356/91‑1 (MQ=255)
                                                             ctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaagaa                               >  1:2219382/1‑91 (MQ=255)
                                                               tttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATc                             >  2:803988/1‑91 (MQ=255)
                                                               tttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATc                             <  1:2679591/91‑1 (MQ=255)
                                                                ttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                            <  2:2404599/91‑1 (MQ=255)
                                                                ttAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                            <  1:2420152/91‑1 (MQ=255)
                                                                ttAACACACCTCAATTGCACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCa                            <  1:2313680/91‑1 (MQ=255)
                                                                  aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                          <  1:736530/91‑1 (MQ=255)
                                                                  aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                          >  1:1229315/1‑91 (MQ=255)
                                                                  aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                          <  2:1709040/91‑1 (MQ=255)
                                                                  aaCACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAAt                          <  1:3590803/91‑1 (MQ=255)
                                                                    cacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAc                        >  1:1776220/1‑91 (MQ=255)
                                                                     acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                       <  1:3739619/91‑1 (MQ=255)
                                                                      cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                      <  1:2862236/91‑1 (MQ=255)
                                                                      cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                      >  2:271254/1‑91 (MQ=255)
                                                                      cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                      >  1:604139/1‑91 (MQ=255)
                                                                       acCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                     >  2:1167949/1‑91 (MQ=255)
                                                                         cTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATc                   <  1:3238945/91‑1 (MQ=255)
                                                                          tCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAGTTCGCACGGGGAAGAATCAATACCTATcc                  <  2:3385053/91‑1 (MQ=255)
                                                                            aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg                <  2:2129234/91‑1 (MQ=255)
                                                                            aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg                <  2:2130124/91‑1 (MQ=255)
                                                                             aTTGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGt               >  2:596711/1‑91 (MQ=255)
                                                                             aTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGt               >  2:1330825/1‑91 (MQ=255)
                                                                              ttGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa              >  1:1700512/1‑91 (MQ=255)
                                                                              ttGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa              <  1:2179232/91‑1 (MQ=255)
                                                                              ttGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa              <  2:1863841/91‑1 (MQ=255)
                                                                               tGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc             >  2:707228/1‑91 (MQ=255)
                                                                               tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc             <  1:256613/91‑1 (MQ=255)
                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt            <  2:3332892/91‑1 (MQ=255)
                                                                                  aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa          >  2:3498574/1‑91 (MQ=255)
                                                                                   cAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGAt         >  2:3519413/1‑91 (MQ=255)
                                                                                     aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAATTACTAATCCAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc       >  1:3092560/1‑91 (MQ=255)
                                                                                     aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc       <  2:3460545/91‑1 (MQ=255)
                                                                                         tCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCcta   >  1:2958148/1‑91 (MQ=255)
                                                                                         tCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCcta   >  2:115088/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTAC  >  USA300TCH1516_ALE/2530701‑2530881

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: ‑