Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,678,358 A→G 5.5% intergenic (+104/+46) glyQS → / ← recO Glycine‑‑tRNA ligase/DNA repair protein RecO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,678,3580AG 5.5% 416.3 / 8.5 123intergenic (+104/+46)glyQS/recOGlycine‑‑tRNA ligase/DNA repair protein RecO
Reads supporting (aligned to +/- strand):  ref base A (66/50);  new base G (5/2);  total (71/52)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

CCACCATATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTT  >  USA300TCH1516_ALE/1678268‑1678447
                                                                                          |                                                                                         
ccaccaCATATGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCaaa                                                                                           >  1:300273/1‑91 (MQ=255)
ccaccaCATATGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCaaa                                                                                           >  1:300141/1‑91 (MQ=255)
ccaccaCATATGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCaaa                                                                                           >  1:478236/1‑91 (MQ=255)
 caccaTATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAAt                                                                                          >  1:3960048/1‑91 (MQ=255)
  accaTATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATa                                                                                         >  1:107840/1‑91 (MQ=255)
   ccaTATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATaa                                                                                        >  1:864164/1‑91 (MQ=255)
    caTATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATaaa                                                                                       >  1:461072/1‑91 (MQ=255)
      tatGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAaca                                                                                     <  1:709831/91‑1 (MQ=255)
      catatggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAaca                                                                                     >  1:1129851/5‑91 (MQ=255)
        tatggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAAcaca                                                                                   >  1:3006875/3‑91 (MQ=255)
         atggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACAt                                                                                  >  1:510000/2‑91 (MQ=255)
          tggtggCTATATAATTGATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATc                                                                                 <  1:2822710/91‑1 (MQ=255)
           ggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCt                                                                                >  1:2019931/1‑91 (MQ=255)
           ggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCt                                                                                >  1:121841/1‑91 (MQ=255)
           ggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCt                                                                                >  1:1132287/1‑91 (MQ=255)
              ggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTAtg                                                                             <  1:1475286/91‑1 (MQ=255)
              ggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTAtg                                                                             <  1:809890/91‑1 (MQ=255)
                cTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTAtgtg                                                                           <  1:1840222/91‑1 (MQ=255)
                  atataATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGc                                                                         <  1:4106852/91‑1 (MQ=255)
                  atataATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGc                                                                         >  1:2175476/1‑91 (MQ=255)
                   tataATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCt                                                                        >  1:2244220/1‑91 (MQ=255)
                    ataATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTa                                                                       <  1:2941516/91‑1 (MQ=255)
                     taATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTAc                                                                      >  1:1748018/1‑91 (MQ=255)
                       aTTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACtt                                                                    <  1:2949389/91‑1 (MQ=255)
                          tataAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGaa                                                                 >  1:3748240/1‑91 (MQ=255)
                             aaaaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCGTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATAc                                                              >  1:2248517/1‑91 (MQ=255)
                              aaaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATAcc                                                             >  1:1558144/1‑91 (MQ=255)
                               aaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAGTAGCCACATAGATGTGTTTATTTGAATACCt                                                            <  1:1691559/91‑1 (MQ=17)
                               aaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAGTAGCCACATAGATGTGTTTATTTGAATACCt                                                            <  1:1933101/91‑1 (MQ=17)
                                aaGGGATGCTGGTTAGTAAAATTGTCAGTATCACTGTTCTATTCAAAAAGAGATTCAAATAAACTCAACGATGTGGCTACTTTAATCCCtt                                                           <  1:1619431/91‑1 (MQ=255)
                                  gagaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCtttt                                                         >  1:727297/1‑91 (MQ=255)
                                   agaTGCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTtta                                                        <  1:1359659/91‑1 (MQ=255)
                                   agaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTtta                                                        <  1:1359732/91‑1 (MQ=255)
                                   agaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTtta                                                        <  1:3624556/91‑1 (MQ=255)
                                    gaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTttat                                                       >  1:710590/1‑91 (MQ=255)
                                     aTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTttatt                                                      >  1:556354/1‑91 (MQ=255)
                                      tACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTttatta                                                     <  1:1616039/91‑1 (MQ=255)
                                       aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                    >  1:2171182/1‑91 (MQ=255)
                                       aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                    >  1:3314180/1‑91 (MQ=255)
                                       aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                    >  1:1086822/1‑91 (MQ=255)
                                       aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                    >  1:2271993/1‑91 (MQ=255)
                                       aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                    >  1:2272172/1‑91 (MQ=255)
                                            ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAATGTATTCAAATAAACACGTCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                               >  1:1625588/1‑91 (MQ=255)
                                            ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAGTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                               >  1:3286144/1‑91 (MQ=255)
                                            ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                               >  1:2173069/1‑91 (MQ=255)
                                             tAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATg                                              <  1:1765000/91‑1 (MQ=255)
                                              aGTAAAATTGCTAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGt                                             <  1:3760988/91‑1 (MQ=255)
                                              aGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACAAATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGt                                             >  1:680679/1‑91 (MQ=255)
                                                tAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTa                                           >  1:2734615/1‑91 (MQ=255)
                                                    aTTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGt                                       >  1:1458751/1‑91 (MQ=255)
                                                    aTTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGt                                       <  1:3838895/91‑1 (MQ=255)
                                                    aTTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGt                                       >  1:2523045/1‑91 (MQ=255)
                                                     ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                      <  1:4079685/91‑1 (MQ=255)
                                                     ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                      <  1:4079713/91‑1 (MQ=255)
                                                     ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                      <  1:2399938/91‑1 (MQ=255)
                                                        ccAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCa                                   >  1:3949755/1‑91 (MQ=255)
                                                         cAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCaa                                  >  1:462393/1‑91 (MQ=255)
                                                           gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                <  1:2659095/91‑1 (MQ=255)
                                                           gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                >  1:3298335/1‑91 (MQ=255)
                                                           gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                <  1:1099656/91‑1 (MQ=255)
                                                           gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                >  1:2422640/1‑91 (MQ=255)
                                                            tATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCt                               <  1:2264363/91‑1 (MQ=255)
                                                               ctctTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCtttt                            <  1:2481051/91‑1 (MQ=255)
                                                                tctTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCttttt                           >  1:947933/1‑91 (MQ=255)
                                                                  ttttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTaa                         <  1:3224095/91‑1 (MQ=255)
                                                                  ttttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTaa                         >  1:2472647/1‑91 (MQ=255)
                                                                   tttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAAt                        <  1:2754042/91‑1 (MQ=255)
                                                                   tttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAAt                        >  1:2658692/1‑91 (MQ=255)
                                                                   tttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAAt                        <  1:2297741/91‑1 (MQ=255)
                                                                    ttCTATTCAAAAGGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAAtt                       <  1:2458216/91‑1 (MQ=255)
                                                                    ttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAAtt                       <  1:57572/91‑1 (MQ=255)
                                                                      cTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTgg                     <  1:2803259/91‑1 (MQ=255)
                                                                        aTTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGtt                   <  1:3121538/91‑1 (MQ=255)
                                                                         ttCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTg                  >  1:1431761/1‑91 (MQ=255)
                                                                          tCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGa                 <  1:1298680/91‑1 (MQ=255)
                                                                           cAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGAt                >  1:808133/1‑91 (MQ=255)
                                                                           cAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGAt                >  1:4010772/1‑91 (MQ=255)
                                                                              aaaGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAg             <  1:1946275/91‑1 (MQ=255)
                                                                               aaGGTATTCAAGTAGCCACATAGATGTGTTTATTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGt            >  1:930192/1‑91 (MQ=17)
                                                                               aaGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGt            >  1:1101666/1‑91 (MQ=255)
                                                                                aGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGtt           <  1:3929997/91‑1 (MQ=255)
                                                                                aGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGtt           <  1:1761560/91‑1 (MQ=255)
                                                                                aGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGtt           <  1:1839943/91‑1 (MQ=255)
                                                                                aGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGtt           <  1:1839922/91‑1 (MQ=255)
                                                                                 ggTATTCAAGTAACCACATAGATGTGTTTATTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGttt          >  1:3180809/1‑91 (MQ=18)
                                                                                 ggTATTCAAGTAACCACATAGATGTGTTTATTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGttt          >  1:3181743/1‑91 (MQ=18)
                                                                                  gTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTc         >  1:1427249/1‑91 (MQ=255)
                                                                                  gTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTc         >  1:2415020/1‑91 (MQ=255)
                                                                                    aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg       >  1:3812932/1‑91 (MQ=255)
                                                                                    aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg       <  1:3480070/91‑1 (MQ=255)
                                                                                    aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg       <  1:2854177/91‑1 (MQ=255)
                                                                                     ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTgg      >  1:1045145/1‑90 (MQ=255)
                                                                                     ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTgg      >  1:1044127/1‑90 (MQ=255)
                                                                                     ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa      <  1:1822224/91‑1 (MQ=255)
                                                                                     ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa      <  1:1822266/91‑1 (MQ=255)
                                                                                     ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa      >  1:204585/1‑91 (MQ=255)
                                                                                       cAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACt    <  1:1826699/91‑1 (MQ=255)
                                                                                       cAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACt    >  1:224941/1‑91 (MQ=255)
                                                                                         aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt  >  1:848933/1‑91 (MQ=255)
                                                                                         aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt  >  1:3734257/1‑91 (MQ=255)
                                                                                          |                                                                                         
CCACCATATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTT  >  USA300TCH1516_ALE/1678268‑1678447

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: