Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,107,402 T→C 12.7% intergenic (+37/‑131) pdhD → / → USA300HOU_RS05475 Dihydrolipoyl dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,107,4020TC12.7% 648.6 / 75.5 243intergenic (+37/‑131)pdhD/USA300HOU_RS05475Dihydrolipoyl dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (108/104);  new base C (5/26);  total (113/130)
Fisher's exact test for biased strand distribution p-value = 2.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 248 total aligned reads displayed.

ATGACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATAC  >  USA300TCH1516_ALE/1107312‑1107492
                                                                                          |                                                                                          
aTGACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCt                                                                                            <  2:2497088/91‑1 (MQ=255)
 tGACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTg                                                                                           <  2:3040390/91‑1 (MQ=255)
 tGACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTg                                                                                           <  1:2824412/91‑1 (MQ=255)
  gACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGt                                                                                          <  2:1113756/91‑1 (MQ=255)
   aCTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTa                                                                                         >  1:2402331/1‑91 (MQ=255)
    cTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAg                                                                                        <  1:1244176/91‑1 (MQ=255)
    cTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAg                                                                                        >  1:927793/1‑91 (MQ=255)
     tATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGc                                                                                       <  2:62347/91‑1 (MQ=255)
     tATGGAAGCAGCAGAAAAAGCTAGCGGATACCCAATCCATACACTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGc                                                                                       <  2:3624370/91‑1 (MQ=255)
        ggAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCata                                                                                    <  1:1794194/91‑1 (MQ=255)
          aaGCAGCAGAAAAAGCTATCGGATACCCAATCGATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATg                                                                                  <  1:3302077/91‑1 (MQ=255)
          aaGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATg                                                                                  <  1:1119342/91‑1 (MQ=255)
           agcagcAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGc                                                                                 <  1:3564041/91‑1 (MQ=255)
           agcagcAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGc                                                                                 <  2:3723579/91‑1 (MQ=255)
           agcagcAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGc                                                                                 <  1:2123734/91‑1 (MQ=255)
            gcagcaGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCt                                                                                <  2:1436366/91‑1 (MQ=255)
            gcagcaGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCt                                                                                <  1:2743886/91‑1 (MQ=255)
              agcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTAc                                                                              <  2:1922793/91‑1 (MQ=255)
              agcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTAc                                                                              <  1:3330171/91‑1 (MQ=255)
              agcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTAc                                                                              <  1:2753625/91‑1 (MQ=255)
              agcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTAc                                                                              <  2:2259644/91‑1 (MQ=255)
               gcagAAAGAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTAca                                                                             <  2:3448972/91‑2 (MQ=255)
               gcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTAca                                                                             <  1:541805/91‑2 (MQ=255)
                cagAAAAAGCTATCGGATACCCTATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                            <  2:2949941/91‑1 (MQ=255)
                cagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACgg                                                                            <  2:3179619/91‑1 (MQ=255)
                 agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                           <  1:1529764/91‑1 (MQ=255)
                 agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTATCATATGCTAcagc                                                                           <  2:1621971/91‑4 (MQ=255)
                  gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAATGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                          <  1:3266716/91‑1 (MQ=255)
                  gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                          <  2:2187790/91‑1 (MQ=255)
                  gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                          <  1:1331437/91‑1 (MQ=255)
                  gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                          <  1:585266/91‑1 (MQ=255)
                  gAAAAAGATATCGGATATGTAATGTATATAATGTAATAACTGATTATGTATAAAGATTAAGTAATTAAAAGATGTAGGATATGGTAAggat                                                                          >  2:358127/1‑89 (MQ=255)
                  gAAAAAGATATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTACAGCt                                                                          <  1:1937320/91‑1 (MQ=39)
                   aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtt                                                                         <  1:639534/91‑1 (MQ=255)
                   aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtt                                                                         <  1:2632374/91‑1 (MQ=255)
                   aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTACAGCtt                                                                         <  2:18394/91‑1 (MQ=255)
                    aaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCttt                                                                        >  2:1652580/1‑91 (MQ=255)
                     aaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttt                                                                       <  1:317562/91‑1 (MQ=255)
                       aGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttttt                                                                     <  2:2301468/91‑1 (MQ=255)
                         cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                   >  1:3461702/1‑91 (MQ=255)
                           aTCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGttt                                                                 >  2:544870/1‑91 (MQ=255)
                           aTCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGttt                                                                 >  1:3806301/1‑91 (MQ=255)
                            tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                <  2:1309458/91‑1 (MQ=255)
                            tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                >  1:1211153/1‑91 (MQ=255)
                            tCGGATACCCAAGCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                <  2:2649176/91‑1 (MQ=255)
                              ggATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAg                                                              <  2:3863701/91‑1 (MQ=255)
                                 tACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTa                                                           <  1:3673102/91‑1 (MQ=255)
                                  aCCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTaa                                                          <  2:1515460/91‑1 (MQ=255)
                                    ccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAg                                                        >  2:2922125/1‑91 (MQ=255)
                                    ccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAg                                                        >  1:2596845/1‑91 (MQ=255)
                                     cAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGt                                                       <  1:342016/91‑1 (MQ=255)
                                         ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                   >  1:689508/1‑91 (MQ=255)
                                         ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                   >  1:1067276/1‑91 (MQ=255)
                                            tACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTaa                                                >  2:2847618/1‑91 (MQ=255)
                                             aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                               >  1:2099276/1‑91 (MQ=255)
                                               aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                             <  1:484794/91‑1 (MQ=255)
                                               aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                             >  2:2342314/1‑91 (MQ=255)
                                               aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTATAGGTAAAGTAACGTAAGGa                                             >  1:2009114/1‑91 (MQ=255)
                                                aTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaa                                            <  2:2119683/91‑1 (MQ=255)
                                                 tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                           <  2:275639/91‑1 (MQ=255)
                                                 tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                           <  2:3160621/91‑1 (MQ=255)
                                                  gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAGGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAt                                          <  1:877624/91‑1 (MQ=255)
                                                  gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAt                                          <  1:1612828/91‑1 (MQ=255)
                                                   taataaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAtt                                         <  1:1189679/91‑1 (MQ=255)
                                                    aataaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTg                                        >  1:1833832/1‑91 (MQ=255)
                                                    aataaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTg                                        <  2:689508/91‑1 (MQ=255)
                                                       aaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTTCGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAtt                                     >  2:1021572/1‑91 (MQ=255)
                                                       aaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAtt                                     <  2:927600/91‑1 (MQ=255)
                                                         cTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTg                                   <  2:1501658/91‑1 (MQ=255)
                                                         cTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTg                                   <  2:1188219/91‑1 (MQ=255)
                                                          tGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTga                                  <  1:853129/91‑1 (MQ=255)
                                                           gATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTgag                                 >  1:2206485/1‑91 (MQ=255)
                                                             ttATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGat                               >  1:326019/1‑91 (MQ=255)
                                                               aTCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGatat                             <  1:2465623/91‑1 (MQ=255)
                                                                tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTTGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                            >  2:2871070/1‑91 (MQ=255)
                                                                tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                            <  1:3772984/91‑1 (MQ=255)
                                                                 cTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCg                           <  1:46247/91‑1 (MQ=255)
                                                                   ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                             >  1:1514722/1‑71 (MQ=255)
                                                                   ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGtt                         <  2:1771462/91‑1 (MQ=255)
                                                                   ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGtt                         >  2:2096177/1‑91 (MQ=255)
                                                                     aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                       <  2:3256150/91‑1 (MQ=255)
                                                                     aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                       >  1:802557/1‑91 (MQ=255)
                                                                     aaaGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTaa                       >  2:245896/1‑91 (MQ=255)
                                                                        gATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAt                    <  2:2817334/91‑1 (MQ=255)
                                                                         aTTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtg                   <  2:3103517/91‑1 (MQ=255)
                                                                          ttCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                           <  1:2147666/66‑1 (MQ=255)
                                                                          ttCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgt                  <  2:2758950/91‑1 (MQ=255)
                                                                           tCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgtg                 >  2:1084991/1‑91 (MQ=255)
                                                                            cAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAAGTAAGGAAATTGATTTGAGATATCGTTAACATGTGa                >  2:3506249/1‑91 (MQ=255)
                                                                              gTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAca              <  1:467018/91‑1 (MQ=255)
                                                                                  ttAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgca          >  1:3390115/1‑91 (MQ=255)
                                                                                   aaaaaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcat         >  2:1574259/2‑91 (MQ=255)
                                                                                    aaaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg        >  1:3630561/1‑91 (MQ=255)
                                                                                    aaaaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg        >  1:1416801/1‑91 (MQ=255)
                                                                                      aaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGtt      >  1:1700165/1‑91 (MQ=255)
                                                                                       aGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTta     <  1:222487/91‑1 (MQ=255)
                                                                                        gCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtat    <  1:2082275/91‑1 (MQ=255)
                                                                                         cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata   >  2:3211494/1‑91 (MQ=255)
                                                                                         cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata   >  1:2980759/1‑91 (MQ=255)
                                                                                          tGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGTCATGCATGTTATAc  <  2:3292389/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATGACTATGGAAGCAGCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATAC  >  USA300TCH1516_ALE/1107312‑1107492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: