Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,825,587 C→T 6.1% intergenic (‑402/+446) cap8A_2 ← / ← icaR Capsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,825,5870CT 6.1% 371.1 / 10.2 114intergenic (‑402/+446)cap8A_2/icaRCapsular polysaccharide type 8 biosynthesis protein cap8A/Biofilm operon icaADBC HTH‑type negative transcriptional regulator IcaR
Reads supporting (aligned to +/- strand):  ref base C (65/42);  new base T (2/5);  total (67/47)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 126 total aligned reads displayed.

TTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATA  >  USA300TCH1516_ALE/2825498‑2825677
                                                                                         |                                                                                           
tttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACAcc                                                                                            <  1:3429525/91‑1 (MQ=255)
tttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACAcc                                                                                            <  1:1621816/91‑1 (MQ=255)
tttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACAcc                                                                                            <  1:3941375/91‑1 (MQ=255)
 ttACAACACCGAGTGTGCCCGTACTCGGTAATTCAATACCTTGGGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCt                                                                                           <  1:4149838/91‑1 (MQ=255)
 ttACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCt                                                                                           >  1:322583/1‑91 (MQ=255)
  tACAACACCTAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCGACACCtt                                                                                          >  1:2085397/1‑91 (MQ=255)
   acaacaCCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCttt                                                                                         <  1:2767655/91‑1 (MQ=255)
    caacaCCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCtttt                                                                                        >  1:2958027/1‑91 (MQ=255)
    caacaCCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCtttt                                                                                        <  1:1676719/91‑1 (MQ=255)
     aacaCCGAGCGTGCCCGTACTCAGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTa                                                                                       <  1:4105360/91‑1 (MQ=255)
      acacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaa                                                                                      <  1:73492/91‑1 (MQ=255)
      acacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaa                                                                                      >  1:3384815/1‑91 (MQ=255)
      acacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaa                                                                                      >  1:3634958/1‑91 (MQ=255)
       cacCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCtttttaa                                                                                     <  1:676983/91‑4 (MQ=255)
        acCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaaa                                                                                    >  1:3259173/1‑91 (MQ=255)
         ccGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTaaaaa                                                                                   >  1:466748/1‑91 (MQ=255)
         ccGAGCGGGCCCGTACGCGGTAATTCAAAACCTTGCGTAACCCGTCACTGCGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTtaaaa                                                                                   <  1:3442963/91‑1 (MQ=255)
          cGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTTAAAAg                                                                                  <  1:2941950/91‑1 (MQ=255)
          cGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTTAAAAg                                                                                  <  1:1507949/91‑1 (MQ=255)
          cGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAg                                                                                  >  1:2771667/1‑91 (MQ=255)
            aGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACATCTTTTTAAAAggt                                                                                <  1:672724/91‑3 (MQ=255)
             gCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtg                                                                               >  1:3089429/1‑91 (MQ=255)
              cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                              >  1:3756630/1‑91 (MQ=255)
              cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                              >  1:3817557/1‑91 (MQ=255)
              cGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgt                                                                              >  1:2352998/1‑91 (MQ=255)
               gTGCCCGTAGTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgtg                                                                             <  1:3976729/91‑1 (MQ=255)
               gTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgtg                                                                             <  1:397334/91‑1 (MQ=255)
               gTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGGTAATGATAATAAAGCCCACACCTTTTAAAAAGAtgtg                                                                             <  1:3502676/91‑1 (MQ=255)
                   ccGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTa                                                                         >  1:1298982/1‑91 (MQ=255)
                     gTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAAt                                                                       >  1:1844019/1‑91 (MQ=255)
                       aCTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAAttt                                                                     <  1:1775390/91‑1 (MQ=255)
                       aCTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAAttt                                                                     <  1:1482289/91‑1 (MQ=255)
                           ggTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatat                                                                 >  1:2923864/1‑91 (MQ=255)
                           ggTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatat                                                                 >  1:322653/1‑91 (MQ=255)
                           ggTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatat                                                                 >  1:322806/1‑91 (MQ=255)
                           ggTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatat                                                                 <  1:583014/91‑1 (MQ=255)
                            gTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatata                                                                <  1:2044330/91‑1 (MQ=255)
                            gTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatata                                                                <  1:2044315/91‑1 (MQ=255)
                             tAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATaa                                                               <  1:1665176/91‑1 (MQ=255)
                              aaTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTtatataga                                                              >  1:1137761/1‑89 (MQ=255)
                              aaTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAt                                                              >  1:2481487/1‑91 (MQ=255)
                              aaTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAt                                                              >  1:2681863/1‑91 (MQ=255)
                              aaTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAt                                                              >  1:1306048/1‑91 (MQ=255)
                               aTTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAtt                                                             <  1:409788/91‑1 (MQ=255)
                                ttCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAAttt                                                            >  1:1968624/1‑91 (MQ=255)
                                  cAATACCTTGCGTAACCCGTCACTGTGAGTTCGGTTAATGATAATAAAGCCCACACCATTTAAAAAGATGTGGGTAATTTATATAAttttt                                                          >  1:3683055/1‑91 (MQ=255)
                                   aaTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTttta                                                         >  1:877096/1‑91 (MQ=255)
                                   aaTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTttta                                                         >  1:877095/1‑91 (MQ=255)
                                    aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                        >  1:2951628/1‑91 (MQ=255)
                                    aTACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttat                                                        >  1:1160969/1‑91 (MQ=255)
                                     tACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttatt                                                       <  1:1954272/91‑1 (MQ=255)
                                     tACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTGATATAATTtttatt                                                       >  1:655023/1‑91 (MQ=255)
                                       ccTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTtttattta                                                     <  1:3115246/91‑1 (MQ=255)
                                          tGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAt                                                  >  1:3094376/1‑91 (MQ=255)
                                           gCGTATCCCGTCACTGTGAGTTGGGCTAATGATAATAAATCCCACACCTTTTATAAAGATGTGGGTAATTTATATAATTTTTATTTACAtt                                                 <  1:3860288/91‑1 (MQ=255)
                                           gCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAtt                                                 >  1:2529085/1‑91 (MQ=255)
                                           gCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAtt                                                 >  1:2528782/1‑91 (MQ=255)
                                           gCGTAACCCGTCACTGTGAGTTGGGTTAAGGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAtt                                                 <  1:3325540/91‑1 (MQ=255)
                                              tAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAttttt                                              >  1:1958578/1‑91 (MQ=255)
                                              tAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACAttttt                                              <  1:2727085/91‑1 (MQ=255)
                                               aaCCCGTCACTGTGAGTTGGGTTAATGATAAAAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTCTATTTACAttatta                                             >  1:3741916/1‑91 (MQ=255)
                                                 cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAc                                           >  1:2242939/1‑91 (MQ=255)
                                                 cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAc                                           >  1:1816277/1‑91 (MQ=255)
                                                 cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAc                                           >  1:2813604/1‑91 (MQ=255)
                                                 cccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAAc                                           >  1:2957179/1‑91 (MQ=255)
                                                  ccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACt                                          <  1:2847198/91‑1 (MQ=255)
                                                  ccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACt                                          >  1:3844769/1‑91 (MQ=255)
                                                  ccGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACt                                          >  1:3669007/1‑91 (MQ=255)
                                                         tgtgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaa                                  >  1:1154739/1‑91 (MQ=255)
                                                           tgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaaa                                >  1:1471767/1‑91 (MQ=255)
                                                           tgAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATaaaaaa                                <  1:1833420/91‑1 (MQ=255)
                                                               ttGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAgcg                            <  1:1773496/91‑1 (MQ=255)
                                                                  ggTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAGTTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCt                         <  1:3491634/91‑1 (MQ=255)
                                                                  ggTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCt                         <  1:67846/91‑1 (MQ=255)
                                                                   gTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTa                        >  1:638531/1‑91 (MQ=255)
                                                                      aaTCATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGt                     <  1:3738249/91‑1 (MQ=255)
                                                                        tGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCa                   >  1:3642527/1‑91 (MQ=255)
                                                                          ataataAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATg                 >  1:124840/1‑91 (MQ=255)
                                                                          ataataAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATg                 <  1:3781882/91‑1 (MQ=255)
                                                                            aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt               >  1:3022555/1‑91 (MQ=255)
                                                                            aataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAt               >  1:1894143/1‑91 (MQ=255)
                                                                             ataaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAtt              >  1:2731999/1‑91 (MQ=255)
                                                                              taaAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGAttt             >  1:3421731/1‑91 (MQ=255)
                                                                                aaGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTAc           <  1:2665451/91‑1 (MQ=255)
                                                                                aaGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTAc           >  1:861792/1‑91 (MQ=255)
                                                                                 aGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTAcc          <  1:3930468/91‑1 (MQ=255)
                                                                                  gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa         >  1:3657653/1‑91 (MQ=255)
                                                                                  gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa         >  1:1537420/1‑91 (MQ=255)
                                                                                  gCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCa         >  1:3389918/1‑91 (MQ=255)
                                                                                   cccACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCAt        <  1:1371490/91‑1 (MQ=255)
                                                                                     cacaTCTTTTTAAAAGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTCTAAAAAAAGCGCCTATGTCATGATTTACCATca      >  1:3681896/1‑91 (MQ=255)
                                                                                     cacaTCTTTTTAAAAGGTGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTCTAAAAAAAGCGCCTATGTCATGATTTACCATca      >  1:3681719/1‑91 (MQ=255)
                                                                                     cacaCCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAAGCGCCTATGTCATGATTTACCATc       <  1:645874/91‑1 (MQ=255)
                                                                                     cacaCCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATca      >  1:2584943/1‑91 (MQ=255)
                                                                                       cacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcaca    >  1:3083288/1‑91 (MQ=255)
                                                                                       cacCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATcaca    <  1:1290887/91‑1 (MQ=255)
                                                                                        acCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACAt   >  1:3038583/1‑91 (MQ=255)
                                                                                         ccTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATa  <  1:2211374/91‑1 (MQ=255)
                                                                                         ccTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATa  <  1:2160827/91‑1 (MQ=255)
                                                                                         ccTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATa  <  1:3652558/91‑1 (MQ=255)
                                                                                         |                                                                                           
TTTACAACACCGAGCGTGCCCGTACTCGGTAATTCAATACCTTGCGTAACCCGTCACTGTGAGTTGGGTTAATGATAATAAAGCCCACACCTTTTAAAAAGATGTGGGTAATTTATATAATTTTTATTTACATTTTTAACTTATAAAAAAAGCGCCTATGTCATGATTTACCATCACATA  >  USA300TCH1516_ALE/2825498‑2825677

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: