Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,678,362 A→C 5.1% intergenic (+108/+42) 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,3620AC 5.1% 408.1 / 7.6 117intergenic (+108/+42)glyQS/recOGlycine‑‑tRNA ligase/DNA repair protein RecO
Reads supporting (aligned to +/- strand):  ref base A (63/48);  new base C (4/2);  total (67/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 130 total aligned reads displayed.

CATATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTAAAA  >  USA300TCH1516_ALE/1678272‑1678452
                                                                                          |                                                                                          
 atatGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAac                                                                                           >  1:63551/1‑91 (MQ=255)
 atatGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAac                                                                                           <  1:2159211/91‑1 (MQ=255)
  tatGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAaca                                                                                          <  1:709831/91‑1 (MQ=255)
  catatggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAaca                                                                                          >  1:1129851/5‑91 (MQ=255)
  catatggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAaca                                                                                          >  1:1130155/5‑91 (MQ=255)
    tgtgGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAAcaca                                                                                        <  1:1573163/91‑1 (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:121841/1‑91 (MQ=255)
       ggtggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCt                                                                                     >  1:2019931/1‑91 (MQ=255)
         tggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTAt                                                                                   >  1:2426192/1‑91 (MQ=255)
          ggCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTAtg                                                                                  <  1:1475286/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)
                 taATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTAc                                                                           >  1:1748018/1‑91 (MQ=255)
                   aTTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACtt                                                                         <  1:2949389/91‑1 (MQ=255)
                     ttATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGa                                                                       >  1:1926663/1‑91 (MQ=255)
                     ttATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGa                                                                       >  1:1661861/1‑91 (MQ=255)
                      tataAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGaa                                                                      >  1:3748240/1‑91 (MQ=255)
                        taAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATa                                                                    >  1:445013/1‑91 (MQ=255)
                         aaaaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCGTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATAc                                                                   >  1:2248517/1‑91 (MQ=255)
                          aaaaGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATAcc                                                                  >  1:1558144/1‑91 (MQ=255)
                            aaGGGATGCTGGTTAGTAAAATTGTCAGTATCACTGTTCTATTCAAAAAGAGATTCAAATAAACTCAACGATGTGGCTACTTTAATCCCtt                                                                <  1:1619431/91‑1 (MQ=255)
                              gagaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCtttt                                                              >  1:727297/1‑91 (MQ=255)
                               agaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTtta                                                             <  1:3624556/91‑1 (MQ=255)
                               agaTACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTtta                                                             <  1:1359732/91‑1 (MQ=255)
                                   aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                         >  1:3314180/1‑91 (MQ=255)
                                   aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                         >  1:2272172/1‑91 (MQ=255)
                                   aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                         >  1:2271993/1‑91 (MQ=255)
                                   aCTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTtattat                                                         >  1:1086822/1‑91 (MQ=255)
                                       gTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTa                                                     <  1:2962895/91‑1 (MQ=255)
                                        ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAATGTATTCAAATAAACACGTCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                                    >  1:1625588/1‑91 (MQ=255)
                                        ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAGTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                                    >  1:3286144/1‑91 (MQ=255)
                                        ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                                    >  1:2173069/1‑91 (MQ=255)
                                        ttAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTAt                                                    <  1:2496029/91‑1 (MQ=255)
                                         tAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATg                                                   <  1:1443775/91‑1 (MQ=255)
                                         tAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATg                                                   <  1:1765000/91‑1 (MQ=255)
                                          aGTAAAATTGCTAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGt                                                  <  1:3760988/91‑1 (MQ=255)
                                          aGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGt                                                  <  1:1089695/91‑1 (MQ=255)
                                          aGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACAAATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGt                                                  >  1:680679/1‑91 (MQ=255)
                                            tAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTa                                                >  1:2734615/1‑91 (MQ=255)
                                                aTTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGt                                            >  1:2523045/1‑91 (MQ=255)
                                                aTTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGt                                            <  1:3838895/91‑1 (MQ=255)
                                                 ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                           <  1:4079713/91‑1 (MQ=255)
                                                 ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                           <  1:2399938/91‑1 (MQ=255)
                                                 ttGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGtt                                           <  1:4079685/91‑1 (MQ=255)
                                                    ccAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCa                                        >  1:3949755/1‑91 (MQ=255)
                                                    ccAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCa                                        >  1:1769023/1‑91 (MQ=255)
                                                     cAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCaa                                       >  1:462393/1‑91 (MQ=255)
                                                       gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                     >  1:3298335/1‑91 (MQ=255)
                                                       gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATc                                     >  1:2422640/1‑91 (MQ=255)
                                                       gTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATCGTTCCAATc                                     >  1:3963669/1‑91 (MQ=255)
                                                        tATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCt                                    <  1:2264363/91‑1 (MQ=255)
                                                           ctctTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCtttt                                 <  1:2481051/91‑1 (MQ=255)
                                                              ttttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTaa                              <  1:3224095/91‑1 (MQ=255)
                                                              ttttCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTaa                              >  1:2178688/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)
                                                                    aTTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGtt                        <  1:3121538/91‑1 (MQ=255)
                                                                     ttCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTg                       >  1:1431761/1‑91 (MQ=255)
                                                                     ttCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTg                       <  1:1301201/91‑1 (MQ=255)
                                                                     ttCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTg                       >  1:734866/1‑91 (MQ=255)
                                                                      tCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGa                      <  1:1298680/91‑1 (MQ=255)
                                                                       cAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGAt                     >  1:808133/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)
                                                                             ggTATTCAAGTAACCACATAGATGTGTTTATTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGttt               >  1:3181743/1‑91 (MQ=18)
                                                                             ggTATTCAAGTAACCACATAGATGTGTTTATTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGttt               >  1:3180809/1‑91 (MQ=18)
                                                                             ggTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGttt               >  1:2744903/1‑91 (MQ=255)
                                                                              gTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTc              >  1:1427249/1‑91 (MQ=255)
                                                                                aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg            <  1:3480070/91‑1 (MQ=255)
                                                                                aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg            >  1:3812932/1‑91 (MQ=255)
                                                                                aTTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTg            <  1:2854177/91‑1 (MQ=255)
                                                                                 ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTgg           >  1:1044127/1‑90 (MQ=255)
                                                                                 ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTgg           >  1:1045145/1‑90 (MQ=255)
                                                                                 ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa           >  1:204585/1‑91 (MQ=255)
                                                                                 ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa           >  1:508507/1‑91 (MQ=255)
                                                                                 ttCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGa           <  1:1822266/91‑1 (MQ=255)
                                                                                   cAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACt         >  1:224941/1‑91 (MQ=255)
                                                                                   cAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACt         <  1:1826829/91‑1 (MQ=255)
                                                                                   cAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACt         <  1:1826699/91‑1 (MQ=255)
                                                                                     aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt       >  1:2827059/1‑91 (MQ=255)
                                                                                     aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt       >  1:848933/1‑91 (MQ=255)
                                                                                     aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt       >  1:848971/1‑91 (MQ=255)
                                                                                     aaTAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACttt       >  1:3734257/1‑91 (MQ=255)
                                                                                          acacATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTaaaa  <  1:460482/91‑1 (MQ=255)
                                                                                          acacATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTaaaa  <  1:460535/91‑1 (MQ=255)
                                                                                          acacATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTaaaa  <  1:2365219/91‑1 (MQ=255)
                                                                                          acacATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTaaaa  <  1:2365546/91‑1 (MQ=255)
                                                                                          acacATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTaaaa  >  1:403304/1‑91 (MQ=255)
                                                                                          |                                                                                          
CATATGTGGTGGCTATATAATTTATAAAAAGAGATACTGGTTAGTAAAATTGCCAGTATCTCTTTTCTATTCAAAAAGGTATTCAAATAAACACATCTATGTGGCTACTTGAATACCTTTTATTATTTTATGTTATTGTTCCAATCTTTTTAATTGGTTGATTAGTTTCTGACTTTTAAAA  >  USA300TCH1516_ALE/1678272‑1678452

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: