Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 659,372 T→G 9.4% intergenic (+178/‑734) USA300TCH1516_00590 → / → USA300HOU_RS03165 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE659,3720TG 9.4% 797.4 / 53.1 249intergenic (+178/‑734)USA300TCH1516_00590/USA300HOU_RS03165hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (110/115);  new base G (11/13);  total (121/128)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

TACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGAT  >  USA300TCH1516_ALE/659282‑659462
                                                                                          |                                                                                          
tACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATtttttattttt                                                                                            >  1:1924819/1‑91 (MQ=255)
  cGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTttttatttttat                                                                                          >  2:1527772/1‑91 (MQ=255)
  cGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTttttatttttat                                                                                          >  1:413901/1‑91 (MQ=255)
  cGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTttttatttttat                                                                                          >  1:1386231/1‑91 (MQ=255)
   ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                         >  2:778772/1‑91 (MQ=255)
   ggTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATc                                                                                         >  1:571647/1‑91 (MQ=255)
    gTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCt                                                                                        >  2:3668263/1‑91 (MQ=255)
     tAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCta                                                                                       >  1:3456690/1‑91 (MQ=255)
      aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                      >  2:645519/1‑91 (MQ=255)
      aaCTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaa                                                                                      <  2:343900/91‑1 (MQ=255)
         ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                   <  2:1383617/91‑1 (MQ=255)
         ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                   <  1:3695096/91‑1 (MQ=255)
         ttcttTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                   <  1:2117516/91‑1 (MQ=255)
         ttctcTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCtaaata                                                                                   <  1:3228460/91‑1 (MQ=255)
             ttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaa                                                                               <  2:1973765/91‑1 (MQ=255)
             ttAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATATAAATaaaaa                                                                               <  2:3528372/91‑1 (MQ=255)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATaaaaaa                                                                              >  1:1589875/1‑91 (MQ=255)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                              <  1:1852231/91‑1 (MQ=37)
              tAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTAGATAAAAATaaaaaa                                                                              <  1:291226/91‑1 (MQ=37)
                aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                            >  2:2912706/1‑91 (MQ=255)
                aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                            >  2:2913318/1‑91 (MQ=255)
                aTGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGACGTTAATTTTTTATTTTTATCTAAATAAAAAAtt                                                                            >  1:3176790/1‑91 (MQ=255)
                  gTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTaa                                                                          >  2:1782966/1‑91 (MQ=255)
                     cAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                       >  2:268208/1‑91 (MQ=255)
                     aaGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAt                                                                       <  2:3624116/90‑1 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                      <  2:2861875/91‑1 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                      >  1:677777/1‑91 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                      >  2:1923058/1‑91 (MQ=255)
                      aGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACAtt                                                                      <  2:1904950/91‑1 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                   <  2:319415/91‑1 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                   <  1:2441082/91‑1 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                   <  2:2087779/91‑1 (MQ=255)
                         ttGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaa                                                                   <  1:812570/91‑1 (MQ=255)
                          tGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTaaaa                                                                  <  1:1976026/91‑1 (MQ=255)
                            aaaaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAt                                                                <  2:1219346/91‑1 (MQ=255)
                              aaaCAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAttt                                                              >  1:1769719/1‑91 (MQ=255)
                                acaacTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAAttttt                                                            <  2:1878224/91‑1 (MQ=255)
                                   acTAGCAATCACTGTATTTTTTATCATGTTAATNGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAAt                                                         <  1:1177323/91‑1 (MQ=255)
                                    cTAGCAATCACTGTATTTTTTATCATGTTAATAGTGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                        <  2:1426097/91‑1 (MQ=255)
                                    cTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                        >  2:1707867/1‑91 (MQ=255)
                                    cTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATa                                                        >  1:2262379/1‑91 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                       <  2:376125/91‑1 (MQ=255)
                                     tAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                       >  1:2072525/1‑91 (MQ=255)
                                     tAGCAATCACTGTATTTTTTAGCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAc                                                       <  2:2335094/91‑1 (MQ=255)
                                      aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                      >  1:3125853/1‑91 (MQ=255)
                                      aGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATAcc                                                      <  2:2448372/91‑1 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                     >  1:326437/1‑91 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                     >  2:1464785/1‑91 (MQ=255)
                                       gCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTGTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCa                                                     >  2:2988685/1‑91 (MQ=255)
                                        cAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCaa                                                    <  1:123245/91‑1 (MQ=255)
                                         aaTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAAc                                                   >  1:868368/1‑91 (MQ=255)
                                         aaTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAAc                                                   <  2:2566773/91‑1 (MQ=255)
                                          aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                  >  2:2200433/1‑91 (MQ=255)
                                          aTCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACg                                                  <  1:3838026/91‑1 (MQ=255)
                                           tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                 <  1:3567912/91‑1 (MQ=255)
                                           tCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTGATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGt                                                 >  1:1155559/1‑91 (MQ=255)
                                            cACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTg                                                <  2:1122597/91‑1 (MQ=255)
                                            cACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTg                                                <  2:1738369/91‑1 (MQ=255)
                                             aCTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTgg                                               <  1:245298/91‑1 (MQ=255)
                                              cTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGt                                              <  1:2329725/91‑1 (MQ=255)
                                                     ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                       >  2:734584/1‑91 (MQ=255)
                                                     ttttATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAg                                       <  1:488657/91‑1 (MQ=255)
                                                          tCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATag                                  >  2:2512393/1‑91 (MQ=255)
                                                          tCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATag                                  >  2:1608178/1‑91 (MQ=255)
                                                          tCATGTTAATAGGGATGTTAATTTTTTATTTTTAGCTAGATAAAAAATGAACATTAAAAATTTTTAATGGCAGCGTGGTTAATAAGAATag                                  <  2:992157/91‑1 (MQ=255)
                                                           cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACTAACGTGGTTAATAAGATTaga                                 >  2:530927/1‑91 (MQ=255)
                                                           cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                 >  1:3464095/1‑91 (MQ=255)
                                                           cATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATaga                                 >  2:497049/1‑91 (MQ=255)
                                                            aTGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATagag                                <  1:2661890/91‑1 (MQ=255)
                                                               ttAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTg                             <  1:2880738/91‑1 (MQ=255)
                                                                 aaTAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTa                           >  1:914063/1‑91 (MQ=255)
                                                                   tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                         <  2:2746909/91‑1 (MQ=255)
                                                                   tAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTAtt                         <  1:1885543/91‑1 (MQ=255)
                                                                    tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTa                        >  1:1683861/5‑91 (MQ=21)
                                                                    tttaaTGTTAATTTTTTATTTAGATAAAAATAAAAAATTAACATTAAAAATTTTCAATACCAACGAGGTTAATAAGAATAGAGTTGTATTa                        >  2:2591309/5‑91 (MQ=18)
                                                                     ttaaTGTTAATTTTATATTTCGATAAAAATAAAAAATTAACATTAATAATTTTTAATACCCCCGTGGTTAATAAGACTAGAGTTGTATTaa                       >  2:469308/4‑91 (MQ=18)
                                                                     gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                       >  2:1575044/1‑91 (MQ=255)
                                                                     gggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTaa                       >  2:324440/1‑91 (MQ=255)
                                                                      ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                      >  1:1269112/1‑91 (MQ=255)
                                                                      ggATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAAt                      >  1:1888457/1‑91 (MQ=255)
                                                                       gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTACTAAGAATAGAGTTGTATTAATa                     >  2:3132301/1‑91 (MQ=255)
                                                                       gATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATa                     >  1:1011792/1‑91 (MQ=255)
                                                                        aTGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAc                    >  2:1269145/1‑91 (MQ=255)
                                                                         tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                   >  2:2488017/1‑91 (MQ=255)
                                                                         tGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATAcc                   <  1:3575481/91‑1 (MQ=255)
                                                                          gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                  >  2:3058715/1‑91 (MQ=255)
                                                                          gTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCg                  >  2:3352905/1‑91 (MQ=255)
                                                                            tAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGAc                >  2:3100658/1‑91 (MQ=255)
                                                                             aaTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACt               <  1:1012548/91‑1 (MQ=255)
                                                                              aTTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACtt              <  1:108714/91‑1 (MQ=255)
                                                                               ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTATTACCGACttt             >  2:412996/1‑91 (MQ=255)
                                                                               ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt             >  1:1817414/1‑91 (MQ=255)
                                                                               ttttttATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACttt             >  1:2844443/1‑91 (MQ=255)
                                                                                   ttatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCaaa         >  1:2789079/1‑91 (MQ=255)
                                                                                    tatttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAAc        <  1:3371396/91‑1 (MQ=255)
                                                                                      tttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACaa      >  2:103980/1‑91 (MQ=255)
                                                                                       ttttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAg     <  2:1989525/91‑1 (MQ=255)
                                                                                        tttatttaGATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg    >  2:197843/1‑91 (MQ=255)
                                                                                        tttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAgg    >  1:3749695/1‑91 (MQ=255)
                                                                                         ttatCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGa   <  2:1172064/91‑1 (MQ=255)
                                                                                          |                                                                                          
TACGGTAACTTCTTTAATGTTCAGTTTGAAAAACAACTAGCAATCACTGTATTTTTTATCATGTTAATAGGGATGTTAATTTTTTATTTTTATCTAAATAAAAAATTAACATTAAAAATTTTTAATACCAACGTGGTTAATAAGAATAGAGTTGTATTAATACCGACTTTCAAACAAGGAT  >  USA300TCH1516_ALE/659282‑659462

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: